FR 3039763 A1 20170210 FR 1557582 A 20150806 1 DOMAINE GENERAL La présente invention concerne un procédé et un système pour caractériser la démarche d'une personne. ETAT DE LA TECHNIQUE Les procédés visant à caractériser la démarche d'une personne sont de deux types : les procédés de simple comptage de pas, et les procédés de détection de pas, plus complexes car visant à déterminer des instants caractéristiques du pas bien précis, tel que l'instant de pose au sol du talon. La majorité des procédés de détection de pas connus est basée sur une succession d'étapes ad hoc de traitement appliquées à un signal de mouvement mesuré sur un sujet cible. On connaît par exemple le prétraitement de Pan-Tomkins, décrit dans l'article « "A real-time QRS detection algorithm", par J. Pan et W. J Tompkins. Ce prétraitement comprend un filtrage passe bande, une dérivation, une mise au carré, une intégration et une recherche de pics appliqués au signal de mouvement mesuré sur le sujet cible. Un inconvénient important de ces procédés est qu'ils requièrent la calibration de nombreux paramètres (tailles des filtres, valeurs des seuils, critères de sélection de pics, etc...) qui sont bien souvent difficiles à estimer ou à apprendre, et sont souvent choisis de façon empirique. De plus, ces procédés reposent sur un a priori très fort sur la forme que doit avoir un pas, ce qui limite clairement leur efficacité sur des sujets cibles dont le pas est de forme non conventionnelle, par exemple des personnes atteintes d'un trouble de la démarche. RESUME DE L'INVENTION Un but de l'invention est de caractériser la démarche d'un sujet 30 sans pour autant requérir d'étape préliminaire de calibration. Cet objectif 3039763 2 peut inclure le fait de détecter des pas ou de déterminer une information caractéristique de la démarche du sujet. Un autre but de l'invention est de pouvoir mettre en oeuvre un procédé de caractérisation de démarche qui puisse fonctionner aussi bien 5 sur des sujets atteints de troubles de la démarche que de sujets sains. Il est donc proposé selon un premier aspect de l'invention un procédé de caractérisation de démarche comprenant des étapes de : - réception d'un signal de mesure de mouvement d'un sujet cible préalablement mesuré, 10 - comparaison d'au moins une portion du signal de mesure avec une pluralité de signaux de référence, la pluralité de signaux de référence comprenant au moins un signal de mouvement préalablement mesuré sur un sujet de référence différent du sujet cible, chaque comparaison du signal de mesure avec un signal de 15 référence produisant un coefficient de similarité représentatif d'une probabilité de similarité entre les signaux comparés, - sélection d'au moins un maximum local d'intérêt parmi les coefficients de similarité produits. Le procédé de caractérisation de démarche ci-dessus proposé peut 20 être complété à l'aide des caractéristiques suivantes prises seules ou bien en combinaison lorsque cela est techniquement possible. Le procédé peut comprendre une étape d'identification, dans le signal de mesure, d'au moins un instant de repère caractéristique d'un pas du sujet cible, en fonction du signal de référence à partir duquel le 25 maximum local d'intérêt sélectionné a été produit. Le procédé peut comprendre une étape de comptage du nombre de maximum local d'intérêt sélectionné correspondant à un signal de référence donné, chaque signal de référence étant associé à une information caractéristique du sujet de référence correspondant. 3039763 3 Le procédé peut comprendre une étape de comptage du nombre de maximum local d'intérêt lié chaque classes de signaux de référence, chaque classes de signaux de référence ayant été générées à partir d'informations caractéristiques des sujets de référence. 5 Au moins un signal de référence peut être comparé avec plusieurs portions du signal de mesure temporellement décalées les unes des autres, de sorte à produire une séquence de coefficients de similarité ordonnée d'après les positions temporelles des portions du signal de mesure ; un coefficient de similarité d'une séquence est alors sélectionné 10 comme maximum local en fonction du résultat d'une comparaison entre la valeur du coefficient avec la valeur du coefficient qui le précède dans la séquence, et du résultat d'une comparaison entre la valeur du coefficient avec la valeur du coefficient qui le suit dans la même séquence. Chaque signal de référence comprenant au moins un instant de 15 repère prédéterminé et caractéristique d'un pas d'un sujet de référence, au moins un signal de référence peut être comparé avec plusieurs portions du signal de mesure temporellement décalées les unes des autres, et l'instant de repère caractéristique d'un pas du sujet cible être identifié à partir de l'instant de repère prédéterminé du signal de référence à partir 20 duquel a été produit le maximum local d'intérêt sélectionné. Un maximum local peut ne pas être sélectionné lorsqu'il a une valeur inférieure à un seuil prédéterminé. Les maximums locaux peuvent être sélectionnés par ordre décroissant de valeur des coefficients de similarités. 25 Un maximum local peut ne pas être sélectionné lorsque la portion du signal mesuré, à partir duquel il a été produit, chevauche temporellement une autre portion du signal mesuré à partir de laquelle a été produit un maximum local déjà sélectionné. Le procédé peut comprendre le calcul d'un rapport entre des 30 caractéristiques dimensionnelles des deux signaux comparés. Dans ce cas, l'étape d'identification n'est pas mise en oeuvre pour un maximum local 3039763 4 produit à partir des signaux comparés dont le rapport de dimension calculé franchit un seuil prédéterminé. Les caractéristiques dimensionnelles peuvent être des amplitudes ou des variances des signaux comparés. 5 L'étape de comparaison peut comprendre une corrélation produisant un coefficient linéaire de Pearson, le coefficient de similarité dépendant de la valeur absolue du coefficient linéaire de Pearson. La pluralité de signaux de référence peut comprendre : - au moins un signal de mouvement préalablement mesuré sur un 10 sujet atteint d'un trouble de la démarche, et - au moins un signal de mouvement préalablement mesuré sur un sujet de référence sain. Le signal de mesure peut être un signal représentatif du mouvement d'un pied du sujet cible. 15 Le signal de mesure peut être ou dépendre d'au moins des signaux suivants : - un signal d'accélération du sujet cible suivant une direction normale au sol, et/ou - un signal d'accélération du sujet cible suivant une direction 20 normale à une surface supérieure d'un pied du sujet cible, et/ou - un signal de vitesse angulaire autour d'un axe de rotation d'un pied de la cible. Le procédé peut comprendre l'acquisition du signal de mesure par au moins un capteur de mouvement attaché au sujet cible. 25 Au moins un signal de référence peut être acquis au cours d'une séquence de marche comprenant une marche en ligne droite d'au moins 5 mètres d'un sujet de référence. Selon un deuxième aspect, il est en outre proposé un système de caractérisation de démarche, comprenant: 30 - une interface de réception configurée pour recevoir un signal de mesure d'accélération d'un sujet cible, et 3039763 5 - au moins un processeur configuré pour : o comparer le signal de mesure avec une pluralité de signaux de référence, la pluralité de signaux de référence comprenant au moins un signal d'accélération préalablement 5 mesuré sur un sujet de référence différent du sujet cible, chaque comparaison du signal de mesure avec un signal de référence produisant un coefficient de similarité représentatif d'une probabilité de signaux comparés, 10 o sélectionner au moins similarité entre les les un maximum local parmi coefficients de similarité produits. Le système, selon ce deuxième aspect, peut comprendre en outre au moins un capteur de mouvement adapté pour fournir le signal de mesure à l'interface de réception. 15 Ce système peut en outre comprendre une mémoire dans laquelle est mémorisée la pluralité de signaux de référence. Il est en outre proposé, selon un troisième aspect de l'invention, un procédé de diagnostic et/ou de pronostic de trouble de la démarche d'un sujet cible, comprenant une mise en oeuvre d'un procédé de 20 caractérisation de démarche selon le premier aspect de l'invention, le signal de mesure de mouvement du sujet cible ayant été préalablement mesuré. DESCRIPTION DES FIGURES 25 La figure 1 représente de façon schématique un système de caractérisation de démarche et le pied d'un sujet vu de profil. La figure 2 est un organigramme d'étapes d'une phase préliminaire à un procédé de caractérisation de démarche, selon un mode de réalisation de l'invention. 30 Les figures 3a et 3b représentent deux signaux de mesure acquis par mise en oeuvre de la phase préliminaire de la figure 2. 3039763 6 La figure 4 est un organigramme d'étapes d'un procédé de caractérisation de démarche selon un mode de réalisation de l'invention. La figure 5 illustre un signal de référence et un signal cible comparés au cours de la mise en oeuvre du procédé de détection de la 5 figure 4. Sur les figures, les éléments identiques sont désignés par des références identiques. DESCRIPTION DETAILLEE 10 En référence à la figure 1, un système de caractérisation de démarche comprend un dispositif de mesure 1, un dispositif de traitement de signaux 2 et un serveur 3. Le dispositif de mesure 1 comprend au moins un capteur de mouvement 10 et une ceinture 12. 15 La ceinture 12 est adaptée pour être attachée autour du pied d'un sujet, par exemple entre le talon et les doigts de pied comme représenté sur la figure 1. Les capteurs de mouvements 10 sont fixés à la ceinture 12 de sorte à ne pas gêner la marche du pied auquel la ceinture 12 est attachée. 20 Par exemple, il est prévu dans le dispositif de mesure 1 trois capteurs de mouvement 10, afin de mesurer des mouvements du pied suivant trois axes X, Y et Z définissant ensemble un repère de mesure attaché au dispositif 1. On utilise, par exemple, au moins un capteur de vitesse angulaire 25 autour d'un des axes précités et/ou au moins un capteur d'accélération le long d'un des axes précités. Le dispositif de mesure peut également comprendre au moins un magnétomètre. Par exemple, il est prévu trois magnétomètres, chaque magnétomètre étant configuré pour acquérir une composante de champ 30 magnétique suivant l'un des respectifs axes X, Y, Z. 3039763 7 Le dispositif de mesure 1 comprend par ailleurs des moyens de traitement de signal internes 14, et une interface de communication 16. Les moyens 14 de traitement de signal sont configurés pour élaborer au moins un signal de mouvement de sortie à partir d'au moins 5 un des signaux élémentaires fournis par chacun des différents capteurs de mouvement. L'interface de communication 16 est configurée pour communiquer avec le dispositif de traitement de signaux 2. L'interface de communication est notamment apte à transmettre au dispositif de 10 traitements de signaux 2 chaque signal de sortie produit par les moyens de traitements internes du dispositif de mesure 1. L'interface de communication 16 est de préférence sans fil pour ne pas gêner la marche du sujet lorsque le dispositif de mesure 1 est attaché à son pied. 15 Dans un mode de réalisation, le dispositif de mesure 1 est un dispositif commercialisé par la société Xsens, lequel est fixé à un pied grâce à une bande à velours-crochet, et les signaux de sortie sont échantillonnés à 100 Hz. Le dispositif 2 de traitements de signaux comprend au moins un 20 processeur 20, une première interface de communications 22, une deuxième interface de communication 24 et une mémoire interne 26. Le processeur 20 est configuré pour mettre en oeuvre un procédé de caractérisation de démarche codé sous la forme d'instructions de code de programme, sur la base de signaux reçus via les deux interfaces de 25 communications 22 et 24. L'interface de communication 22 est adaptée pour recevoir les signaux de sortie émis par le dispositif de mesure 1, par exemple au moyen du protocole Wi-Fi ou du protocole Bluetooth. L'interface de communication 24 est adaptée pour communiquer 30 avec le serveur distant 3, via un réseau public tel qu'Internet. 3039763 8 Le dispositif 2 de traitements de signaux peut par exemple se présenter sous la forme d'un terminal mobile, d'un ordinateur portable, ou d'un PC. Le serveur 3 comprend une interface de communication 34 pour 5 communiquer avec le dispositif de traitements de signaux 2, un processeur 30 et une mémoire 32 adaptée pour mémoriser un ensemble de signaux de mesures, et des données attenantes à ces signaux, sous la forme d'une base de données. On va à présent détailler les étapes d'un procédé de caractérisation 10 de démarche mis en oeuvre au moyen du système représenté sur la figure Phase préliminaire de collecte de signaux de référence En référence à la figure 2, une phase préliminaire de collecte de 15 signaux de référence est mise en oeuvre. Dans une première étape 102 de cette phase préliminaire, chaque signal de référence est mesuré sur un sujet de référence. L'acquisition 102 d'un signal est mise en oeuvre au moyen du système d'acquisition illustré sur les figures 1 et 2. Pour ce faire, on 20 attache un dispositif de mesure 1 à un pied d'un sujet de référence (voir deux dispositifs 1, l'un pour le pied gauche et l'autre pour le pied droit). Le sujet de référence est ensuite invité à marcher selon un protocole clinique prédéterminé. Au cours de la marche du sujet de référence, les capteurs de 25 mouvements 10 acquièrent des signaux représentatifs d'un paramètre de mouvement (accélération et/ou vitesse, par exemple vitesse angulaire) du pied selon les axes X, Y, Z, qui sont ensuite traités par les moyens de traitements internes du dispositif de mesure 1 pour générer au moins un signal d'accélération de référence du pied correspondant du sujet de 30 référence. 3039763 9 Les moyens de traitement internes peuvent également utiliser les mesures de champ magnétiques acquises par les magnétomètres, si présents dans les dispositifs 1 de mesure. Chacun des deux dispositifs de mesure 1 (pieds gauche et droit) 5 acquiert neuf signaux bruts (accélérations 3D, vitesses angulaires 3D et champs magnétiques 3D) dans le repère lié au capteur défini par les axes X, Y et Z, ces signaux pouvant éventuellement être recalibrés ou calibrés dans un repère fixe par rapport au sol via les moyens 14 de traitements internes des dispositifs 1 de mesure, et ce afin d'obtenir les signaux de 10 référence. De préférence, on utilise comme signaux de référence certaines des composantes pertinentes pour caractériser une démarche. Trois composantes pertinentes pour caractériser la démarche, et représentées sur la figure 1, sont: 15 - un signal d'accélération du sujet de référence mesuré suivant une direction G normale au sol, - un signal d'accélération du sujet de référence mesuré suivant l'axe Z lié au dispositif 1, lequel axe est normal à la surface supérieure d'un pied du sujet de référence, 20 - un signal de vitesse angulaire de rotation d'un pied du sujet de référence autour de l'axe Y au cours d'une marche en ligne droite du sujet de référence (cet axe Y étant perpendiculaire à l'axe Z et parallèle à un axe d'inclinaison du pied au niveau de la cheville). Chaque signal de référence est ensuite transmis au dispositif de 25 traitement 2, via l'interface 22, puis au serveur 3 via les interfaces 24 et 30. Dans une étape 104, le signal de référence mesuré est mémorisé dans la mémoire 32 du serveur 3. Un signal de référence peut être mémorisé dans la mémoire interne 26 du dispositif de traitement 2, à la 30 place de ou en complément de la mémoire 32 du serveur distant 3. 3039763 10 Dans une étape 106, chaque signal de référence mémorisé fait l'objet d'une annotation, visant à déterminer au moins un instant de repère caractéristique du pas du sujet de référence. Cette étape d'annotation peut être mise en oeuvre de façon 5 empirique par un praticien qui scrute la démarche du sujet de référence pendant l'acquisition 102 du signal de référence. Par ailleurs, peut être associé à chaque signal de référence au moins une information caractéristique du sujet de référence correspondant, préalablement déterminée. 10 Une information caractéristique de démarche peut par exemple identifier une classe particulière de sujet, parmi plusieurs classes prédéterminées. Une classe peut ainsi se rapporter à un type de démarche (normale, boitement, ...), une sévérité de l'atteinte à la démarche, une vitesse de 15 démarche, l'âge du sujet de référence (jeune, âgé, ...), etc. Par exemple, cette information peut se rapporter à une plainte exprimée par des sujets (sensoriel, moteur, équilibre, troubles cognitifs ou de l'ouïe, vertige, étourdissement ou douleur), à un déficit constaté et caractérisé (latéralisation, occurence chronique ou aiguë) du sujet 20 (proprioceptif, vestibulaire, visuel, moteur : hémiparésie, paraparésie, syndrome parkinsonien...). Cette information peut inclure le niveau anatomique de la lésion : système nerveux central ou périphérique ou dans les organes sensoriels (vestibule central, vestibule périphérique, les nerfs périphériques, la moelle épinière, le tronc, le cervelet, articulations 25 des membres inférieurs, Hémisphère cérébral: cortex, ganglions de la base ou de la matière blanche). Ces informations peuvent également porter sur l'étiologie de l'atteinte (vasculaire, inflammation, traumatisme, infection, tumeur avant ou après la chirurgie, dégénérative, toxique, métabolique, psychosomatique...). 30 Les mêmes étapes 102, 104, 106 sont mises en oeuvre de façon répétées pour plusieurs sujets de référence. 3039763 11 Au terme de la phase préliminaire, sont ainsi mémorisées dans la mémoire 32 une pluralité de signaux de référence, et chaque signal de référence est annoté avec au moins un instant caractéristique du pas du sujet de référence sur lequel le signal de référence a été mesuré, et/ou 5 au moins une information caractéristique du sujet de référence. Il est possible de mémoriser un signal dont la durée couvre au moins une partie d'un pas ou plusieurs pas du sujet de référence correspondant. Néanmoins, il est préférable de mémoriser un signal de référence présentant une durée qui couvre au moins une partie d'un seul et unique 10 pas du sujet de référence, de sorte à limiter la consommation de mémoire par le signal de référence, et d'éviter une redondance potentielle d'informations. Le ou chaque instant caractéristique mémorisé en sus du signal de référence lui-même se rapport(ent) à cet unique pas. De plus, l'utilisation d'un seul et unique pas du sujet de référence 15 ou d'une partie de ce pas permet d'augmenter la sensibilité et la spécificité de cette méthode. En outre elle permet de meilleurs résultats lorsque le pas du sujet cible est de forme non conventionnelle. L'utilisation d'au moins une partie d'un seul et unique pas du sujet de référence correspond par exemple à l'utilisation d'une partie ou de 20 plusieurs parties des phases du pas telle que la phase d'appui, le contact initial, le balancement... En pratique, cette annotation est matérialisée par la mémorisation de données permettant de localiser les instants caractéristiques dans le 25 signal de référence correspondant, ces données étant liées de façon logique avec ce signal de référence dans la mémoire 32 ou 26. Les sujets de référence auprès desquels les signaux d'accélération de référence sont mesurés présentent de préférence des profils de démarche différents. De manière générale, les sujets de référence auprès 30 desquels les signaux d'accélération de référence sont mesurés peuvent 3039763 12 être des sujets sains ou des patients (sujet atteint d'une affection médicale). Par exemple, - au moins un des signaux de référence est mesuré sur un sujet 5 atteint de trouble de la démarche (trois exemples de signaux de référence mesurés sur un patient atteint d'un problème à la hanche sont illustrés en figure 3b) ; et - au moins un des signaux de référence est mesuré sur un sujet sain, c'est-à-dire un sujet qui n'est pas atteint de troubles de la 10 démarche (tel que les trois signaux de référence représentés en figure 3a). Ainsi, l'invention porte également sur un procédé de collecte d'informations sur la démarche telles que le style de marche ou une 15 pathologie éventuelle. Ainsi, il est par exemple possible de distinguer de légères boiteries, des façons particulières de lever le pied, des jambes arquées, un chaloupement durant la marche. Même si une structure commune peut éventuellement être décelée entre les signaux de référence représentés sur les figures 3a et 3b, il convient de constater que 20 leurs amplitudes sont très différentes (sur la figure 3b, le sujet a en effet tendance à traîner des pieds, ce qui modifie la structure et l'amplitude du pas comparativement au sujet de référence de la figure 3a). Par ailleurs, les différents signaux de référence collectés ne sont pas forcément de même durée, ni associés à un nombre identique 25 d'instant caractéristiques de pas. On comprend en effet, que, d'un sujet de référence à l'autre, un pas peut être de durée différente. Les signaux de référence collectés constituent ainsi une bibliothèque de pas couvrant un large panel de démarches, que l'on va utiliser pour caractériser le pas d'un sujet qualifié dans la suite de sujet 30 « cible », ce sujet cible étant de préférence différent des sujets de référence. 3039763 13 Phase de caractérisation de démarche du sujet cible En référence à la figure 4, un procédé de caractérisation de démarche du sujet cible comprend les étapes suivantes. 5 Dans une étape 202, le dispositif de mesure 1, préalablement fixé à un pied du sujet cible, acquiert au moins un signal de mouvement du sujet cible, ce signal étant qualifié dans la suite de « signal cible ». Le procédé selon l'invention peut être mis en oeuvre sur un signal cible ayant été préalablement mesuré. 10 Un signal cible est par exemple un signal d'accélération du pied le long de l'axe X, Y et/ou Z, ou un signal de vitesse angulaire du pied autour de l'axe X, Y et/ou Z. Le signal cible et les signaux de référence sont donc mesurés sur des sujets différents. Afin d'améliorer la détection, le protocole 15 d'acquisition des signaux de référence est représentatif des conditions dans lesquels le signal cible est généré. Le signal cible est ensuite transmis au dispositif de traitement 2. Dans une étape 204, une comparaison entre le signal cible et un signal de référence est mise en oeuvre. 20 Cette comparaison 204 peut être mise en oeuvre par le processeur 20 du dispositif de traitement 1 ; dans ce cas, le dispositif de traitement 1 envoie une requête au serveur 3 visant à télécharger un signal de référence en vue de cette comparaison 204. En variante, la comparaison 204 est mise en oeuvre par le serveur 3, le dispositif de traitement servant 25 alors uniquement de relai de transmission du signal cible entre le dispositif de mesure 1 et le serveur 3. La comparaison 204 est plus précisément mise en oeuvre pour différentes portions temporelles du signal cible décalées temporellement les unes par rapport aux autres dans le signal cible. On comprend alors 30 que sont mises en oeuvre autant de comparaisons que de portions temporelles définies dans le signal cible. 3039763 14 Chaque portion temporelle du signal cible est de même durée T que le signal de référence auquel cette portion va être comparée. Les différentes portions du signal cible peuvent être décalées les unes des autres d'un pas temporel constant. 5 Les différentes portions du signal cible peuvent se chevaucher. La comparaison 204 d'une portion temporelle donnée du signal cible avec un signal de référence produit un coefficient représentatif d'une probabilité de similarité entre les deux signaux comparés, ce coefficient étant appelé dans la suite « coefficient de similarité ». 10 La comparaison 204 comprend par exemple une corrélation de la portion temporelle du signal cible avec le signal de référence ; dans ce cas, le coefficient de similarité dépend du coefficient linéaire de Pearson entre les deux signaux corrélés. Le coefficient de similarité p(x, y) entre deux signaux x et y peut 15 plus précisément être exprimé comme la valeur absolue du coefficient linéaire de Pearson, comme suit : p(x, y) = std(x)std(y) où cov(x, y) désigne la covariance ces signaux x et y, std(x) et std(y) désignent les écarts-type respectifs des signaux x et y. Dans ce cas, une valeur élevée de coefficient de similarité 20 indiquera une forte similitude entre la portion temporelle du signal cible et le signal de référence comparés. En revanche, une valeur faible de ce coefficient de similarité indiquera une faible similitude entre les deux signaux comparés. Pour N portions temporelles du signal cible, sont donc produits N 25 coefficients de similarités ayant chacun une valeur spécifique, via N comparaisons 204. Les N comparaisons opérées balayent l'intégralité de la durée du signal cible. L'ensemble des N coefficients de similarité produits à partir du signal cible et d'un même signal de référence donné est dans la suite 30 appelé séquence de coefficients. I cov(x, y)I 3039763 15 Les coefficients d'une séquence sont classés par ordre de position temporelle des N différentes portions temporelles du signal cible considérées lors de l'étape de comparaison 204. Un sens de parcours d'une séquence de coefficients correspond donc à un sens de parcours 5 temporel du signal cible. L'étape de comparaison 204 est répétée pour chaque signal de référence préalablement mémorisé lors de la phase préliminaire. On notera que, chaque signal de référence étant de durée T potentiellement distincte de celle des autres, les portions temporelles du signal cible qui 10 lui sont comparées peuvent être également redéfinies différemment pour chaque signal de référence. En croisant chaque signal de référence avec le signal cible, le nombre total de séquences de coefficients de similarité produit à l'issue de ces répétitions est égal au nombre de signaux cibles acquis 15 simultanément sur le sujet cible, multiplié par le nombre de signaux de référence mémorisés. Chaque séquence peut avoir un nombre de coefficients différent, en fonction du pas choisi pour décaler les portions temporelles les unes par rapport aux autres dans le signal cible correspondant, et en fonction 20 de la longueur (durée) du signal de référence considéré. Dans une étape 206, le processeur met en oeuvre une identification de maxima locaux dans une séquence donnée de coefficients de similarités. Dans une séquence donnée, un coefficient courant constitue un 25 maximum local si le coefficient suivant dans la séquence et le coefficient précédent dans la séquence ont des valeurs qui sont inférieures à celle du coefficient courant. La détection de maxima locaux procède typiquement par une série de comparaisons entre valeurs de coefficients adjacents dans une séquence. 3039763 16 Chaque maximum local est alors un coefficient de similarité produit à partir d'une portion du signal qui ressemble davantage au signal de référence que les portions adjacentes du signal cible. Une séquence donnée peut ainsi comprendre zéro, un ou plusieurs 5 maxima locaux. L'étape 206 est répétée pour chaque séquence de coefficients produite. Dans une étape de sélection 208, le processeur sélectionne, parmi plusieurs maxima locaux identifiés, certains maxima locaux d'intérêt. 10 Les maxima locaux d'intérêt sont sélectionnés sur la base de la valeur des coefficients de similarité : un maximum local n'est pas sélectionné comme maximum local d'intérêt si sa valeur est inférieure à un seuil prédéterminé À, de sorte à éliminer des maxima locaux ayant des valeurs de probabilité de similarité faibles. 15 Par exemple, dans le cas de coefficients de similarités définis comme la valeur absolue de coefficients linéaires de Pearson, seuls les maxima locaux ayant une valeur supérieure ou égale à 0,5 sont retenus comme maxima locaux d'intérêt. La sélection 208 peut être avantageusement être mise en oeuvre 20 par ordre décroissant de valeur des coefficients de similarité ; de la sorte, le nombre de comparaisons avec le seuil À est réduit. Un autre critère de sélection 208 d'un maximum local est le chevauchement de la portion du signal cible correspondant avec la portion temporelle du même signal cible à partir duquel a été produit un autre 25 coefficient de similarité, constituant non seulement un autre maximum local, et en plus ayant déjà été sélectionné comme maximum local d'intérêt parmi les maxima locaux de la séquence. Autrement dit, sachant qu'un premier maximum local est produit à partir d'une comparaison entre une première portion du signal cible et un 30 signal de référence, puis sélectionné comme maximum local d'intérêt, 3039763 17 - un deuxième maximum local produit à partir d'une comparaison entre le signal de référence et une deuxième portion qui chevauche la première portion n'est pas sélectionné comme maximum local d'intérêt, et 5 - un troisième maximum local produit à partir d'une comparaison entre le signal de référence et une troisième portion qui ne chevauche pas la première portion est bien sélectionné comme maximum local d'intérêt. De préférence, l'étape de sélection 208 des maxima locaux 10 d'intérêt est mise en oeuvre par ordre de valeur décroissante de valeur de coefficient ; on est de la sorte certain que les maxima locaux ayant les plus grandes valeurs dans la séquence seront sélectionnées comme maximum local d'intérêt. Dans la mesure où autant de séquences que de signaux de référence 15 peuvent être produites, l'ordre de sélection 208 peut se faire après un tri de l'ensemble des maxima locaux répertoriés dans l'ensemble des séquences. A l'issue de l'étape de sélection 208, on dispose d'au moins un maximum local d'intérêt, lequel est associé à une portion spécifique du 20 signal cible ; et quand une pluralité de maxima locaux d'intérêts sont produits, il est certain que les portions spécifiques du signal cible associées sont temporellement disjointes dans le signal cible (c'est-à-dire ne se chevauchent pas). Dans une étape d'identification 210, le processeur identifie dans le 25 signal cible au moins un instant caractéristique d'un pas du sujet cible sur la base du ou des maxima locaux d'intérêt sélectionnés. Cette étape d'identification 210 comprend le téléchargement depuis le serveur 3 vers le dispositif de traitement 2, du ou des instant(s) de repère caractéristique(s) annotés dans le signal de référence ici 30 considéré, lors de la phase préliminaire de collecte des signaux de 3039763 18 référence. En variante, l'étape d'identification est mise en oeuvre par le serveur 3. Un instant de repère caractéristique du sujet cible est identifié dans le signal cible en examinant un maximum local d'intérêt sélectionné 5 dans l'étape 208. Pour rappel, le signal de référence ayant conduit à l'obtention d'un maximum local d'intérêt donné a été annoté avec au moins un instant caractéristique au cours de la phase préliminaire et/ou avec au moins une caractéristique du sujet de référence correspondant. 10 On a par exemple représenté en figure 5 un signal de référence p comprenant deux instants caractéristiques : un instant de début de pas ta et un instant de fin de pas tb. L'instant caractéristique ta est reporté dans la portion temporelle du signal cible à la position temporelle t, + ta, où ti désigne l'instant de 15 début de la portion temporelle dans le signal cible. Similairement, l'instant caractéristique tb est reporté dans la portion temporelle du signal cible à la position temporelle t, + ta, où ti désigne l'instant de début de la portion temporelle dans le signal cible. Chaque instant de repère prédéterminé au cours de la phase 20 préliminaire est ainsi calqué dans le signal cible. Le processeur 20 procède de la même manière pour tout autre maximum local d'intérêt sélectionné. L'étape d'indentification 210 permet ainsi dans une application particulière de compter les pas du sujet cible. 25 Dans un mode de réalisation, l'identification 210 d'instant(s) de repère caractéristique(s) n'est pas forcément mise en oeuvre pour tous les maxima locaux sélectionnés au cours de l'étape de sélection 208. Pour chaque maximum local d'intérêt, le processeur calcule dans ce cas un rapport entre des caractéristiques dimensionnelles des deux signaux 30 comparés sur la base desquels le maxima a été identifié 206, puis sélectionné 208. 3039763 19 Si ce rapport de dimension franchit un seuil prédéterminé p, alors l'étape d'identification 210 d'instant caractéristique n'est pas faite dans la portion temporelle du signal cible. Par exemple, si le numérateur du rapport est la caractéristique 5 dimensionnelle se rapportant au signal cible et si le dénominateur du rapport est la caractéristique dimensionnelle se rapportant à un signal de référence, alors on vérifie si le rapport est inférieur au seuil p, et si tel est le cas, alors l'étape d'identification 210 d'instant caractéristique n'est pas faite dans la portion temporelle du signal cible. 10 Sinon, l'étape d'identification 210 selon les modalités qui précèdent est mise en oeuvre sur cette portion temporelle du signal cible. Le rapport de dimension peut être par exemple un rapport d'amplitude entre les signaux comparés ou un rapport de variance entre les signaux comparés. 15 Plusieurs rapports peuvent être calculés pour une paire de signaux comparés ; il peut alors être prévu de ne pas déclencher l'étape d'identification d'instant de repère si au moins un des rapports est supérieur à un seuil prédéfini. De préférence, le calcul des rapports de dimension est mis en 20 oeuvre après l'identification 206 des maxima locaux, très préférentiellement après la sélection des maxima locaux d'intérêt, de sorte à limiter le nombre de rapports à calculer. Les étapes 202, 204, 206, 208 et 210 peuvent être répétées pour plusieurs signaux cibles acquis simultanément via des capteurs différents 25 du dispositif de mesure 1 (par exemple des signaux d'accélération suivant des axes différents parmi ceux évoqués plus haut). D'autres traitements que l'étape d'identification 210 peuvent être mis en oeuvre une fois que des maxima locaux ont été sélectionnés dans l'étape 208 de sélection. 3039763 20 Le processeur 20 peut par exemple compter 212 le nombre de fois qu'un signal de référence a permis de générer un maximum local d'intérêt de référence. Le processeur 20 télécharge ensuite les informations 5 caractéristiques des sujets de référence (recensées au cours de la phase préliminaire) qui sont associés aux signaux de référence ayant générés au moins un des maxima locaux d'intérêt. Le processeur 20 peut par exemple calculer un pourcentage de représentation de chaque information caractéristique dans l'ensemble des 10 maxima locaux d'intérêt et donc du signal cible. Il est tout à faire possible qu'un signal de référence (ou plusieurs signaux de référence porteur d'une même information caractéristique de démarche), soit représenté de façon majoritaire dans les maxima locaux d'intérêts. 15 L'information caractéristique majoritaire peut ainsi être attribuée au sujet cible. Le processeur 20 peut aussi, pour chaque classe de signaux de référence générées à partir d'informations caractéristiques des sujets de référence, compter 212 le nombre de maximum local d'intérêt lié à cette 20 classe et donc ayant été associé au signal cible. Le processeur 20 peut ensuite par exemple calculer un pourcentage de représentation de chaque information caractéristique dans le signal cible. L'information caractéristique majoritaire peut ainsi être attribuée 25 au sujet cible. Exemple d'implémentation algorithmique du procédé de détection Le procédé de détection peut être implémenté sous la forme d'un 30 programme d'ordinateur comprenant des instructions de codes aptes à être exécutées par le processeur 20. 3039763 21 Il va être détaillé ci-après un algorithme consistant en un mode de réalisation possible de ce programme d'ordinateur. On considère dans la suite que le dispositif de mesure 1 fournit un k-uplet de de signaux discrets mesurés sur un même sujet. Ce k-uplet peut 5 être représenté par un vecteur de longueur I p I et de dimension k. On introduit les notations suivantes : - x est un vecteur représentatif d'un k-uplet de signaux cibles discrets mesurés le sujet cible. - p est un vecteur représentatif d'un k-uplet de signaux de référence 10 discrets mesurés sur un même sujet de référence. - I P I est le nombre de vecteurs p. - I x I (respectivement I p I) est la longueur du vecteur x (respectivement p), comptée en nombre d'échantillons. - x(k) (respectivement p(k)) est la kième composante du vecteur x 15 - x,(:kiest la portion temporelle de x(k) comprise entre les échantillons d'indices i et j ; on a donc xi(k.i)x1 = x(k). Les coefficients de similarités, notés r, sont calculés entre les signaux de référence et le signal cible au cours de l'étape de comparaison 204, pour toutes les portions temporelles possibles de longueur I p I , ce 20 que fait un total de 3 x IPI x lx' coefficients de similarité : Vk, Vp E P, V1 x PC 1 r(k, p, = p (p0c),;(:1?+Ipi_i) où r(k, p, i) est le coefficient de similarité issu de la comparaison entre la kième composante du signal de référence p et la kième composante du signal cible à l'instant i, et où p est un coefficient linéaire de Pearson en valeur 25 absolue. Au cours de l'étape 206, un coefficient r(k, p, 0 est désigné comme étant un maximum local s'il est plus grand que ses voisins temporaux directs. L'ensemble L des maxima locaux est ainsi l'ensemble 3039763 22 des coefficients de similarité r(k, p, 0 vérifiant cumulativement les deux conditions suivantes: r(k,p,i) > r(k,p,i - 1) r(k,p,i) > r(k,p, i + 1) L'ensemble L contient toutes les positions acceptable pour des pas du sujet cible ; chaque coefficient r(k,p,i) peut être interprété comme la 5 vraisemblance de trouver, dans le signal cible à l'instant i, un pas similaire à celui du signal de référence d'indice p. L'étape de sélection 208 est codée à la façon d'un processus glouton. A chaque itération de l'algorithme, on choisit la plus grande valeur r(k* ,p*, j*) de coefficient dans l'ensemble L : 10 - Si le pas p* positionné à l'instant i* chevauche un pas qui a déjà été détecté, on le rejette et on recommence le processus ; - Si le pas p* peut être positionné à l'instant i* sans chevaucher les détections précédentes, il est considéré comme détecté et les indices compris entre i* et i* + Ip*I - 1 sont interdits dans les 15 itérations prochaines. L'algorithme glouton de l'étape 208 de sélection prend fin lorsque tous les indices sont interdits, ou lorsque tous les maxima locaux de l'ensemble L ont été considérés, ou lorsque les valeurs des coefficients restant dans l'ensemble L sont toutes plus petites que le seuil À 20 prédéterminé. Ce seuil À agit comme un paramètre d'accélération de l'algorithme ; on évite en effet de faire tourner l'algorithme trop longtemps. L'algorithme glouton de l'étape 208 produit un ensemble S de maxima locaux d'intérêts. 25 L'étape d'identification 210 des instants caractéristiques est ensuite mise en oeuvre sur la base de cet ensemble S et du seuil p. Comme indiqué plus haut, le seuil p est utilisé pour rejeter certains instants caractéristiques, en cas d'écart dimensionnel trop grand entre une portion 3039763 23 du signal cible et une portion de signal de référence associé à un maximum local d'intérêt inclus dans l'ensemble S. Résultats 5 Les inventeurs ont pu obtenir des résultats particulièrement satisfaisants après l'établissement d'une base de données de signaux contenant 295 signaux mesurés auprès de 74 sujets dont 23 sains et 51 atteints de troubles orthopédiques divers au genou, à la hanche ou à la cheville, certains des sujets ayant réalisé le protocole plusieurs fois. 10 Le protocole clinique d'acquisition des signaux de référence, auquel ont été soumis les sujets de référence portant le dispositif de mesure 1 est le suivant : - rester immobile pendant 6 secondes, - marcher tout droit pendant 10 mètres à une vitesse de confort, 15 - faire demi-tour, - revenir au point de départ, - rester immobile pendant 2 secondes. Tous les signaux de la base de données ont été annotés manuellement avec deux instants de repère caractéristiques : 20 - un instant de début de pas ta, auquel le talon décolle du sol, - un instant de fin de pas tb, auquel les orteils sont posés sur le sol et stabilisés. Au total, ont été mémorisés 9357 pas (4674 sur le pied droit et 4683 sur le pied gauche). 25 Le fonctionnement du procédé de détection a été testé au moyen de la méthode d'évaluation qui suit. - On a choisi aléatoirement I P I pas dans la base de données qui constituent une base d'apprentissage. - On a choisi également des pas dans la base de données constituant 30 une base de test. Afin d'éviter le sur-apprentissage, tous les 3039763 24 signaux mesurés sur des sujets présents dans la base d'apprentissage sont écartés de la base de test. - On a appliqué le procédé de caractérisation de démarche à chaque signal de la base de test en utilisant les I P I signaux de la base 5 d'apprentissage comme signaux de référence. - On a ensuite calculé un score de précision et un score de rappel des instants caractéristiques de début et de fin de pas identifiés dans chaque signal de la base de test par application du procédé selon l'invention, en considérant qu'un pas est correct si la moyenne des 10 instants détectés de début et de fin de pas tombe à un instant où un pas est effectivement présent au vu des annotations collectées par les praticiens au cours de la phase préliminaire de collecte décrite plus haut. - On a moyenné ensuite les scores de précision et de rappel pour tous 15 les signaux de la base de test. Dans les résultats présentés dans la suite, ce procédé a été répété 100 fois et on s'est intéressé à la moyenne et à l'écart type des scores de précision et de rappel. Le tableau ci-dessous présente les précisions et rappels (moyenne, 20 et écart-type entre parenthèses sur 100 simulations) pour différentes nombre de I P I de signaux de référence. Valeur de I P I Précision Rappel 5 96.6 (5.04) 96.6 (1.04) 10 96.1 (3.76) 96.5 (3.42) 15 95.7 (3.73) 96.4 (3.25) 20 96.0 (2.81) 97.0 (1.76) 25 95.8 (2.67) 97.0 (1.43) 30 95.7 (2.46) 97.3 (1.09) Le seuil À a été fixé à 0,6 afin de réduire le temps de calcul. La première observation est que tous ces scores semblent assez similaires : 3039763 25 tous les scores moyens sont compris entre 95% et 97%, et cela même lorsque uniquement 5 signaux de référence sont utilisés. Tous les écarts types sont plus petits que 5%, ce qui suggère une grande stabilité de la méthode. 5 Si le nombre de signaux de référence ne semble pas affecter de façon significative les performances moyennes, il influence principalement la robustesse de la méthode. On peut observer que plus le nombre de signaux de référence augmente, plus les écarts-types diminuent, et donc plus la méthode est robuste. Cette constatation fait sens car il est intuitif 10 que plus on a de signaux de référence, plus on a de chance que toutes les pathologies et types de pas soient considérés, ce qui a tendance à uniformiser les résultats pour tous les signaux de la base de test. Lorsque I P I = 5, on observe une dissymétrie entre les écarts-type de la précision et du rappel. Dû au faible nombre de signaux de référence, on force la 15 correspondance entre des pas pathologiques de petite taille et des pas sains de plus grande taille, ce qui peut provoquer une sur-détection et qui a tendance à faire chuter la précision. Le procédé de détection selon l'invention a été comparé avec un procédé de prétraitement de Pan-Tomkins cité en introduction, configuré 20 avec les paramètres suggérés dans l'article « "Automatic step detection in the accelerometer signal,", par H. Ying, C. Silex et al. Les résultats fournis par ce procédé antérieur sont de faible précision (de l'ordre de 25.6%). En définitive, le procédé de détection mis en oeuvre conformément 25 à l'invention obtient de robustes performances, même en utilisant un petit nombre de signaux de référence. Avec seulement 20 signaux de référence, on obtient une précision de 96% et un rappel de 97%. Le procédé de détection selon l'invention peut trouver de multiples applications. Le procédé de détection peut en particulier servir de phase 30 préliminaire au diagnostic ou au pronostic de trouble de la démarche. Étant donné qu'il utilise des données ayant été préalablement mesurées, 3039763 26 il ne nécessite pas d'interaction avec le sujet cible pour être mis en oeuvre. 5 REVENDICATIONS1. Procédé de caractérisation de démarche comprenant des étapes de : - réception d'un signal de mesure de mouvement d'un sujet cible préalablement mesuré, - comparaison (204) d'au moins une portion du signal de mesure avec une pluralité de signaux de référence, la pluralité de signaux de référence comprenant au moins un signal de mouvement préalablement mesuré sur un sujet de référence différent du sujet cible, chaque comparaison du signal de mesure avec un signal de référence produisant un coefficient de similarité représentatif d'une probabilité de similarité entre les signaux comparés, - sélection (208) d'au moins un maximum local d'intérêt parmi les coefficients de similarité produits. 15 2. Procédé selon ta revendication précédente, comprenant en outre une étape d'identification (210), dans le signal de mesure, d'au moins un instant de repère caractéristique d'un pas du sujet cible, en fonction du signal de référence à partir duquel le maximum local d'intérêt sélectionné 20 a été produit. 3. Procédé selon la revendication précédente, dans lequel - chaque signal de référence comprend au moins un instant de repère prédéterminé (ta) et caractéristique d'un pas d'un sujet de 25 référence, - au moins un signal de référence est comparé avec plusieurs portions du signal de mesure temporellement décalées les unes des autres,, - l'instant de repère caractéristique d'un pas du sujet cible est identifié (210) à partir de l'instant de repère prédéterminé du 30 signal de référence à partir duquel a été produit le maximum local d'intérêt sélectionné. 3039763 28 4. Procédé selon l'une des revendications 2 à 3, comprenant le calcul d'un rapport entre des caractéristiques dimensionnelles des deux signaux comparés, et dans lequel l'étape d'identification (210) n'est pas mise en oeuvre pour un maximum local produit à partir des signaux comparés dont 5 le rapport de dimension calculé franchit un seuil prédéterminé. 5. Procédé selon l'une des revendications précédentes, comprenant en outre une étape de comptage (212) du nombre de maximum local d'intérêt sélectionné correspondant à un signal de référence donné, 10 chaque signal de référence étant associé à une information caractéristique du sujet de référence correspondant. 6. Procédé selon l'une des revendications précédentes, dans lequel : - au moins un signal de référence est comparé (204) avec plusieurs portions du signal de mesure temporellement décalées les unes des autres, de sorte à produire une séquence de coefficients de similarité ordonnée d'après les positions temporelles des portions du signal de mesure, - un coefficient de similarité d'une séquence est sélectionné comme maximum local (206) en fonction du résultat d'une comparaison entre ta valeur du coefficient avec la valeur du coefficient qui te précède dans la séquence, et du résultat d'une comparaison entre la valeur du coefficient avec ta valeur du coefficient qui te suit dans la même séquence. 7. Procédé selon l'une des revendications précédentes, dans lequel, un maximum local n'est pas sélectionné lorsque la portion du signal mesuré, à partir duquel il a été produit, chevauche temporellement une autre portion du signal mesuré à partir de laquelle a été produit un maximum local déjà sélectionné. 3039763 29 8. Procédé selon l'une des revendications précédentes, dans lequel la pluralité de signaux de référence comprend : - au moins un signal de mouvement préalablement mesuré sur un 5 sujet de référence atteint d'un trouble de la démarche, et - au moins un signal de mouvement préalablement mesuré sur un sujet de référence sain. 9. Procédé selon l'une des revendications précédentes, dans lequel le 10 signal de mesure est représentatif du mouvement d'un pied du sujet cible. 10. Système de caractérisation de démarche, comprenant: - une interface de réception configurée pour recevoir un signal de mesure d'accélération d'un sujet cible, et 15 - au moins un processeur configuré pour : o comparer le signal de mesure avec une pluralité de signaux de référence, la pluralité de signaux de référence comprenant au moins un signal d'accélération préalablement mesuré sur un sujet de référence différent du sujet cible, 20 chaque comparaison du signal de mesure avec un signal de référence produisant un coefficient de similarité représentatif d'une probabilité de similarité entre les signaux comparés, o sélectionner au moins un maximum local parmi les 25 coefficients de similarité produits.