Brevet FR2644262A1 — publié le 19900914

Procédé d'analyse d'un signal, notamment en manométrie digestive. L'invention concerne un système d'acquisition automatique de données de manométrie digestive, et un système informatisé d'interprétation de ces données. Le système a pour objectif d'être susceptible d'une mise en oeuvre dans un matériel commercialisable, compatible, conversationnel et configurable. On peut distinguer deux familles principales d'analyse en motricité digestive, sous deux types de modalités (statique ou ambulatoire): L'analyse de la motricité oesophagienne et anorectale, réalisée par acquisition dynamique (déplacement d'une sonde sur 15 à 60 cm environ) au cours d'un examen statique de 20 minutes à une heure environ; 10. L'analyse de la motricité colique et de la motricité de l'intestin grêle (duodénale, jéjunale et iléale), au cours d'un examen statique de 3 à 24 heures environ en relevés continus et sonde fixe. Les versions ambulatoires de ces analyses de motricité s'effectuent par cycles de 24 heures. Les sondes utilisées sont constituées de 3 à 8 voies d'acquisition, à partir de capteurs de pressions, in situ (jauges de contraintes) ou externes et connectés à des cathéters perfusés à orifices étagés. Les signaux provenant de chacune des voies d'acquisition sont analysés séparément, puis éventuellement en corrélation de façon à détecter les phénomènes de propagation (propulsif, stationnaire, rétropropulsif). L'acquisition des données suppose la mise au point de protocoles relativement normalisés quant aux conditions de réalisation des examens, notamment en condition de durée, et d'évènements (par exemple pour l'oesophage déglutition sèche, déglutition d'eau, douleurs, perfusion d'HCL, distension de ballonnets...) . L'objectif d'une analyse de motricité digestive est d'obtenir des signaux représentatifs du déplacement des ondes mécaniques de propagation dans le circuit digestif, et de leurs caractéristiques, ainsi que l'étude des sphincters. L'étude de la motricité digestive basée sur des mesures de manométrie, est un outil complémentaire des examens morphologiques que sont la radiologie et l'endoscopie. Outre son intérêt diagnostic, elle permet une meilleure caractérisation de la physio-pathologie et de la pharmacologie digestives. La motricité oesophagienne permet notamment de préciser les causes des mysphalgies, des douleurs thoraciques et du reflux gastro-oesophagien. La motricité ano-rectale s'applique notamment à l'étude des phénomènes de constipation, ou d'incontinence anale, et à leur rééducation. Les motricités colique et de l'intestin grêle relèvent encore essentiellement du domaine de la recherche. De façon connue, un examen classique s'effectue en reliant un enregistreur polygraphe en sortie des capteurs de pression dont les signaux ont été préamplifiés. La manipulation comprend par exemple les phases suivantes: phases d'observation des signaux reçus des capteurs, pour valider lrapela,, 10. phase d'enregistrement en continu des signaux de motricité, phase d'enregistrement successivement sur n portions d'acquisition des signaux de motricité, phase de dépouillement des résultats. Dans le cas d'un dépouillement manuel, il faut compter une demi-heure à trois quarts d'heure environ en manométrie de type oesophagienne et environ une après- midi en manométrie de type colique. On connaît également des dispositifs visant à automatiser partiellement l'opération de dépouillement des résultats, comme décrits dans l'article "systèmes automatisés de traitement des signaux de motricité digestive" C ROZE,. Gasto- entérol. Clin. Biol., 1988, 12 368-375. Ainsi, un logiciel a été conçu par la Société Sensor Médics pour le traitement des signaux de manométrie oesophagienne en sortie de polygraphe. Les signaux sont stockés dans un micro-ordinateur, et peuvent être ensuite affichés à l'écran par séquences, de façon à pouvoir effectuer diverses mesures en utilisant un système de repérage sur écran à l'aide de curseurs mobiles. Ce système connu n'a toutefois qu'un fonctionnement semi-automatique, qui est lent et délicat à manier. On connaît également un système de traitement mis au point par le laboratoire de pharmacologie de l'INRA, Toulouse (Thierry HACHET) utilisé pour la mesure de l'index moteur gastro-intestinal pendant des durées de l'ordre de 24 heures chez le chien. L'index moteur est défini comme l'aire comprise entre la courbe d'activité contractile et la ligne de base. Cet index est déterminé et stocké par tranches de 30 ou 60 minutes, sur huit pistes maximum. Les signaux recueillis en motricité digestive, et singulièrement en manométrie digestive, présentent un certain nombre de caractéristiques qui rendent délicates les opérations d'analyse automatique. Parmi ces caractéristiques, on peut noter: l'apparition d'artefacts dans le signal, le plus souvent dus à des mouvements intempestifs du patient (toux,...), la dérive de la ligne de base niveau zéro de pression du signal obtenu, lié par exemple aux effets de paroi apparaissant entre les cathéters alimentés en perfusion et la paroi de la portion de tube digestif dans laquelle ils sont introduits. Le système connu réalisé par HACHET comporte un système de correction de dérive de la ligne de base, au moyen d'un algorithme de reconnaissance basé sur la mesure de 150 valeurs comprises dans une fourchette de plus ou moins 2,3 % de la ligne de base précédente (cas d'un échantillonnage à une fréquence de 0,5 Hz). La nouvelle valeur calculée remplace la valeur précédente et ainsi de suite. Toutefois, ce principe de correction de dérive de la ligne de base ne permet pas de rendre compte, de façon suffisamment précise, des phénomènes générateurs du signal, ce qui fausse les analyses quantitatives basées sur la courbe rectifiée. L'objectif de l'invention est de fournir un procédé d'analyse qui pallie ces inconvénients de l'état de la technique. Plus précisément, un premier objectif de l'invention est de fournir un système automatisé, visant à permettre la réduction du temps de dépouillement des résultats, et l'obtention de fonctions performantes de représentation graphique, d'analyse et d'archivage des données. Un autre objectif essentiel de l'invention est de fournir un tel système qui permette de rendre compte de façon plus fidèle des phénomènes générateurs du signal analysé, notamment au moyen d'une succession adéquate d'opérations de traitement spécifiquement adaptées. Ces opérations doivent notamment permettre de discriminer, dans le signal d'origine, les caractéristiques significatives pour le traitement et l'obtention des résultats. Un objectif complémentaire est de fournir un procédé permettant une mise en oeuvre aisée et rapide, notamment au moyen d'une compression de données adéquate. Ces objectifs, ainsi que d'autres qui apparaîtont par la suite, sont atteints à l'aide d'un procédé d'analyse d'un signal du type représentable par une courbe portant des trains d'ondulations et présentant une dérive de la ligne de base, notamment pour l'analyse d'un signal de manométrie digestive, par exemple en manométrie colique, procédé caractérisé en ce qu'il comprend une première étape de correction de la dérive de la ligne de base, consistant à identifier des minima significatifs dans la courbe à trains d'ondulations, et à réaligner la dite courbe en relation avec les dits minima significatifs identifiés, et une seconde étape d'analyse quantitative de la courbe ainsi réalignée. De façon préférentielle, la dite étape d'identification des minima significatifs consiste à identifier itérativement chaque nouveau minimum significatif par explorations successives de points de la courbe suivant le dernier minimum significatif identifié, un nouveau minimum significatif étant identifié en fonction de la valeur des pentes de droites de liaison tang mem et tang s, tang mem étant défini comme la droite de liaison du dernier minimum significatif au point précédent exploré, et tang s la droite de liaison du minimum significatif au point courant exploré. Avantageusement, ceci est réalisé en mettant en oeuvre au moins un des deux processus suivants: un nouveau minimum significatif est identifié lorsque tang s et tang mem sont de signes contraires, et tang s est négative ou nulle. un nouveau minimum significatif est identifié lorsque tang s et tang mem sont de signes contraires, et que la droite de liaison des deux derniers minima significatifs tang p est strictement positive. Selon une autre caractéristique avantageuse de l'invention, la dite étape de sélection des minima significatifs est précédée d'une étape de sélection de minima candidats dans la courbe, consistant par exemple à découper la courbe en tronçons successifs et à retenir le point d'ordonnée minimale pour chacun des dits tronçons. Selon une autre caractéristique essentielle de l'invention, le nouveau minimum significatif est choisi en fonction de la durée d'un tonus musculaire. Dans ce cas, le minimum significatif est avantageusement choisi comme étant le premier minimum candidat sélectionné au-delà de la durée d'un tonus musculaire, lorsqu'aucun nouveau minimum significatif n'a été identifié pendant la dite durée d'un tonus musculaire, à partir du dernier minimum significatif identifié. On notera que cette caractéristique de l'invention est fondamentale, puisqu'elle permet d'éviter d'absorber, lors de l'analyse de la courbe, les phénomènes de tonus musculaire. Une correction maladroite de la ligne de base peut en effet conduire à fonctionner sous la forme d'un filtre passe-haut absorbant l'onde de tonus qui est d'une durée comprise entre environ 5 et environ 15 minutes, et généralement aux alentours de 10 minutes. L'absorption de ce tonus fausse ensuite les résultats et notamment le calcul de l'index moteur. Selon l'invention, l'initialisation de l'étape d'identification des minima significatifs consiste à sélectionner comme première ligne de base la droite réunissant les deux premiers minima candidats sélectionnés. Ce processus d'initialisation provoque une période transitoire de mise en route du traitement, qui converge rapidement vers un régime stable de correction de la ligne de base. Avantageusement, la dite étape de correction de la dérive de la ligne de base est précédée d'une étape de filtrage, comprenant une opération d'élimination des artefacts de la courbe et/ou une opération de filtrage à fenêtre. Selon une autre caractéristique de l'invention, la dite opération d'élimination des artefacts consiste à effectuer pour chaque point de la courbe, un calcul de la pente de la droite le liant aux points suivants de la courbe, et à initialiser le traitement d'élimination d'artefacts en cas de pente supérieure à un seuil prédéterminé, le dit traitement d'élimination des artefacts consistant à éliminer tous les points de la courbe suivant le point d'initialisation et dont la droite les liant aux dits points d'initialisation présente une pente supérieure au dit seuil prédéterminé. Préférentiellement, la dite étape de correction de la dérive de la ligne de base est suivie d'une détection des zones d'activités, dans la courbe, comprenant une opération de détection des points de la courbe réalignée d'ordonnée inférieure à un seuil prédéterminé pour la détermination des zones de début et de fin d'activité. Dans ce cas, la dite étape de détection des zones d'activité comprend avantageuse- ment une opération de détermination de l'abcisse précise de début et de fin respectivement des dites zones de début et de fin d'activité, par seuillage des points environnants au moyen d'un second seuil prédéterminé, inférieur au dit seuil de détermination des zones. La dite étape de détection des zones d'activité est préférentiellement suivie par une étape d'approximation linéaire des zones d'activité. Avantageusement, les seuils utilisés sont ajustables par l'utilisateur. D'autres caractéristiques des avantages de l'invention apparaîtront à la lecture de la description suivante d'un mode de réalisation préférentiel du procédé de l'invention, donnée à titre illustratif et non limitatif, et des dessins annexés dans lesquels: la figure 1 est un schéma synoptique illustrant un mode d'utilisation d'une chaîne d'acquisition et de traitement de signaux en manométrie digestive, selon l'invention la figure 2 est un schéma fonctionnel présentant les modules principaux d'une chaîne d'acquisition compatible avec le système d'analyse de l'invention. les figures 3A-3B, 4A-4B, 5A-5B, et 6A-6B représentent des paires de forme d'onde correspondant respectivement au signal original en sortie des capteurs de pression, aux mêmes signaux après filtrage, aux mêmes signaux après alignement, et aux mêmes signaux après identification des zones d'activité et linéarisation; les figures 7A, 7B, 7C illustrent la mise en oeuvre de l'étape de correction de la ligne de base suivant l'invention, par identification de minima significatifs au sein d'un ensemble de minima candidats; la figure 8 est un organigramme schématisant un mode de réalisation avantageux du processus d'identification des minima significatifs, pour l'étape de correction de la ligne de base, selon l'invention; les figures 9A, 9B, 9C schématisent les cas de branchement sélectif de l'organigramme de la figure 8, conduisant à l'identification de minima significatifs par des opérations d'analyse de pentes de droites de liaison entre les points explorés dans la courbe la figure 10 est une courbe schématisant une présentation des résultats, de l'index moteur et de l'activité relevée au cours de l'analyse de signaux de manométrie colique selon l'invention. Le mode de réalisation décrit ici concerne essentiellement l'acquisition et le traitement de signaux de manométrie colique. Bien entendu, les principes exposés ici sont transposables aux autres signaux de manométrie digestive, ainsi que d'une manière générale à des signaux du type constitué d'une courbe portant des trains d'ondulation et présentant une dérive de la ligne de base. Comme représenté figure 1, un mode de réalisation avantageux de la chaîne d'acquisition et de traitement comprend successivement un jeu de tubes perfusés 10 munis d'orifices étagés 11 par exemple de 5 cm en 5 cm. Chacun de ces tubes 10 est relié d'une part à un perfuseur 12 et d'autre part à un capteur de pression 13 comprenant par exemple des jauges de contraintes. Le signal de sortie du capteur de pression 13 alimente un préamplificateur 14. De façon connue, le signal de sortie du préamplificateur 14 peut être utilisé pour le traçage des courbes représentatives de l'activité de motricité digestive, sur un enregistreur polygraphe 15. L'enregistreur comporte autant de voies d'entrée 16 que de capteurs, et donc de tubes perfusés 10. Selon l'invention, les signaux pré-amplifiés obtenus en sortie du préamplificateur 14 sont acheminés vers une unité 17 d'enregistrement numérique, et d'exploitation automatique. L'unité 17 réalise d'une part l'opération d'acquisition des signaux provenant des capteurs. Cette opération d'acquisition s'effectue avantageusement de façon interactive, au moyen d'un clavier de commande 18 permettant d'une part de fournir à l'appareil les paramètres d'acquisition de traitement, et d'autre part d'assurer la saisie d'informations complémentaires à l'acquisition des signaux captés (tels que les conditions de l'analyse: apparition de douleurs, prise de repas,...). L'unité 17 comprend d'autre part des moyens de traitement des signaux enregistrés sous forme de données numériques, les résultats de traitement pouvant par exemple être sortis vers des terminaux tels qu'un écran de visualisation 19, des moyens de traçage de courbes 20, et des moyens 21 d'impression de résultats qualitatifs et/ou chiffrés. Le système d'acquisition schématisé en figure 2 consiste en un ensemble de deux cartes compatibles PC XT, insérées dans un micro ordinateur 28. Les signaux prélevés en sortie des préamplificateurs (par exemple de type BAROVAR) attaquent, après adaptation 22 en amplitude (0 à 1 V correspondant à 0 à 100 mm Hg), les entrées du convertisseur analogique - numérique 23. La logique de commande 24 associée et reliée au micro ordinateur 28 par un bus de contrôle 31 permet un échantillonage systématique des 8 voies d'entrée, à une fréquence programmable par l'utilisateur (entre 0,5 et 10 Hz), le temps de conversion du C A N 23 (80 ps) assurant une simultanéité quasi parfaite du prélèvement des échantillons sur les voies d'entrée 25. L'utilisation d'une mémoire tampon 26 (de type FIFO) de taille suffisante, contrôlée par un contrôleur de FIFO 29 en coopération avec deux tampons de bus 30, permet d'assurer les transferts de la mémoire système 26, vers la mémoire de masse 27 (disquette ou disque dur) d'un micro-ordinateur 28 sans perte d'information. Le mode de fonctionnement retenu, comporte une initialisation du nombre de voies utilisées 25 et de la fréquence d'échantillonnage, puis le lancement d'une séquence d'enregistrement. Le programme d'acquisition permet avantageusement d'interrompre une séquence d'enregistrement et d'en lancer une autre à volonté. Les communications, en phase d'acquisition, entre les cartes d'acquisition et le micro ordinateur 28 sont gérées au travers des ports et interruptions de celui-ci (Interruptions matérielles). Ces interruptions interviennent lorsque la mémoire FIFO 26 est pleine pour déclencher un transfert de la mémoire 26 vers la mémoire centrale du micro ordinateur 28, ou, au cours de ce transfert, pour répondre à toute demande de mémorisation de données émanant du convertisseur 23. Ainsi, on peut aussi bien envisager: 15. soit une procédure de vidage de la mémoire 26 lorsqu'elle est pleine, procédure acceptable quand la fréquence d'échantillonnage est faible (par exemple 10 Hz maximum) ce qui autorise un vidage complet de la mémoire 26 entre l'acquisition de deux échantillons consécutifs, soit une procédure de vidage de la mémoire systématique, tant qu'aucune interruption en écriture n'est émise par les moyens d'acquisition. Cette procédure convient par exemple pour la récupération en mode accéléré de signaux enregistrés sur bande magnétique. Le logiciel d'acquisition permet de sauvegarder le contexte de la saisie, à savoir: heure et date de début et de fin de chaque séquence d'enregistrement, fréquence d'échantillonnage, chronologie des événements extérieurs (ex: début et fin de douleurs ressenties en cours d'examen...), informations concernant le patient (identité, âge,...), et tout autre message utile pour l'exploitation ultérieure. Les données enregistrées sur disquette ou disque dur 27 sont exploitables par tout autre micro ordinateur compatible, au cours d'une étape de traitement 32 permettant d'obtenir des résultats 33 destinés à l'opérateur. Selon l'invention, on réalise avantageusement en parallèle l'acquisition des données manométriques avec des données de contexte, entrées par exemple au clavier 18 (Fig 1), telles que les indications de prise de repas, de douleurs, de prise de médicaments, de manoeuvres spécifiques (gonflage de ballonnets, avec dimension, position et durée,...) au moment de leur survenance. Ces indications seront prises en compte et restituées lors du traitement et de l'émission des résultats 33. La procédure de traitement des données stockées est illustrée par les Fig 3A,3B à 6A,6B. Le traitement du signal original (Fig 3A, 3B) comprend avantageusement les étapes suivantes: Filtrage du signal numérisé (Fig 4A, 4B) Recherche de la ligne de base (Fig5A, 5B et Fig 7A, 7B, 7C) Détection des zones d'activité et linéarisation (Fig 6A, 6B). Le signal d'origine est représenté, à titre d'exemple, sous la forme de deux portions de courbes 35, 36, correspondant à deux périodes quelconques distinctes d'acquisition. Les mêmes courbes 35, 36 se retrouvent, telles qu'elles apparaissent après les traitements successifs, en Figures 4A, SA, 6A, et 4B, 5B, 6B respectivement. Le filtrage du signal numérisé dont les effets sont schématisés en Fig 4A, 4B, se compose de deux parties: Il élimine les artefacts 37, 38, 39 en se basant sur le calcul de la dérive. Partant d'un point d'indice n on incrémente i tant que I| Y(n+i)-Y(n) || / i > seuil. Une interpolation linéaire entre les points d'indice n et n+i permet alors d'éliminer les artefacts dus à des parasites extérieurs ou à des mouvements brusques du patient au cours de l'examen. Il élimine les ondes respiratoires par utilisation d'un filtre médian de largeur de fenêtre 5 points. Ce traitement présente les propriétés intéressantes d'éliminer les bruits de type impulsionnel dont l'étendue spatiale est inférieure à la moitié de la fenêtre (cas des perturbations liées à la respiration), tout en préservant les transitions du signal en rampe et en échelon. C'est un filtre non linéaire, facilement implantable et parmi ceux nécéssitant les temps de calcul les plus courts. D'une manière générale, un filtrage à fenêtre s'exprime par m2 = (1-1/N) m1 +1/N dn(i) avec: ml = moyenne précédente (initiale) m2 = moyenne courante (nouvelle moyenne) N = taille de la fenêtre dn(i) = dernier échantillon retenu. A titre d'exemple, une loi de filtrage pourrait s'exprimer, pour une fréquence d'échantillonnage de 3 Hz, par: m2 = m1 + (mi + Xmir)C avec. N > 1820 (par exemple N=2048) pour prendre une valeur supérieure à la durée d'un tonus (voir plus loin) 5. xmn = minimum candidat (voir plus loin) La recherche de la ligne de base correspond à l'étape la plus délicate à mettre au point du fait de la complexité du signal de manométrie colique (périodes d'activité irrégulières). Pour la mise en oeuvre de la procédure, on prend en compte trois types d'ondes distinctes, en manométrie colique, à savoir: Monophasiques (1 seul pic), À Polyphasiques (au moins 2 pics), "Onde(s) + Tonus", le(s) onde(s) étant de l'un ou l'autre des types cidessus. La durée maximale d'une "onde + Tonus" est d'environ 10 min (soit 1800 points pour une fréquence d'échantillonnage de 3 Hz). Un tonus correspond à une ondulation basse fréquence (période 10 min environ) modulée par les ondes mono ou polyphasiques. La ligne de base doit respecter la forme des ondes, et donc ne pas absorber ni les ondes mono ou polyphasiques, ni les ondulations de tonus. La récupération des données (après élimination des artefacts et filtrage médian) se fait par plages de 6144 points ( soit 6Ko), mais le déplacement sur le fichier est de 4096 (soit 4Ko). La zone excédentaire de 2Ko permet un recouvrement, entre deux traitements consécutifs, d'une durée supérieure à la période d'un tonus. L'organigramme de la figure 7 présente le détail de la procédure de correction de la ligne de base. Les étapes sont les suivantes Détermination de minima locaux dans la courbe (étape 70) Initialisations des variables (étapes 71, 72, 73); Exploration itérative des minima locaux sélectionnés, de façon à identifier les minima significatifs, par tests et branchements (étapes 74 à 82), jusqu'à traitement de la totalité du fichier (étape 83). Pour chaque plage de 6144 points, les minima locaux sont déterminés de la manière suivante Dans chaque tranche de 50 points (par exemple), on extrait l'échantillon ayant la plus petite ordonnée et se trouvant dans une fenêtre autour de la moyenne récursive. Cette dernière est réactualisée après l'extraction d'un minimum local. L'abcisse et l'ordonnée du minimum local courant sont respectivement notées minxx[i] et a[minxx[i]]. On a représenté en figures 8A, 8B, 8C, l'étape de sélection des minima locaux, candidats à la sélection des minima significatifs, par des traits verticaux 85, 86, de deux tailles différentes: Les traits verticaux 85 de taille courte correspondent à des minima locaux qui n'ont pas été sélectionnés comme minima significatifs; Les traits de grande taille 86 correspondent à des minima locaux qui ont été sélectionnés comme minima significatifs. Les figures 8A, 8B, 8C correspondent à trois portions consécutives de signaux numérisés, et filtrés, sur les plages 0 à 600 points, 600 à 1200 points et 1200 à 1800 points respectivement. La détermination des minima significatifs, parmi l'ensemble des minima locaux candidats déterminés, s'effectue dans les conditions suivantes: a) La procédure précédente, qui détecte les minima locaux permet de remplir le tableau suivant: (minxx[0],minxx[1],....., minxx[hmax]) hmax étant le nombre de minima locaux extraits dans une plage de 6144 points. Le calcul des minima significatifs se fera à partir de ce tableau des minxx. La dite étape d'identification des minima significatifs consiste à identifier itérativement chaque nouveau minimum significatif 86 par explorations successives des minima candidats 85 de la courbe suivant le dernier minimum significatif identifié, un nouveau minimum significatif 86 étant identifié en fonction de la valeur des pentes de droites de liaison tang mem et tang s, tang mem étant défini comme la droite de liaison du dernier minimum significatif au point précédent exploré, et tang s la droite de liaison du minimum significatif au point courant exploré. Avantageusement, ceci est réalisé en mettant en oeuvre au moins un des deux processus suivants (étapes 76 et 77): un nouveau minimum significatif est identifié lorsque tang s et tang mem sont de signes contraires, et tang s est négative ou nulle. un nouveau minimum significatif est identifié lorsque tang s et tang mem sont de signes contraires, et que la droite de liaison des deux derniers minima significatifs tang p est strictement positive. Selon une autre caractéristique essentielle de l'invention, le nouveau minimum significatif est choisi en fonction de la durée du tonus musculaire. Dans ce cas, le minimum significatif est avantageusement choisi comme étant le premier minimum candidat 85 sélectionné au-delà de la durée d'un tonus musculaire, lorsqu'aucun nouveau minimum significatif 86 n'a été identifié pendant la dite durée d'un tonus musculaire, à partir du dernier minimum significatif 86 identifié (étape 74). b) Définition des tangentes utilisées dans le programme: (1) tangp = (a[minsj]] - a[mins[j-1]] / (mins[j] - mins[j-1]) mins[j] étant l'abcisse du minimum significatif courant (le dernier min significatif qui a été retenu). tangp est la tangente entre le minimum significatif courant et le minimum significatif du précédent. (2) tangs = (a[minxx[i]] - a[mins[jll) / (minxx[i] - mins[j]) est la tangente entre le minimum local courant et le minimum significatif courant. (3) tangmem = tangs (étape 80) Cette relation permet de mémoriser la valeur précédente de tangs. (4) Tangp et tangmem sont initialisés au début de chaque nouvelle plage, en fonction des valeurs finales de la plage précédente (étape 72). Pour estimer qu'un minimum est significatif, l'une des conditions suivantes doit être satisfaite: Si la distance entre le minxx[i], minimumlocal, et le mins[j], minimum significatif courant, est supérieure ou égale à 1800 points, alors: minxx[i] est significatif. Si tangs et tangmem ont le même signe alors: minxx[i] n'est pas significatif (84). Si tangs et tangmem n'ont pas le même signe alors: minxx[i] peut être significatif si au moins l'une des deux conditions suivantes est vérifiée (étape 77) a) tangs ≤0, pour n'importe quel signe de tangp b) tangp >0, pour n'importe quel signe de tangs. Les trois solutions correspondant à ces cas d'identification d'un minimum significatif 86i,. à partir des valeurs des tangp, tang mem et tangs et du minimum significatif précédemment sélectionné 86i sont schématisées en Fig. 9A, 9B, 9C. Dans tous les cas, on vérifie que tangs et tang mem ne sont pas de même signe. La Fig 9A illustre le cas o tang p ≤ 0 et tang s ≤0 La Fig 9B illustre le cas o tang p >0 et tang s ≤0 La Fig 9C illustre le cas o tang p >0 et tangs >0. L'étape de correction de la ligne de base est suivie d'une étape de détection des zones d'activité, qui repose sur l'initialisation de deux seuils de détection d'amplitude, ajustables par l'utilisateur, l'un assez élevé S1 (toute activité inférieure à S1 sera exclue de toute analyse ultérieure) pour une bonne efficacité de détection, l'autre plus bas S2 (par exemple S2=S1/3) pour l'estimation des abcisses du début et de la fin de chaque période d'activité et le calcul de sa durée. La courbe, bien qu'inférieure au seuil S2, peut ne pas retourner à zéro. Tout échantillon n tel que Yn totale du signal. Enfin, l'approximation linéaire des zones d'activité optionnelle associe au signal filtré et aligné une série temporelle définie par une suite de segments adjacents. Ceci permet la compression des données et l'élimination des bruits résiduels subsistant après filtrage. Calcul de la pente: pente(n) = (Yn-Yk) / (Xn-Xk) avec Xk, Yk les coordonnées de l'origine du segment k et Xn, Yn les coordonnées du point courant. Partant d'une pente pour le segment "en cours" (pente (k)) on se fixe une plage de recherche: pente (k) + o à pente (k) - o. Si la pente(n) n'appartient pas à cette plage, le point courant n'appartient pas au segment initial (indice k). Il est dit significatif (initialisation d'un nouveau segment k+l dont l'origine sera X(n-1), Y(n1) sinon on passe au point suivant (suite du segment "en cours" k). Les étapes de détection des zones d'activité et de linéarisation sont conçues pour permettre de conserver l'énergie de la courbe. Les Fig 6A, 6B fournissent l'état des courbes après la succession de traitements. Un classement des ondes est effectué qui consiste en la séparation des divers types d'ondes constitutives du signal à savoir: ondes monophasiques 62 (1 seul pic) ondes polyphasiques 63 (au moins deux pics et une durée inférieure à un seuil de durée Sd). "ondes + tonus "61 (ondes de durée supérieure ou égale au seuil Sd). Pour chaque type d'onde de nombreux paramètres peuvent être extraits dont les plus importants sont: le nombre de monotonies (nombre de segments obtenus lors de la linéarisation de l'onde) la durée Ti = Tfin - T début la valeur de crête l'amplitude moyenne Am = (SieI'i) (Si étant la surface de l'onde, Ti sa durée) 15.la fréquence des ondes élémentaires superposées aux ondes polyphasiques. Pour chaque plage du signal (typiquement 10 minutes, mais cette durée est régable par l'utilisateur) les caractéristiques suivantes sont déterminées: pourcentage d'activité pour chaque type d'ondes ( des durées des ondes d'un type donné / durée de la plage) 20. pourcentage d'activité totale pour une plage index moteur (amplitude moyenne x pourcentage d'activité totale pour une plage). Pour chaque enregistrement les paramètres retenus sont: activité totale 25. index moteur. Des représentations graphiques en fonction du temps, histogrammes et impression des résultats sous forme de tableaux peuvent être extraits, comme représenté en Fig. 10. Dans l'exemple représenté, la courbe 90 correspond par exemple à l'activité globale calculée, et la courbe 91 a la valeur de l'index moteur, sur une période de 100 minutes. L'instant d'une prise de repas est marqué. REVENDICATIONS 1. Procédé d'analyse d'un signal du type représentable par une courbe portant des trains d'ondulations et présentant une dérive de la ligne de base, notamment pour l'analyse d'un signal de manométrie digestive, par exemple en manométrie colique, procédé caractérisé en ce qu'il comprend une première étape de correction de la dérive de la ligne de base, consistant à identifier des minima significatifs (86) dans la courbe à trains d'ondulations, et à réaligner la dite courbe en relation avec les dits minima significatifs identifiés (86), et une seconde étape d'analyse quantitative de la courbe ainsi réalignée. 2. Procédé selon la revendication 1 caractérisé en ce que la dite étape d'identification des minima significatifs (86) consiste à identifier itérativement chaque nouveau minimum significatif (86i.,) par explorations successives de points de la courbe suivant le dernier minimum significatif identifié (86i), un nouveau minimum significatif (86i.1) étant identifié en fonction de la valeur des pentes de droites de liaison tang mem et tang s, tang mem étant défini comme la droite de liaison du dernier minimum significatif (86i) au point précédent exploré (85), et tang s la droite de liaison du dernier minimum significatif (86i) au point courant exploré. 3. Procédé selon la revendication 2 caractérisé en ce qu'un nouveau minimum significatif (86i+.) est identifié lorsque tang s et tang mem sont de signes contraires, et tang s est négative ou nulle. 4. Procédé selon l'une quelconque des revendications 2 et 3 caractérisé en ce qu'un nouveau minimum significatif (86i.+) est identifié lorsque tang s et tang mem sont de signes contraires, et que la droite de liaison des deux derniers minima significatifs (tang p) est strictement positive. 5. Procédé selon l'une quelconque des revendications 1 à 4 caractérisé en ce que la dite étape de sélection des minima significatifs (86) est précédée d'une étape de sélection de minima candidats (85) dans la courbe représentative du signal. 6. Procédé selon la revendication 5 caractérisé en ce que la dite sélection de minima candidats (85) consiste à découper la courbe en tronçons consécutifs, et à retenir au moins un point d'ordonnée minimale pour chacun des dits tronçons. 7. Procédé selon l'une quelconque des revendications 1 à 6 caractérisé en ce que le minimum significatif est choisi en fonction de la durée d'un tonus musculaire. 8. Procédé selon la revendication 7 caractérisé en ce que la dite durée d'un tonus musculaire correspond à une durée comprise entre en environ 5 et environ 15 minutes, préférentiellement aux environs de 10 minutes. 9. Procédé selon l'une des revendications 7 et 8 et l'une des revendications ou 6 caractérisé en ce que le nouveau minimum significatif (86i.1) est choisi comme étant le premier minimum candidat sélectionné au-delà de la durée d'un tonus musculaire, lorsqu'aucun nouveau minimum significatif (86i+i1) n'a été identifié pendant la dite durée d'un tonus musculaire, à partir du dernier minimum significatif identifié (86i). 10. Procédé selon la revendication 2 et l'une des revendications 5 ou 6 caractérisé en ce que l'initialisation de la dite étape d'identification des minima significatifs consiste à sélectionner comme première ligne de base, la droite réunissant les 2 premiers minima candidats sélectionnés. 11. Procédé selon la revendication 5 et l'une quelconque des revendications 6 à 10 caractérisé en ce que la durée d'un tronçon de sélection d'au moins un minimum condidat est comprise entre 10 secondes environ et 20 secondes environ, de préférence aux environs de 16 secondes. 12. Procédé selon l'une quelconque des revendications 1 à 11 caractérisé en ce que la dite étape de correction de la dérive de la ligne de base est précédée d'une étape de filtrage, comprenant une opération d'élimination des artefacts de la courbe et/ou une opération de filtrage à fenêtre. 13. Procédé selon la revendication 12 caractérisé en ce que la dite opération d'élimination des artefacts consiste à effectuer pour chaque point de la courbe, un calcul de la pente de la droite le liant aux points suivants de la courbe, et à initialiser le traitement d'élimination d'artefacts en cas de pente supérieure à un seuil prédéterminé, le dit traitement d'élimination des artefacts consistant à éliminer tous les points de la courbe suivant le point d'initialisation et dont la droite les liant aux dits points d'initialisation présente une pente supérieure au dit seuil prédéterminé. 14. Procédé selon l'une quelconque des revendications 1 à 13 caractérisé en ce que la dite étape de correction de la dérive de la ligne de base est suivie d'une étape de détection des zones d'activités, dans la courbe, comprenant une opération de détection des points de la courbe réalignée d'ordonnée inférieure à un seuil prédéterminé pour la détermination des zones de début et de fin d'activité. 15. Procédé selon la revendication 14 caractérisé en ce que la dite détection des zones d'activité comprend une opération de détermination de l'abcisse précise de début et de fin respectivement des dites zones de début et de fin d'activité, par seuillage des points environnants au moyen d'un second seuil prédéterminé, inférieur au dit seuil de détermination des zones. 16. Procédé selon l'une quelconque des revendications 14 et 15 caractérisé en ce que le dit seuil de détermination des zones de début et de fin d'activité est d'environ 1/32 de la dynamique de la courbe, et/ou le dit second seuil est d'environ 1/3 de la valeur du dit seuil de détermination des dites zones. 17. Procédé selon l'une quelconque des revendications 14 à 16 caractérisé en ce que la dite étape de détection des zones d'activité est suivie par une étape d'approximation linéaire des zones d'activité détectées. 18. Procédé selon l'une quelconque des revendications 13 à 15 caractérisé en ce que les dites valeurs de seuil prédéterminé sont ajustables.