CA2797886C - Procede et systeme d'analyse de l'activite respiratoire d'un patient et applications correspondantes - Google Patents

Procede et systeme d'analyse de l'activite respiratoire d'un patient et applications correspondantes Download PDF

Info

Publication number
CA2797886C
CA2797886C CA2797886A CA2797886A CA2797886C CA 2797886 C CA2797886 C CA 2797886C CA 2797886 A CA2797886 A CA 2797886A CA 2797886 A CA2797886 A CA 2797886A CA 2797886 C CA2797886 C CA 2797886C
Authority
CA
Canada
Prior art keywords
signal
phi
elementary
respiratory
phase
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CA2797886A
Other languages
English (en)
Other versions
CA2797886A1 (fr
Inventor
Patrick Hanusse
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Centre National de la Recherche Scientifique CNRS
Original Assignee
Centre National de la Recherche Scientifique CNRS
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Centre National de la Recherche Scientifique CNRS filed Critical Centre National de la Recherche Scientifique CNRS
Publication of CA2797886A1 publication Critical patent/CA2797886A1/fr
Application granted granted Critical
Publication of CA2797886C publication Critical patent/CA2797886C/fr
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/08Detecting, measuring or recording devices for evaluating the respiratory organs
    • A61B5/087Measuring breath flow
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/72Signal processing specially adapted for physiological signals or for diagnostic purposes
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/72Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7225Details of analog processing, e.g. isolation amplifier, gain or sensitivity adjustment, filtering, baseline or drift compensation
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/72Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7271Specific aspects of physiological measurement analysis
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/72Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7271Specific aspects of physiological measurement analysis
    • A61B5/7275Determining trends in physiological measurement data; Predicting development of a medical condition based on physiological measurements, e.g. determining a risk factor
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/72Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7271Specific aspects of physiological measurement analysis
    • A61B5/7278Artificial waveform generation or derivation, e.g. synthesising signals from measured signals
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61MDEVICES FOR INTRODUCING MEDIA INTO, OR ONTO, THE BODY; DEVICES FOR TRANSDUCING BODY MEDIA OR FOR TAKING MEDIA FROM THE BODY; DEVICES FOR PRODUCING OR ENDING SLEEP OR STUPOR
    • A61M16/00Devices for influencing the respiratory system of patients by gas treatment, e.g. mouth-to-mouth respiration; Tracheal tubes
    • A61M16/021Devices for influencing the respiratory system of patients by gas treatment, e.g. mouth-to-mouth respiration; Tracheal tubes operated by electrical means
    • A61M16/022Control means therefor
    • A61M16/024Control means therefor including calculation means, e.g. using a processor
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61MDEVICES FOR INTRODUCING MEDIA INTO, OR ONTO, THE BODY; DEVICES FOR TRANSDUCING BODY MEDIA OR FOR TAKING MEDIA FROM THE BODY; DEVICES FOR PRODUCING OR ENDING SLEEP OR STUPOR
    • A61M16/00Devices for influencing the respiratory system of patients by gas treatment, e.g. mouth-to-mouth respiration; Tracheal tubes
    • A61M16/06Respiratory or anaesthetic masks
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/14Fourier, Walsh or analogous domain transformations, e.g. Laplace, Hilbert, Karhunen-Loeve, transforms
    • G06F17/141Discrete Fourier transforms

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Animal Behavior & Ethology (AREA)
  • Veterinary Medicine (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Public Health (AREA)
  • General Health & Medical Sciences (AREA)
  • Mathematical Physics (AREA)
  • General Physics & Mathematics (AREA)
  • Biophysics (AREA)
  • Surgery (AREA)
  • Molecular Biology (AREA)
  • Medical Informatics (AREA)
  • Pathology (AREA)
  • Physiology (AREA)
  • Signal Processing (AREA)
  • Artificial Intelligence (AREA)
  • Psychiatry (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Pulmonology (AREA)
  • Computational Mathematics (AREA)
  • Data Mining & Analysis (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Hematology (AREA)
  • Emergency Medicine (AREA)
  • Anesthesiology (AREA)
  • General Engineering & Computer Science (AREA)
  • Algebra (AREA)
  • Software Systems (AREA)
  • Databases & Information Systems (AREA)
  • Discrete Mathematics (AREA)
  • Power Engineering (AREA)
  • Measurement Of The Respiration, Hearing Ability, Form, And Blood Characteristics Of Living Organisms (AREA)
  • Measuring Or Testing Involving Enzymes Or Micro-Organisms (AREA)

Abstract

Ce procédé d'analyse de l'activité respiratoire d'un patient comprend les étapes d'acquisition (20) d'au moins un signal d'activité respiratoire comprenant au moins un signal élémentaire correspondant à un cycle respiratoire, dont la forme générale peut être exprimée par x(t) = x 0 + x 1 cos(F(t)), où F(t) est la phase dudit signal élémentaire, et d'analyse (24, 26, 28) dudit signal d'activité respiratoire. Cette analyse (24, 26, 28) comprend les étapes d'extraction (24), à partir dudit signal d'activité respiratoire, dudit signal élémentaire, de détermination (26) d'une expression d'une équation de phase (I) dudit signal élémentaire et de détermination (26) d'une expression de la phase F(t) dudit signal élémentaire en fonction de paramètres mesurant l'anharmonicité dudit signal élémentaire et sa morphologie, à partir des fonctions pcosn et psinn définies par : (II)

Description

Procédé et système d'analyse de l'activité respiratoire d'un patient et applications correspondantes La présente invention concerne un procédé et un système d'analyse de l'activité
respiratoire d'un patient, le procédé comprenant les étapes d'acquisition d'au moins un signal d'activité respiratoire comprenant au moins un signal élémentaire correspondant à
un cycle respiratoire, dont la forme générale peut être exprimée par x(t)= x0 + x1 cos(c13(t)), où ci(t) est la phase dudit signal élémentaire, et d'analyse dudit signal d'activité respiratoire.
Elle concerne également des applications de ceux-ci à la commande d'un dispositif d'assistance respiratoire et à un dispositif de surveillance respiratoire.
Elle s'applique en particulier à la détection de troubles respiratoires ou à
la commande de dispositifs de ventilation artificielle.
Un signal d'activité respiratoire est un signal de mesure de la variation d'une grandeur liée à l'activité respiratoire du patient, telle que le débit et la pression d'air ou la concentration en oxygène et en dioxyde de carbone à l'entrée de ses voies respiratoires, ou la concentration en oxygène dans le sang. Ces grandeurs peuvent être mesurées par des appareils de mesure non invasifs, par exemple un capteur de débit ou de pression intégré à un masque placé devant la bouche du patient ou un oxymètre, ou par des capteurs internes, par exemple des capteurs de pression placés dans le circuit respiratoire du patient. De tels signaux peuvent également être déduits de signaux d'électrocardiogramme.
L'activité respiratoire est constituée d'une succession de cycles respiratoires, comprenant une phase inspiratoire et une phase expiratoire, à une fréquence appelée fréquence respiratoire. Par conséquent, les signaux d'activité respiratoire sont des signaux quasi-périodiques, comprenant une succession de signaux élémentaires, chacun de ces signaux élémentaires étant caractéristique d'un cycle respiratoire.
L'analyse de ces signaux permet de détecter des troubles ou anomalies respiratoires, tels que l'apnée du sommeil ou l'asthme. Cependant, cette analyse est généralement limitée à la détermination de la fréquence respiratoire et de sa variabilité et de l'amplitude de ces signaux, et aucune analyse de la forme d'onde de ces signaux n'est réalisée.
Or, la forme d'onde de signaux d'activité respiratoire est caractéristique de cette activité respiratoire, et leur analyse peut permettre de détecter efficacement d'éventuelles anomalies respiratoires.
On connaît de nombreuses méthodes d'analyse et de caractérisation d'un signal périodique. En particulier, l'analyse fréquentielle d'un signal permet de décrire ce signal
2 dans l'espace de Fourier. La décomposition de Fourier consiste en effet à
décomposer un signal périodique de fréquence f en une somme infinie de fonctions sinusoïdales de fréquences multiples de f, pondérées par les coefficients de Fourier. Ces coefficients de Fourier, qui constituent un codage du signal analysé, sont des paramètres caractéristiques de ce signal. En pratique, le nombre de coefficients de Fourier conservés est limité, et seuls les premiers termes de la décomposition de Fourier sont gardés. Ces termes doivent cependant être en nombre suffisant pour caractériser efficacement le signal.
Or, les signaux d'activité respiratoire sont des signaux anharmoniques, c'est-à-dire non-linéaires, et la décomposition de Fourier de tels signaux nécessite de conserver un grand nombre de coefficients, coefficients auxquels il est difficile de donner un sens physique. La décomposition de Fourier est donc inadaptée pour l'analyse de ces signaux.
L'invention a donc pour but de permettre l'analyse des formes d'onde de signaux d'activité respiratoire au moyen d'un petit nombre de paramètres porteurs d'un sens physique et constituant une signature simple et explicite de la forme de ces signaux.
A cet effet, l'invention a pour objet un procédé d'analyse du type précité, caractérisé en ce que l'analyse du signal d'activité respiratoire comprend les étapes suivantes :
- extraction, à partir dudit signal d'activité respiratoire, dudit signal élémentaire;
ck13 - détermination d'une expression d'une équation de phase F(c13) = ¨dt dudit signal élémentaire ; et - détermination d'une expression de la phase c13(t) dudit signal élémentaire en fonction de paramètres (r, rk, 430, pk) mesurant l'anharmonicité dudit signal élémentaire et sa morphologie, à partir des fonctions pcos, et psinr, définies paroc :
pcos, (t, r)= Ecos(e)¨rk k et psin, (t, r)= sin(kt)¨r .
k=1 kn k=1 k"
Le procédé selon l'invention comporte également les caractéristiques suivantes, prises séparément ou en combinaison :
- l'équation de phase est exprimée sous la forme :
dc13 1+ r2 + 2r cos() dt 1¨r2 dans laquelle r, variant dans [0,1], est un paramètre mesurant l'anharmonicité
dudit signal élémentaire ;
- le signal élémentaire est exprimé au moyen de deux paramètres r et (130, sous la forme :
3 x(t)= x, a,h sin(t,r)+ cos(t,r) où a, = x, cos(c1)0) et b, = -x sin(d) 0) , les fonctions hsin et hcos étant définies par:
\ (1+ r2)cos(t)- 2r (1- r2)sin(t) h cos :(t,r --> et h sin : (t, ->
1 + r2 - 2r cosU 1+r2 - 2r cos(t) - l'équation de phase est exprimée sous la forme :
P4) Fc1))=
Q((1)) dans laquelle P(0) et Q((I)) sont des polynômes trigonométriques ;
- l'expression de la phase 1100 est déterminée sous la forme :
t (C) = '51) akp (cb -Pk,rk )- bk P cosi (cp -Pk ,rk) k dans laquelle les fonctions psini et pcosi sont définies par:
r p cosi(tr) cos((t)¨ rk et p sin, (t, r) = sinkt)¨

k =1 k=1 Ainsi réalisé, le procédé selon l'invention permet d'analyser les signaux d'activité
respiratoire et de caractériser ces signaux au moyen d'un petit nombre de paramètres, par rapport aux procédés d'analyse de signaux périodiques selon l'état de la technique.
De plus, ces paramètres ont un sens physique, et sont caractéristiques des formes d'onde de ces signaux.
Selon un autre aspect, l'invention a également pour objet un système d'analyse de l'activité respiratoire d'un patient comprenant des moyens pour acquérir un signal d'activité respiratoire comprenant au moins un signal élémentaire correspondant à un cycle respiratoire dont la forme générale peut être exprimée par x(t)- x, + xl cos(a)(0), où
(DO est la phase dudit signal élémentaire, et des moyens pour analyser ledit signal d'activité respiratoire, caractérisé en ce que les moyens pour analyser ledit signal d'activité respiratoire comprennent :
- des moyens pour extraire, à partir dudit signal d'activité respiratoire, ledit signal élémentaire ;
dict.
- des moyens pour déterminer une expression d'une équation de phase F()=
di dudit signal élémentaire ; et - des moyens pour déterminer une expression de la phase D(t) dudit signal élémentaire en fonction de paramètres (r, rk, (1)0, pk) mesurant l'anharmonicité dudit signal élémentaire et sa morphologie, à partir des fonctions pcosn et psin, définies par:
4 rk rk p cos, (t, = )¨ et p sin, (t, r) E sin ---k k kn Le système selon l'invention comporte également les caractéristiques suivantes, prises séparément ou en combinaison :
- le système comporte des moyens pour exprimer l'équation de phase sous la forme :
el) 1 -L r 2 + 2r cos (cl)) dt 1-r2 dans laquelle r, variant dans [0,1], est un paramètre mesurant l'anharmonicité
dudit signal élémentaire ;
- le système comporte des moyens pour exprimer ledit signal élémentaire au moyen de deux paramètres r et (Do, sous la forme :
x(t)-= x0 + ai h sin(t, r)+ cos(t, r) où al =x1 cos(c1)0) et b1 = -xl sin(11)0), les fonctions hsin et hcos étant définies par:
\ + r )cos(t) - 2r , - T2 )sin(t) h cos : t, r õ et h sin : r) 1 + r2 -2r cos(t ) 1 + r 2 - 2r cos(t) - le système comporte des moyens) pour exprimer l'équation de phase sous la forme Q() dans laquelle P(1:1)) et Q(cI)) sont des polynômes trigonométriques ;
- le système comporte des moyens pour exprimer la phase (1)(t) sous la forme :
01= (D+ ak p sin, (.1) - pk ,rk)- b k p cosi(a) - pk , rk) k=1 dans laquelle les fonctions psini et pcosi sont définies par:
rk p cos (t, r) cos(kt)¨ et p sin, (t,T)-sin(kt)¨rk k =1 k=1 Selon d'autres aspects, l'invention a également pour objet un dispositif d'assistance respiratoire et un système d'analyse de l'activité respiratoire comprenant un système d'analyse de l'activité respiratoire selon l'invention.
L'invention sera mieux comprise au regard d'un exemple de réalisation de l'invention qui va maintenant être décrit en faisant référence aux figures annexées parmi lesquelles :

