Brevet FR2662808A1 — publié le 19911206

DESCRIPTION PROCEDE D'ANALYSE AUTOMATIQUE DE SIGNAUX PAR SEGMENTATION ET CLASSIFICATION La présente invention a pour objet un procédé d'analyse automatique de signaux par segmentation et classification. L'invention trouve une application privi Lé- giée dans le domaine biomédical, notamment dans le traitement des signaux d'origine physiologique La première application visée est l'analyse des signaux de contraction oesophagiens Mais le procédé peut très bien être appliqué à tout signal impulsionnel et bruité. La figure I annexée rappelle la forme des ondes péristaltiques à l'origine des signaux à traiter (dans l'application visée). En (a), on voit le phénomène de déglutition qui déclenche une variation de pression (P) donnant naissance d'abord à une onde de proximité (b) qui se déplace et se déforme (c, d). L'étude des signaux de contraction oesopha- giens amène à mesurer diverses grandeurs comme l'ampli- tude A du pic, la durée D de l'impulsion, les pentes maximales des fronts de montée ou de descente (d P/dt)M, l'aire sous la courbe S, etc Des analyses statistiques sont également entreprises pour obtenir le nombre et les pourcentages 3 U des différents types d'impulsions apparaissant dans l'ensemble des signaux mesurés (ondes péristaltiques, tertiaires, polyphasiques, etc). L'ensemble de ces résultats présente une signification pour le clinicien. L'entité physique à traiter (pression, PH, etc) est mesurée par diverses sondes (oesophagien- nes, intestinales, de P Hmétrie) ou divers capteurs. Les signaux traduisant cette entité physique sont ensuite amplifiés et enregistrés Avec l'avènement de l'ordinateur, le dépouillement des mesures ne se fait p Lus manuellement, mais par des procédés de traitement informatique Pour cela, le signal est échantillonné et sauvegardé dans un ordinateur muni de moyens de visualisation (écran) ou d'impression. Cette technique est décrite notamment dans "Computer-Aided Analysis of Human Esophageal Peristal- sis-(I) Technical Description and Comparison with Manual Analysis" de Donald O CASTELL et al publié dans la revue "Digestive Diseases and Sciences", vol 29, n 1, January 1984, pp 65-72 et dans "Computer Analysis of Human Esophageal Peristalsis and Lower Esophageal Sphincter Pressure (II) An Interactive System for On-Line Data Collection and Analysis" de June A CASTELL et a I, publié dans "Digestive Diseases and Sciences", vol 31, no 11, November 1986, pp 1211-1216. Dans de tels procédés de traitement, chaque impulsion est définie non pas par l'ensemble des points qui la constitue mais par quelques paramètres, comme les échantillons de début et de fin et l'échantillon de pic (ou de sommet), tous échantillons définis par leur amplitude et leur instant d'apparition. Les deux articles cites montrent essentiel- lement que le recours à l'ordinateur est fondé et fiable, dans la mesure o les résultats obtenus sont en accord très satisfaisant avec les résultats obtenus manuellement. Cette technique, bien que marquant un pro- grès, présente encore des inconvénients, notamment en ce qu'elle est semiautomatique et nécessite l'intervention du praticien en cours d'analyse pour initialiser le processus (principalement au commencement de l'analyse) La nature foncièrement non stationnaire des phénomènes étudiés complique singulièrement la tâche. Le but de l'invention est précisément de remédier à cet inconvénient en proposant un procédé d'analyse entièrement automatique permettant de déchar- ger le médecin du travail fastidieux du repérage des signaux significatifs Non seulement le procédé de l'invention est entièrement automatique, mais il ne nécessite que très peu d'informations a priori sur le signal, ce qui le rend particulièrement robuste pour des applications dans le domaine médical. Le procédé de l'invention se divise essentiellement en deux parties: segmentation et classification: pour la segmentation, on détecte les impul- sions liées au phénomène étudié à l'aide d'une méthode de croissance d'impulsions utilisant des opérateurs de fusion mettant en oeuvre des critères liés à la forme du signal; pour la classification, l'ensemble des impulsions obtenues est classé par mise en oeuvre de méthodes d'analyse de données permettant de séparer les parties significatives du bruit. Dans la première partie du procédé, relative à la détection des impulsions par segmentation, on opère par itérations: on effectue une première segmentation en fixant des seuils arbitraires pour les opérations de fusions d'impulsions On obtient alors une sur-segmentation du signal A partir du résultat obtenu, on se fixe de nouveaux seuils et l'on réitère la segmentation; on procède ainsi jusqu'à ce que le nombre de fusions ne croisse plus Le signal finalement obtenu est alors soumis à la seconde partie du traitement, à savoir la classification. De toute façon, les caractéristiques et avantages de l'invention apparaîtront mieux à La lumiè- re de la description qui va suivre Cette description porte sur des exemples de réalisation donnés à titre explicatif et nullement limitatif et elle se réfère à des dessins annexes, sur lesquels: la figure 1, déjà décrite, montre l'allure des signaux à traiter; la figure 2 montre diverses étapes dans l'opération de segmentation; la figure 3 illustre la forme d'un histogramme obtenu après segmentation; la figure 4 montre schématiquement un exemple de nuage de points finalement obtenu; la figure 5 est un organigramme général du procédé de l'invention; la figure 6 illustre une installation permettant la mise en oeuvre du procédé de l'invention. La première phase du procédé de l'invention, destinée à segmenter le signal en tronçons correspon- dants chacun à une impulsion, est illustrée par la figure 2 (lignes (a) à (e). a) Sur la ligne (a) de la figure 2 on voit une partie d'un signal impulsionnel enre- gistré On se fixe tout d'abord une amplitude d'analyse ayant une valeur A 1 égale à l'amplitude maximum des échantillons mémorisés L'échantillon correspondant Ep est pris comme échantillon unique marquant à la fois le début, le pic et la fin d'une première impulsion à constituer. b) On abaisse ensuite l'amplitude d'analyse d'une quantité déterminée et l'on recher- che les échantillons mémorisés ayant la nouvelle amplitude abaissée A 2, ligne (b). c) Pour chaque échantillon trouvé (il y en a deux, en l'occurrence soit, E Cd et E Cf), on détermine si l'échantillon est connexe à l'échantillon Ep déjà trouvé dans l'opération précédente a) et: c 1) si l'échantillon est connexe, il est pris comme nouvel échantillon de début E Cd ou de fin E Cf de la première impulsion en cours de consti- tution, c Z) si l'échantillon n'est pas connexe, il est pris comme nouvel échantillon de début, de pic et de fin d'une nouvelle impulsion à constituer. d) on réitère les opérations b) et c) en abaissant progressivement l'amplitude d'analyse A 2,, Ai-1, Ai (ligne c). Les échantillons connexes de début et de fin d'impulsion se trouvent ainsi aggrégés peu à peu aux impulsions en cours de constitution Mais de nouveaux échantillons apparaissent, comme E'p sur la ligne c); ils sont pris comme pics de nouvelles impulsions à constituer. On voit donc apparaître peu à peu des impulsions. e) A chaque itération de l'opération d) deux impulsions I 1, I 2 en cours de consti- tution sont fusionnées en une seule (cf ligne d) si: el) elles possèdent en commun un échantil- lon Ev, cet échantillon définissant une vallée entre les deux pics Epl, Ep 2, e 2) ayant calculé pour chaque impulsion I 1, I 2 la différence entre l'amplitude de l'échantillon de pic Epl, Ep 2 et l'amplitude de l'échantillon commun de vallée Ev, le rapport de ces diffé- rences est supérieur à un premier seuil T 1, e 3) la durée de l'impulsion ayant la plus petite différence d'amplitude entre pic Ep 2 et vallée Ev (dans le cas illustré sur la ligne d, c'est l'impulsion I 2) est inférieure à un deuxième seuil T 2, e 4) les deux différences d'amplitude pic-vallée sont inférieures à un troisième seuil T 3. La fusion de deux impulsions I 1, I 2 en une impulsion unique I (cf ligne e) consiste à prendre, pour échantillon de début de l'impulsion unique, l'échantillon de début Ed de la première impulsion I 2, pour échantillon de fin, l'impulsion de fin Ef de la seconde impulsion I 1 et pour échantillon de pic, l'échantillon de pic Ep de la plus haute des deux impulsions I 1. Dans cette opération de fusion, les trois seuils T 1, T 2, T 3 ont, dans cette première phase (A), trois valeurs arbitraires. f) On abaisse ainsi l'amplitude d'analyse jusqu'à atteindre une valeur minimale prédéterminée A(min) (ligne f), ce qui laisse alors, à la dernière itération correspondant à cette valeur minimale, une suite de segments S définissant chacun une impulsion I avec un échantillon de pic Ep, un échantillon de début Ed et un échantillon de fin Ef. Ayant obtenu cette première segmentation, à partir de trois seuils choisis arbitrairement, on va affiner la segmentation en modifiant ces seuils. Trois nouveaux seuils vont être définis à partir des résultats obtenus à la suite de la première segmenta- tion Pour cela, on effectue une analyse statistique des impulsions obtenues (ligne f de la figure 2) pour établir trois histogrammes, relatifs respectivement aux différences pic-vallée, aux durées et aux amplitu- des des impulsions obtenues. Ces histogrammes ont l'allure de la courbe de la figure 3 qui représente le nombre N d'impulsions ayant pour paramètre p une certaine valeur Un premier maximum M correspond aux impulsions de bruit IL est suivi d'un minimum m qui donne la valeur du nouveau seuil à prendre pour le paramètre en question (différence pic-vallée, durée, amplitude). On réitère alors les opérations d) à f) de la première phase avec les nouveaux seuils T'11, T'2, T'3 définis, ce qui permet de resegmenter le signal en refusionnant de nouvelles impulsions. On réitère les opérations d) à f) de la première phase, et ainsi de suite, jusqu'à ce que le traitement ne permette plus d'effectuer de nouvelles fusions d'impulsions et ne réduise plus le nombre d'impulsions obtenues Le signal est alors définitivement segmenté et la série d'impulsions obtenue est celle sur laquelle l'analyse peut alors être entreprise. Cette opération de segmentation définit une série d'impulsions caractérisées chacune par leur début, leur fin et leur pic Mais seul un sous-ensemble de ces impulsions est significatif, c'est-à-dire correspond à des signaux de contraction oesophagiens. Les autres impulsions correspondent au bruit dû à la respiration, aux battements cardiaques et à la toux du patient Il est donc nécessaire de poursuivre l'analyse du signal par une classification permettant de séparer les impulsions significatives du bruit. De préférence, on met en oeuvre une méthode dite des nuées dynamiques Pour cela, on définit quantitativement chaque impulsion par des paramètres comme: la différence entre l'amplitude du pic et l'amplitude de l'extrémité la plus basse, la largeur de l'impulsion, l'amplitude moyenne de l'impulsion, la largeur de l'impulsion à l'amplitude moyenne, la variance de l'amplitude, la moyenne des pentes du front de montée et du front de descente de l'impulsion, la moyenne des différences entre l'amplitu- de d'un pic et les amplitudes des deux pics voisins. On représente l'ensemble des impulsions obtenues dans un espace à autant de dimensions que de paramètres choisis, ce qui conduit à un nuage de points dans cet espace, chaque point correspondant à une impulsion. Dans le cas o les sept paramètres précédents sont retenus, on obtient un nuage de points dans un espace à sept dimensions. De préférence, on partitionne le nuage de points obtenus en deux classes, la première correspon- dant aux impulsions significatives plutôt courtes et de grande amplitude et la seconde correspondant aux impulsions plutôt longues et de faible amplitude. Selon une autre variante, on partitionne le nuage de points obtenu en trois classes, la première correspondant aux impulsions de forte amplitude, la seconde correspondant aux impulsions de très longue durée, la troisième correspondant aux impulsions inter- médiaires, et l'on partitionne ensuite la troisième classe en deux sous-classes. De préférence encore, pour adapter le poids relatif de chaque paramètre aux caractéristiques spéci- fiques de chaque enregistrement, on effectue sur le nuage de points obtenu, une analyse en composante principale, ce qui permet d'analyser le nuage de points selon ses axes principaux Ce traitement permet de séparer plus nettement les différentes classes rencon- trées. La figure 4 illustre schématiquement l'allure du nuage de points affiché sur un écran d'un calcula- teur après le traitement qui vient d'être décrit. Le nuage correspond à la distribution des impulsions selon deux paramètres P 1, P 2 Les croix sont censées représenter les impulsions non significatives et les points les impulsions significatives. La figure 5 montre l'organigramme général du procédé qui vient d'être décrit La signification des blocs représentés est la suivante: : déterminer L'amplitude maximum (A 1) et l'amplitude minimum (Amin), 12: donner à l'amplitude d'analyse A La valeur du maximum, Phase A: 14: opération de croissance des impulsions, 16: y-a-t-il connexité ? 18: s'il n'y a pas connexité, formation d'une nouvelle impulsion, : s'il y a connexité, aggrégation du nouvel échantillon à L'impulsion en cours de croissance, 22: l'amplitude d'analyse est-elle égale à l'amplitude minimum ? 24: si la réponse à la question 22 est négative, faire A=A-1, c'est-à-dire diminuer l'amplitude d'analyse, Phase B: 26: si la réponse à la question 22 est affirmative, effectuer l'analyse statis- tique des impulsions obtenues, 28: fixer de nouveaux seuils, 30: des impulsions ont-elles été fusionnées depuis la dernière itération ? Dans l'affirmative, on réitère les opérations de fusion d'impulsions avec les nouveaux seuils, dans la négative, l'analyse est terminée, 32: classification des impulsions. La figure 6, enfin, montre schématiquement un dispositif permettant de mettre en oeuvre le procédé qui vient d'être décrit Ce dispositif comprend un 1 1 calculateur 50 (qui peut être un calculateur de type PC), un clavier 52 et une imprimante 54 Le calculateur comprend une mémoire 62 contenant les données à traiter (sous forme de disquette par exemple), une mémoire 64 de type ROM contenant le programme de traitement, une unité centrale de traitement 66 et un écran 68. La figure 6 montre encore un patient 70 sur l'oesophage duquel on a posé une sonde 72 capable de mesurer la pression Un amplificateur 74 permet de délivrer un signal électrique qui est ensuite échan- tillonné puis converti en numérique et stocké dans la mémoire 62 du calculateur. Naturellement, comme le procédé qui vient d'être décrit s'effectue en différé, on peut enregis- trer préalablement le signal sur un support approprié (disquette par exemple), indépendamment de tout calculateur, puis introduire, le moment voulu, ce support dans un calculateur, pour effectuer le traitement du signal. REVENDICATIONS 1 Procédé d'analyse automatique d'un signal impulsionnel, notamment d'origine physiologique, dans lequel on mesure une entité physique, on convertit cette entité en données numériques que l'on mémorise et que l'on traite en différé pour en extraire des impulsions significatives et pour obtenir certaines caractéristiques statistiques de l'entité physique mesurée, ce procédé étant caractérisé par le fait que, le traitement destiné à extraire les impulsions significatives comprend les opérations suivantes: A) dans une première phase: on définit, dans le si- gnal, des segments correspondant chacun à une impulsion ayant un échantillon de début, marquant une extrémité du segment, un échantillon de pic marquant le sommet de l'impulsion et un échantillon de fin marquant l'au- tre extrémité du segment, cette première phase dite de segmentation comprenant les opérations suivantes: a) on fixe une amplitude d'analyse de valeur (A 1) égale à l'amplitude maximum des échantillons mémorisés, l'échantillon correspondant (Ep) étant pris comme échan- tillon de début, de pic et de fin d'une première impulsion à constituer, b) on abaisse l'amplitude d'analyse d'une quantité déterminée et l'on recherche les échantillons mémorisés ayant cette nouvelle amplitude abaissée (A 2), c) pour chaque échantillon trouvé, on détermine s'il est connexe à l'échantillon déjà trouvé dans l'opération précédente a) et: cl) si l'échantillon est connexe, il est pris comme nouvel échantillon de début (E Cd) ou de fin (E Cf) de l'impulsion en cours de constitution, c 2) si l'échantillon n'est pas connexe, il est pris comme nouvel échantillon de début, de pic et de fin d'une nouvelle impulsion à constituer, d) on réitère les opérations b) et c) en abaissant progressivement l'amplitude d'analyse (A 1, Ai-1, Ai), ce qui aggrège peu à peu des échantillons connexes de début et de fin d'impulsion aux impulsions en cours de constitution et fait apparaî- tre de nouveaux échantillons pris comme pics de nouvelles impulsions à constituer, e) à chaque itération de l'opération d), deux impulsions (I 1, I 2) en cours de constitution sont fusionnées en une seule si el) elles possèdent en commun un échantil- lon (Ev) définissant une vallée entre les deux pics (Epl, Ep 2) des deux impulsions (I 1, I 2), e 2) ayant calculé pour chaque impulsion la différence entre l'amplitude de l'échantillon de pic (Ep 1, Ep 2) et l'amplitude de l'échantillon commun de vallée (Ev), le rapport de ces différences est supérieur à un premier seuil (T 1), e 3) la durée (D 2) de l'impulsion (I 2) ayant la plus petite différence d'amplitude entre pic (Ep 2) et vallée (Ev) est inférieure à un deuxième seuil (T 2), e 4) les deux différences d'amplitude pic-vallée sont inférieures à un troisième seuil (T 3), la fusion de deux impulsions en une impul- sion unique consistant à prendre, pour échantillon de début de l'impulsion uni- que, l'échantillon de début (Ed) de la première impulsion, pour échantillon de fin, l'impulsion de fin (Ef) de la seconde impulsion et pour échantillon de pic, l'échantillon de pic (Ep) de la plus haute des deux impulsions, les trois seuils (T 1, T 2, T 3) ayant, dans cette première phase, trois valeurs arbi- traires, f) on abaisse l'amplitude d'analyse jusqu'à atteindre une valeur minimale prédétermi- née (Amin), ce qui laisse alors, à la dernière itération correspondant à cette valeur minimale, une suite de segments (S) définisssant chacun une impulsion avec une amplitude de pic, une durée et une différence picvallée, B) dans une seconde phase: g) on effectue une analyse statistique des impulsions obtenues après la première phase de segmentation et l'on détermi- ne trois histogrammes relatifs respective- ment aux différences pic-vallée, aux durées et aux amplitudes des impulsions obtenues, h) à partir des trois histogrammes obtenus, on définit trois nouveaux seuils respecti- vement pour le rapport des différences pic-vallée (T'1), pour la durée (T'2) et pour les différences pic-vallée (T'3), i) on réitère l'opération e) de la première phase avec ces nouveaux seuils (T'1, T'2, T'3), ce qui permet de resegmenter le signal en reconstituant et en refusionnant de nouvelles impulsions, j) avec les nouvelles impulsions obtenues après l'opération i), on réitère l'opération e) de la première phase, et ainsi de suite, jusqu'à ce que le traitement ne permette plus d'effectuer de nouvelles fusions d'impulsions et ne réduise plus le nombre d'impulsions obtenues, ce qui laisse un signal définitivement segmenté avec une série définitive d'impulsions sur laquelle l'analyse est alors entreprise. 2 Procédé selon la revendication 1, caracté- risé par le fait qu'il comprend, après la deuxième phase, une opération de classification des impulsions, dans laquelle: a) on définit quantitativement chaque impulsion par des paramètres comme: la différence entre l'amplitude du pic et l'amplitude de l'extrémité la plus basse, la largeur de l'impulsion, l'amplitude moyenne de l'impulsion, la largeur de l'impulsion à l'amplitude moyenne, la variance de l'amplitude, la moyenne des pentes du front de montée et du front de descente de l'impulsion, la moyenne des différences entre l'amplitu- de d'un pic et les amplitudes des deux pics voisins, b) on représente l'ensemble des impulsions obtenues dans un espace à autant de dimensions que de paramètres choisis, ce qui conduit à un nuage de points dans cet espace, chaque point correspondant à une impulsion et on analyse le nuage obtenu. 3 Procédé selon la revendication 2, caracté- risé par le fait qu'on partitionne le nuage de points obtenu en deux classes, la première correspondant aux impulsions significatives plutôt courtes et de grande amplitude et la seconde correspondant aux impul- sions plutôt longues et de faible amplitude. 4 Procédé selon la revendication 2, caracté- risé par le fait qu'on partitionne d'abord le nuage de points obtenu en trois classes, la première corres- pondant aux impulsions de forte amplitude, la seconde correspondant aux impulsions de très longue durée, la troisième correspondant aux impulsions intermédiai- res, et qu'on partitionne ensuite la troisième classe en deux sous-classes. Procédé selon la revendication 2, caracté- risé par le fait que l'on effectue, sur le nuage de points obtenu, une analyse en composante principale.