- la figure 1 représente de manière schématique un système d'analyse selon un mode de réalisation de l'invention ;
- la figure 2 représente un signal d'activité respiratoire ; et - la figure 3 est un schéma synoptique illustrant le procédé selon un mode de
5 réalisation de l'invention.
On a représenté sur la figure 1, un système d'acquisition et d'analyse de signaux d'activité respiratoire.
Ce système comprend des moyens 1 d'acquisition d'un signal d'activité
respiratoire, un boîtier 3 d'acquisition, relié aux moyens 1 d'acquisition, et des moyens 5 d'analyse d'un signal d'activité respiratoire, par exemple un processeur, reliés au boîtier 3 d'acquisition.
Les moyens 1 d'acquisition sont aptes à recueillir un signal d'activité
respiratoire. A
titre d'exemple, ces moyens 1 comprennent un pneumotachographe 7, placé à la sortie d'un masque respiratoire 9 couvrant le nez et la bouche d'un patient 11. Le pneumotachographe 7 est apte à mesurer de manière continue le débit d'air inspiré et expiré par le patient 11 et à transmettre au boîtier 3 d'acquisition un signal électrique analogique, caractéristique de ce débit.
Le boîtier 3 d'acquisition comprend un convertisseur analogique/numérique, apte à
convertir un signal analogique en un signal numérique, appelé par la suite signal d'activité
respiratoire, par échantillonnage et quantification du signal analogique.
Le processeur 5 est apte à analyser un signal d'activité respiratoire numérique de manière à en extraire des paramètres caractéristiques de la forme de ce signal.
On a représenté sur la figure 2 un tracé illustrant la forme d'un signal 13 d'activité
respiratoire tel que mesuré par le pneumotachographe 7. Sur ce tracé, le temps est représenté en abscisse, et le débit d'air, depuis l'extérieur vers l'appareil respiratoire du patient 11, en ordonnée. On reconnaît sur ce tracé quatre signaux élémentaires 15, correspondant chacun à un cycle respiratoire, comprenant une phase inspiratoire 16, durant laquelle le débit est positif, et une phase expiratoire 17, durant laquelle le débit est négatif.
La figure 3 est un schéma synoptique illustrant l'acquisition et l'analyse d'un signal d'activité respiratoire au moyen du système décrit en référence à la figure 1, selon un mode de réalisation de l'invention.
Dans une première étape 20 d'acquisition, le débit d'air inspiré et expiré par le patient 11 est mesuré par le pneumotachographe 7, qui transmet de manière continue au boîtier 3 d'acquisition un signal électrique analogique, dont l'amplitude instantanée est proportionnelle au débit mesuré.
6 Ce signal analogique est numérisé dans une étape 22 par le convertisseur analogique/numérique du boîtier 3 d'acquisition, par échantillonnage et quantification, et le signal ainsi numérisé, ou signal d'activité respiratoire, est transmis au processeur 5.
Ce signal 13 d'activité respiratoire est composé d'une succession de signaux élémentaires 15, correspondant chacun à un cycle respiratoire. Il n'est cependant pas rigoureusement périodique, en raison notamment de la variabilité de la fréquence respiratoire et du débit d'air inspiré et expiré pendant chacun des cycles. On définit donc la période respiratoire instantanée comme l'intervalle de temps entre le début d'un cycle respiratoire, c'est-à-dire le début de sa phase d'inspiration, et le début du cycle respiratoire précédent.
Dans une étape 24, le processeur décompose le signal 13 d'activité
respiratoire en signaux élémentaires 15. Cette décomposition peut par exemple être effectuée en détectant les instants auxquels le signal passe d'une valeur négative à une valeur positive, ce qui correspond à une transition entre la phase expiratoire d'un cycle et la phase inspiratoire du cycle suivant. Ainsi, lors de l'étape 24, le processeur détermine la fréquence respiratoire instantanée du patient et sa variabilité.
Puis, le processeur 5 analyse dans une étape 26 chacun des signaux élémentaires 15 issus de cette décomposition.
Chaque signal élémentaire x(t) est un signal anharmonique, qui peut être décrit sous la forme suivante :
x(t) = xo + X, cos(cP(t)) (1) dans laquelle toute la dépendance temporelle est contenue dans la fonction de phase (1).
Ce signal élémentaire x(t) est considéré comme un signal périodique de période T, T étant égal à l'inverse de la fréquence respiratoire instantanée.
Or, dans un signal anharmonique, la principale contribution à l'anharmonicité
provient de la brisure de symétrie de la dynamique de phase. Ainsi, toute l'information dynamique pertinente est exprimée par la dynamique de phase. Lors de l'analyse du signal x(t), il convient donc d'étudier cette phase 43(t), et en particulier la dynamique de phase exprimée par la fonction F, dérivée de la fonction (13 par rapport au temps t:
F(43) = cl(1dt 3 (2) Ainsi, la morphologie du signal x(t) est complètement déterminée par la connaissance de F.
L'étape d'analyse 26 du procédé selon l'invention consiste donc à décrire cette fonction F au moyen d'un très petit nombre de paramètres. On entendra par petit nombre
7 de paramètres un nombre de paramètres réduit par rapport au nombre de paramètres nécessaire à la décomposition de la même fonction, au moyen des séries de Fourier, avec un niveau de précision équivalent.
Cette étape d'analyse 26 comprend ainsi une première étape consistant à
exprimer la phase cl), et en particulier la fonction F, dérivée de (I) par rapport au temps.
Dans le cas le plus simple, et pour un signal de période 27c, la dynamique de phase peut être écrite sous la forme :
F(1)) = doI) = 1+r2 + 2r cos(d)) (3) dt 1- r2 appelée équation de phase.
La fonction F présente dans ce cas une symétrie de réflexion par rapport à
l'axe GD=0. Cette expression de la dynamique de phase ne contient qu'un seul paramètre, r, qui varie dans l'intervalle [0,1]. La limite r=0 correspond à un signal harmonique, la limite r=1 à un signal infiniment anharmonique.
Le signal x(t), qui peut être écrit sous la forme :
x(t)-x0 + x1 cos(d)(t, r) - (Do ) (4) où (1)0 est une origine de phase, est décomposé et réécrit sous une forme faisant intervenir les paramètres r et (Do :
x(t) = x0 + aih sin(t, r)+ b cos(t, r) (5) avec ai = x1 cos(d0) et b, = 1 sin(d),:,), et dans laquelle on a défini les fonctions hcos et hsin suivantes :
(1+ r2)cos(t)- 2r h cos : r) --->(6) 1+r - 2r cos(t) (1- r2 )sin(t) hsin: r)-> (7) 1+ r 2 - 2r cos(t) Ainsi, la décomposition du signal x(t) fait intervenir seulement deux paramètres, r et 1.
r, appelé paramètre d'anharmonicité, mesure le degré d'anharmonicité du signal, la limite r=0 correspond à un signal harmonique, la limite r=1 à un signal infiniment anharmonique. Par ailleurs, le paramètre (1)2, qui définit la composition du signal dans les deux fonctions hcos et hsin, est un paramètre de morphologie, qui correspond à
l'angle de symétrie de réflexion de la dynamique de phase.
Dans le cas général, c'est-à-dire pour un signal périodique quelconque, l'équation de phase peut s'écrire sous la forme :
8 F()= Pr, (4.) Q ni (43) (8) dans laquelle Pr, et Qm sont des polynômes trigonométriques de degrés respectifs n et m.
La forme générale d'un polynôme trigonométrique de degré n est:
13, (cl)) = a 0 En ak cos,k(13, + bk sin(kcP) (9) k=1 L'analyse du signal x(t) comprend alors la détermination d'une expression de c13 faisant intervenir un petit nombre de paramètres, ce qui permet de déterminer une expression du signal x(t) en fonction de ces paramètres.
Avantageusement, l'équation de phase (2) peut être réécrite sous la forme :
1 dt Qm (cD) F(cI)) P n (cI)) (10) La factorisation du polynôme F',.,() permet de transformer F(1,) en une somme de termes simples, ce qui permet de réécrire l'équation de phase sous la forme :
d t ak cos(c13 ¨ pk)+ bk sin(c13 + pk) (11) cel) k=1 (1+ rk2 ¨ 2rk cos(c1) + pk )) dans laquelle les paramètres rk, compris entre 0 et 1, mesurent l'anharmonicité du signal x(t), et les paramètres pk caractérisent sa morphologie.
La période T du signal peut être déterminée en intégrant cette équation par rapport à (13, entre 0 et 2t:
<I>=27c ce ( r ak T= F(E13.) = 2A- a + Ek 1 ¨k __________ (12) ,D=o A partir de ce résultat, et des contraintes selon lesquelles la période est égale à 27c et le signal est harmonique lorsque les coefficients rk sont tous nuls, l'équation de phase peut être exprimée ainsi :
d t =1+EDk(cP¨pk) (13) k=1 Où la fonction Dk est définie par:
rk(ak cos(*) + bk sin()¨ ak) Dk : C13 (14) (1+ ¨ 2rk cos(<13)) Et vérifie :
<I>=2z f D k(C0C1c13 = 0 (15) d¨o
9 La définition des fonctions des fonctions polycos et polysin, notées pcosn et psinn, qui s'expriment par:
rk pcosn(t,r)=cos(kt)¨ (16) k=1 kn psinn(t,r), sin(kt)4- (17) k=1 kn et possèdent entre autre les propriétés suivantes :
r(cos(t)¨r) pcoso(t, r) = (18) 1+r 2 ¨ 2r cos(t) r sin(t) psino(t,r)= (19) 1+r 2 ¨ 2r cos(t) pcosl (t, = --1In(1+ r2 ¨ 2r cos(t)) (20) psini(t,r), tan'r rsin(t) (21) r cos(t)) permet de réécrire l'équation de phase sous la forme :
¨dt=1 + En akpcoso (c13¨ pk , rk )+ bk psin 0 (.13 ¨ pk , rk ) (22) k=1 La résolution de cette équation permet d'accéder à une expression analytique de t(<13), qui s'exprime par:
t(c1)), + En akpsinl(d)¨ pk,rk )¨ bkpcosi(c1D¨ pk, rk) (23) k=1 Le temps t est donc exprimé en fonction de la phase <13, et de manière duale la phase <D est exprimée en fonction du temps t, à l'aide de paramètres indépendants clairement définis, qui mesurent l'anharmonicité (paramètres r ou rk), et la morphologie (paramètres c1)0 ou pk).
Ainsi, lors de l'étape 26 d'analyse, le processeur 5 code chaque signal élémentaire x(t) au moyen d'un petit nombre de paramètres. Selon un mode de réalisation, chaque signal élémentaire x(t) est décrite de manière quasi-exacte par une amplitude, une harmonicité r et une morphologie .1)0. Selon un autre mode de réalisation, chaque signal élémentaire x(t) est décrit de manière encore plus précise par deux couples de paramètres (ri, Pi) et (r2, P2), complétés de leurs poids respectifs.
Chacun des signaux élémentaires, donc chaque cycle respiratoire, est donc caractérisé par un nombre restreint de paramètres, porteurs d'un sens physique car représentatifs de la non-linéarité et de la morphologie de ce signal.

Puis, dans une étape 28, le processeur compare les valeurs des paramètres déterminés pour chaque signal élémentaire, c'est-à-dire la fréquence respiratoire instantanée, l'amplitude et les paramètres de morphologie et d'harmonicité du signal, à
des valeurs tabulées de paramètres préalablement enregistrées, de manière à
détecter 5 d'éventuelles anomalies respiratoires.
Le procédé selon l'invention permet ainsi d'analyser l'activité respiratoire, et d'extraire d'un signal d'activité respiratoire un nombre restreint de paramètres, permettant une représentation compacte et pertinente de la forme d'onde de ce signal, et la détection précise d'éventuelles anomalies respiratoires.
10 Les étapes 20, 22, 24, 26 et 28 peuvent être réalisées au fur et à
mesure de l'acquisition du signal d'activité respiratoire, de manière à surveiller en continu l'activité
respiratoire du patient.
Selon un autre mode de réalisation, le système et le procédé selon l'invention peuvent être mis en oeuvre dans un respirateur artificiel.
Un respirateur artificiel est un dispositif visant à pallier l'insuffisance respiratoire d'un patient, en insufflant de l'air vers les poumons du patient par l'intermédiaire d'un masque, au moyen d'un ventilateur. Pour que cette ventilation soit efficace, il est nécessaire que les efforts inspiratoires du patient et le déclenchement du ventilateur soient synchronisés, c'est-à-dire que l'air soit insufflé dans les poumons par le ventilateur seulement pendant la phase inspiratoire du cycle respiratoire du patient.
La mesure du débit d'air inspiré et expiré par le patient permet de détecter le début et la fin de cette phase inspiratoire, et ainsi de synchroniser le fonctionnement du ventilateur sur le cycle respiratoire du patient.
Cette détection est généralement réalisée en fixant deux seuils de débit ou de pression, le ventilateur étant actionné dès que le débit ou la pression d'air inspiré par le patient dépasse un premier seuil, et arrêté dès que ce débit ou cette pression devient inférieur à un deuxième seuil.
Cette méthode de détection peut s'avérer inefficace, en raison notamment de la variabilité des cycles respiratoires, conduisant à une inadéquation entre les seuils de débit fixés et ces cycles respiratoires. Cette inadéquation peut entraîner une mauvaise synchronisation entre le ventilateur et le patient, ou même un non déclenchement du ventilateur.
Ainsi, selon un mode de réalisation de l'invention, le procédé d'analyse de l'activité
respiratoire selon l'invention est mis en oeuvre pour déterminer, préalablement à
l'utilisation du respirateur, les valeurs de seuils des paramètres respiratoires les mieux adaptés pour un fonctionnement optimal du respirateur. Cette analyse est également
11 réalisée de manière continue pendant le fonctionnement du respirateur, de manière à
adapter ces seuils à la respiration du patient.
Selon un autre mode de réalisation, le système et le procédé selon l'invention peuvent être mis en oeuvre en imagerie tomodensitométrique (TMD), l'analyse des signaux d'activités respiratoires permettant de corriger, notamment sur les images des poumons, les artéfacts dus aux mouvements respiratoires.
Selon un autre mode de réalisation, le procédé de mesure de l'activité
respiratoire d'un patient comprend les étapes suivantes :
- placer une interface sur les voies respiratoires du patient, ladite interface comprenant un pneumotachographe apte à mesurer de manière continue le débit d'air inspiré et expiré par le patient ;
- acquérir au moins un signal d'activité respiratoire du patient comprenant au moins un signal élémentaire correspondant à un cycle respiratoire, dont la forme générale peut être exprimée par x(t) = xo +x1 cos(cI)(t)), où t est le temps, IO est la phase dudit signal élémentaire, xo et x, des coefficients, - analyser ledit signal d'activité respiratoire, où l'analyse dudit signal d'activité
respiratoire comprend les étapes suivantes :
- extraction, à partir dudit signal d'activité respiratoire, dudit signal élémentaire ;
\ de) - détermination d'une expression d'une équation de phaseF(c13) = ¨t dudit d signal élémentaire ; et - détermination d'une expression de la phase MO dudit signal élémentaire en fonction de paramètres mesurant l'anharmonicité dudit signal élémentaire et sa morphologie, à partir des fonctions pcos, et psin,, définies par:
p cos n(t , r) .--Icos(kt) ¨r et p sin, (t, r) sin (k) rk , k" k"
k=1 k =1 où r est un paramètre d'anharmonicité mesurant le degré d'anharmonicité du signal, qui varie dans l'intervalle [0,1], la limite r=0 correspondant à un signal harmonique, la limite r=1 à un signal infiniment anharmonique, et où k et n sont des entiers.
Selon un autre mode de réalisation, le système de mesure de l'activité
respiratoire d'un patient comprend une interface pour recouvrir les voies respiratoires du patient,comprenant un pneumotachographe apte à mesurer de manière continue le débit 1 1 a d'air inspiré et expiré par le patient ; des moyens pour acquérir un signal d'activité
respiratoire du patient comprenant au moins un signal élémentaire correspondant à un cycle respiratoire dont la forme générale peut être exprimée par x(t)=. xo +
x1 cos(1)(0), où
t est le temps, (1)(t) est la phase dudit signal élémentaire, xo et xi des coefficients ; des moyens pour analyser ledit signal d'activité respiratoire, où les moyens pour analyser ledit signal d'activité respiratoire comprennent : des moyens pour extraire, à
partir dudit signal d'activité respiratoire, ledit signal élémentaire ; des moyens pour déterminer une \ ch) expression d'une équation de phase F((1)r ________________________________ dudit signal élémentaire ; et des moyens dt pour déterminer une expression de la phase (DO dudit signal élémentaire en fonction de paramètres (r, rk, *0, pk) mesurant l'anharmonicité dudit signal élémentaire et sa morphologie, à partir des fonctions pcos, et psin, définies par:
p cos, (t, r) cos(kt) ' et p sin, (t,r)= sin(kt)¨rk , k=1 kn k =1 kn où r est un paramètre d'anharmonicité mesurant le degré d'anharmonicité du signal, qui varie dans l'intervalle [0,1], la limite r=0 correspondant à un signal harmonique, la limite r=1 à un signal infiniment anharmonique, et où k et n sont des entiers.
Il devra toutefois être compris que les exemples de réalisation présentés ci-dessus ne sont pas limitatifs.
Notamment, le signal d'activité respiratoire n'est pas nécessairement relatif au débit d'air inspiré ou expiré.

Claims (14)

REVENDICATIONS
1.- Procédé de mesure de l'activité respiratoire d'un patient, comprenant :
- placer une interface sur les voies respiratoires du patient, ladite interface comprenant un pneumotachographe apte à mesurer de manière continue le débit d'air inspiré et expiré par le patient ;
- acquérir au moins un signal d'activité respiratoire du patient comprenant au moins un signal élémentaire correspondant à un cycle respiratoire, dont la forme générale peut être exprimée par x(t)=x0+ x1, cos(.PHI.(t)), où t est le temps, (.PHI.(t) est la phase dudit signal élémentaire, x0 et x1 des coefficients, - analyser ledit signal d'activité respiratoire, où l'analyse dudit signal d'activité
respiratoire comprend les étapes suivantes :
- extraction, à partir dudit signal d'activité respiratoire, dudit signal élémentaire ;
- détermination d'une expression d'une équation de phase F(.PHI.)=~ dudit signal élémentaire ; et - détermination d'une expression de la phase .PHI.(t) dudit signal élémentaire en fonction de paramètres mesurant l'anharmonicité dudit signal élémentaire et sa morphologie, à partir des fonctions pcosn et psinn définies par :
où r est un paramètre d'anharmonicité mesurant le degré d'anharmonicité du signal, qui varie dans l'intervalle [0,1], la limite r=0 correspondant à un signal harmonique, la limite r=1 à un signal infiniment anharmonique, et où k et n sont des entiers.
2.- Procédé selon la revendication 1, dans lequel l'équation de phase est exprimée sous la forme : dans laquelle r, variant dans [0,1], est un paramètre mesurant l'anharmonicité
dudit signal élémentaire.
3.- Procédé selon la revendication 2 dans lequel le signal élémentaire est exprimé
au moyen de deux paramètres r et .PHI.0, sous la forme :
x(t)= x0+ a1h sin(t,r) + b1h cos(t,r) où a1 = x1 cos(.PHI.0 ) et b1 = ¨x1 sin(.PHI.0), les fonctions hsin et hcos étant définies par :
4.- Procédé selon la revendication 1, dans lequel l'équation de phase est exprimée sous la forme :
dans laquelle P(.PHI.) et Q(.PHI.) sont des polynômes trigonométriques.
5.- Procédé selon la revendication 4, dans lequel l'expression de la phase .PHI.(t) est déterminée sous la forme :
dans laquelle les fonctions psin1 et pcos1 sont définies par :
6.- Procédé selon l'une quelconque des revendications 1 à 5, dans lequel ladite interface comprend un masque respiratoire et qui est placé sur le nez et la bouche du patient.
7.- Système de mesure de l'activité respiratoire d'un patient comprenant :
- une interface pour recouvrir les voies respiratoires du patient, comprenant un pneumotachographe apte à mesurer de manière continue le débit d'air inspiré et expiré
par le patient ;
- des moyens pour acquérir un signal d'activité respiratoire du patient comprenant au moins un signal élémentaire correspondant à un cycle respiratoire dont la forme générale peut être exprimée par x(t) = x0 + x1 cos(.PHI.(t)), où t est le temps, .PHI.(t) est la phase dudit signal élémentaire, x0 et x1 des coefficients ;

- des moyens pour analyser ledit signal d'activité respiratoire, où les moyens pour analyser ledit signal d'activité respiratoire comprennent :
- des moyens pour extraire, à partir dudit signal d'activité respiratoire, ledit signal élémentaire ;
- des moyens pour déterminer une expression d'une équation de phase F(.PHI.) = ~ dudit signal élémentaire ; et - des moyens pour déterminer une expression de la phase (.PHI.)(t) dudit signal élémentaire en fonction de paramètres (r, rk, .PHI.0, pk) mesurant l'anharmonicité dudit signal élémentaire et sa morphologie, à partir des fonctions pcosn et psinn définies par :
où r est un paramètre d'anharmonicité mesurant le degré d'anharmonicité du signal, qui varie dans l'intervalle [0,1], la limite r=0 correspondant à un signal harmonique, la limite r=1 à un signal infiniment anharmonique, et où k et n sont des entiers.
8.- Système selon la revendication 7, qui comporte des moyens pour exprimer l'équation de phase sous la forme :
dans laquelle r, variant dans [0,1], est un paramètre mesurant l'anharmonicité
dudit signal élémentaire.
9.- Système selon la revendication 8, qui comporte des moyens pour exprimer ledit signal élémentaire au moyen de deux paramètres r et .PHI.0, sous la forme :
x(t) = x0 + a1hsin(t,r)+ b1h cos(t,r) où a1 = x1 cos(.PHI.0) et b1 = ¨x1 sin(.PHI.0), les fonctions hsin et hcos étant définies par :
10.- Système selon la revendication 9, qui comporte des moyens pour exprimer l'équation de phase sous la forme:
dans laquelle P(.PHI.) et Q(.PHI.) sont des polynômes trigonométriques.
11.- Système selon la revendication 10, qui comporte des moyens pour exprimer la phase .PHI.(t) sous la forme :
dans laquelle les fonctions psin1 et pcos1 sont définies par :
12. ¨ Système selon l'une quelconque des revendications 7 à 11, dans lequel ladite interface comprend un masque respiratoire pour être placé sur le nez et la bouche du patient.
13.- Dispositif d'assistance respiratoire comprenant un système de mesure de l'activité respiratoire selon l'une quelconque des revendications 7 à 12.
14.- Dispositif de surveillance respiratoire comprenant un système de mesure de l'activité respiratoire selon l'une quelconque des revendications 7 à 12.
CA2797886A 2010-04-30 2011-04-28 Procede et systeme d'analyse de l'activite respiratoire d'un patient et applications correspondantes Active CA2797886C (fr)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
FR1053357A FR2959407A1 (fr) 2010-04-30 2010-04-30 Procede et systeme d'analyse de l'activite respiratoire d'un patient et applications correspondantes
FR1053357 2010-04-30
PCT/FR2011/050957 WO2011135257A2 (fr) 2010-04-30 2011-04-28 Procédé et système d'analyse de l'activité respiratoire d'un patient et applications correspondantes

Publications (2)

Publication Number Publication Date
CA2797886A1 CA2797886A1 (fr) 2011-11-03
CA2797886C true CA2797886C (fr) 2018-07-03

Family

ID=43661798

Family Applications (1)

Application Number Title Priority Date Filing Date
CA2797886A Active CA2797886C (fr) 2010-04-30 2011-04-28 Procede et systeme d'analyse de l'activite respiratoire d'un patient et applications correspondantes

Country Status (6)

Country Link
US (1) US9649051B2 (fr)
EP (1) EP2564326A2 (fr)
JP (1) JP5932769B2 (fr)
CA (1) CA2797886C (fr)
FR (1) FR2959407A1 (fr)
WO (1) WO2011135257A2 (fr)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
NZ736250A (en) * 2012-04-13 2019-04-26 ResMed Pty Ltd Apparatus and methods for ventilatory treatment

Family Cites Families (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5020540A (en) 1987-10-09 1991-06-04 Biometrak Corporation Cardiac biopotential analysis system and method
US5490502A (en) * 1992-05-07 1996-02-13 New York University Method and apparatus for optimizing the continuous positive airway pressure for treating obstructive sleep apnea
US6306088B1 (en) * 1998-10-03 2001-10-23 Individual Monitoring Systems, Inc. Ambulatory distributed recorders system for diagnosing medical disorders
US6644312B2 (en) * 2000-03-07 2003-11-11 Resmed Limited Determining suitable ventilator settings for patients with alveolar hypoventilation during sleep
JP5451013B2 (ja) * 2007-09-06 2014-03-26 英次 麻野井 睡眠の質を評価するために用いる表示又は印刷を行う装置、方法、コンピュータプログラム、呼吸補助装置、慢性心疾患患者を対象とした呼吸補助装置、睡眠導入装置、マッサージ装置、検査装置
JP5207172B2 (ja) * 2008-03-21 2013-06-12 国立大学法人 大分大学 波形解析装置及び波形解析プログラム
FR2955187A1 (fr) * 2010-01-08 2011-07-15 Centre Nat Rech Scient Procede de decomposition d'un signal periodique anharmonique et programme d'ordinateur correspondant
FR2958529A1 (fr) * 2010-04-13 2011-10-14 Centre Nat Rech Scient Procede et systeme d'analyse de l'activite cardiaque d'un patient et applications correspondantes

Also Published As

Publication number Publication date
WO2011135257A2 (fr) 2011-11-03
EP2564326A2 (fr) 2013-03-06
FR2959407A1 (fr) 2011-11-04
CA2797886A1 (fr) 2011-11-03
US20130096452A1 (en) 2013-04-18
JP5932769B2 (ja) 2016-06-08
WO2011135257A3 (fr) 2012-01-12
US9649051B2 (en) 2017-05-16
JP2013528415A (ja) 2013-07-11

Similar Documents

Publication Publication Date Title
CA2719169C (fr) Dispositif de detection de l&#39;observance par un patient d&#39;un traitement d&#39;oxygenotherapie
US9706946B2 (en) Spirometer system and methods of data analysis
US8425428B2 (en) Nitric oxide measurements in patients using flowfeedback
JP6199312B2 (ja) チェーンストークス呼吸パターンの判別
JP2008525060A5 (fr)
US20130190642A1 (en) Apparatus and method for diagnosing obstructive sleep apnea
EP3110321A1 (fr) Procede, dispositif, systeme et programme informatique de filtrage d&#39;une serie rr obtenue a partir d&#39;un signal cardiaque avec controle automatique de la qualite de la serie rr
EP2590701B1 (fr) Système de détection électroencéphalographique d&#39;une inadéquation entre l&#39;état d&#39;un patient placé sous assistance ventilatoire et le réglage de la machine utilisée pour cette assistance, et utilisation de cette détection pour l&#39;adaptation du réglage
EP0892618A1 (fr) Dispositif de determination de la profondeur d&#39;anesthesie
EP2364108A1 (fr) Systeme de detection de l&#39;activite musculaire respiratoire d&#39;un patient sous assistance respiratoire
CA2797886C (fr) Procede et systeme d&#39;analyse de l&#39;activite respiratoire d&#39;un patient et applications correspondantes
EP2364174B1 (fr) Système de quantification du desaccord entre un patient sous assistance respiratoire et un appareil d&#39;assistance correspondant
JP2007190082A (ja) 肺音診断装置及び肺音診断方法
EP3600010A1 (fr) Systeme de determination d&#39;un ensemble d&#39;au moins un descripteur cardio-respiratoire d&#39;un individu pendant son sommeil
Sa et al. Validation of the Forced Oscillation Technique in the diagnostic of respiratory changes in patients with silicosis
WO2015128568A1 (fr) Procede et dispositif de contrôle automatique de la qualite d&#39;une serie rr obtenue a partir d&#39;un signal cardiaque
WO2007097634A1 (fr) Procédé et appareil de surveillance d&#39;un patient sous sédation
CA2796172C (fr) Procede et systeme d&#39;analyse de l&#39;activite cardiaque d&#39;un patient et applications correspondantes
Galli et al. Accurate ECG monitoring by Gaussian feature streaming
Terrill et al. Investigating parameters participating in the infant respiratory control system attractor
Uchida et al. An ICA approach to reject noise from pressure changes of pillow
Ko et al. An EEG Spike Prescreening Algorithm Using RBF ANN with Multi-Channel Correlation

Legal Events

Date Code Title Description
EEER Examination request

Effective date: 20160321