FR 3080929 A1 20191108 FR 1853796 A 20180502 PROCEDE DE TRAITEMENT D'UN SIGNAL ACCELEROMETRIQUE. DOMAINE DE L'INVENTION La présente invention concerne le domaine du suivi respiratoire (indistinctement ventilatoire) d'un individu. Plus précisément, elle se rapporte dans le domaine du bien être et de la santé, à la mesure et à la caractérisation de la fonction ventilatoire d'un individu grâce à un ensemble 10 d'accéléromètres en réseau. La présente invention permet le recueil de données relatives à la respiration d'un individu. L'analyse ultérieure de ces données peut être utile par exemple pour la caractérisation d'une dynamique ventilatoire ou d'un diagnostic médical. 15 La présente invention trouve des applications par exemple dans : l'analyse du sommeil et la recherche de troubles respiratoires du sommeil, le sport et l'exploration des fonctions respiratoires, en particulier au regard des performances sportives aérobies de l'individu, mais également pour l'entraînement à la 20 plongée en apnée, le bien-être du grand public, par exemple pour l'évaluation de la dynamique respiratoire en position statique, le suivi de la dynamique respiratoire dans les pratiques de relaxation comme le yoga, l'analyse de la réponse d'un individu à un stress donné, par exemple pour une campagne 25 vidéo, un film, une musique, etc. le domaine de la recherche clinique et pré-clinique pour l'évaluation et la validation de l'effet de certains médicaments sur les fonctions respiratoires, le domaine de médecine générale et hospitalière pour le dépistage rapide de pathologies respiratoire. 30 A l'heure actuelle, il existe des bandes de pléthysmographie respiratoire qui se fixent autour du tronc d'un individu. Toutefois, de telles bandes sont peu pratiques à l'utilisation car elles peuvent bouger et/ou se détériorer pendant leur utilisation, ce qui rend alors inexploitables leurs données, et leur fragilité intrinsèque rend leur utilisation délicate. 35 Il apparaît donc souhaitable d'améliorer cette solution existante. 2 RESUME DE L'INVENTION Plus précisément, l'invention concerne un procédé de traitement de signal, comprenant des étapes consistant à: 5 - enregistrer préalablement dans une mémoire un premier signal initial issu d'un premier ensemble d'au moins un accéléromètre (110) tridimensionnel positionné en position thoracique d'un individu. Il est essentiellement caractérisé en ce qu'il comprend en outre : - enregistrer préalablement dans une mémoire un deuxième signal initial issu d'un deuxième ensemble d'au moins un accéléromètre (120) tridimensionnel positionné en position abdominale de l'individu, le deuxième signal étant synchronisé avec le premier signal, - traiter les données du premier signal initial pour calculer un premier vecteur final, représentatif des efforts thoraciques subis par le premier ensemble d'au moins un accéléromètre 15 tridimensionnel, - traiter les données du deuxième signal initial pour calculer un deuxième vecteur final, représentatif des efforts abdominaux subis par le deuxième ensemble d'au moins un accéléromètre tridimensionnel. 20 De préférence, l'étape de traitement des données du premier signal initial, respectivement du deuxième signal initial, comprend une étape de détection d'une fenêtre d'instabilité temporelle dudit premier signal initial, respectivement deuxième signal initial et, en cas de détection d'une fenêtre d'instabilité temporelle, le procédé comprenant alors en outre une étape d'inhibition temporaire du calcul dudit premier vecteur final, respectivement dudit deuxième vecteur final, de 25 préférence jusqu'à retrouver une fenêtre de stabilité temporelle. On peut prévoir que les données du premier signal initial, respectivement du deuxième signal, sont enregistrées sous forme d'une première matrice, respectivement d'une deuxième matrice ; et suivi par au moins une des étapes parmi : 30 - traiter les données du premier signal initial et - traiter les données du deuxième signal initial, comprenant une étape d'analyse en composantes principales des vecteurs de la première matrice, respectivement deuxième matrice, pour obtenir ledit premier vecteur final, respectivement ledit deuxième vecteur final ; 35 le procédé comprenant de préférence en outre une étape de filtrage du premier vecteur final. 3 On peut prévoir en outre une étape consistant à calculer à chaque temps de mesure l'angle de rotation par rapport à la gravité dudit premier, respectivement deuxième, ensemble d'au moins un accéléromètre tridimensionnel pour obtenir ledit premier vecteur final, respectivement deuxième vecteur final. 5 On peut prévoir une étape initiale de prétraitement consistant à sous échantillonner les signaux de chaque voie de chaque accéléromètre. On peut prévoir une étape de filtrage passe-bande de l'un au moins parmi ledit premier 10 signal et ledit deuxième signal, de préférence avec fréquence de coupure basse égale à 0,05 Hz et une fréquence de coupure haute égale à 0,8 Hz ; et optionnellement avec un filtre à réponse impulsionnelle réelle. On peut prévoir une étape consistant à qualifier l'effort respiratoire global de l'individu. De préférence, la qualification de l'effort comprend une combinaison linéaire du premier vecteur final et du deuxième vecteur final. La qualification de l'effort peut aussi comprendre une étape consistant à calculer, pour l'un 20 au moins parmi le premier ensemble (110) et le deuxième ensemble (120) d'au moins un accéléromètre tridimensionnel, une corrélation croisée normalisée, pour caractériser le caractère périodique de l'effort respiratoire. On peut aussi prévoir que la qualification de l'effort comprend une étape consistant à 25 calculer au moins l'un des descripteurs parmi : - le déphasage éventuel entre le signal issu du premier ensemble d'au moins un accéléromètre et le signal issu du deuxième ensemble d'au moins un accéléromètre, - le couplage entre le signal issu du premier ensemble d'au moins un accéléromètre et le signal issu du deuxième ensemble d'au moins un accéléromètre par intercorrélation ou mutuelle 30 information croisée, et - un index de désynchronisation thoraco-abdominal. On peut prévoir en outre une étape de détection de changement de position de l'individu. 35 Selon l'invention, le couple d'accéléromètres permet à lui seul, de reconstituer l'équivalent d'une Pléthysmographie Respiratoire par Inductance (PRI). En ce sens, l'invention peut être 4 considérée comme un pléthysmographe respiratoire virtuel car composé au minimum uniquement de deux accéléromètres, un posé sur le thorax du sujet et l'autre sur l'abdomen du sujet. L'invention trouve un intérêt notamment dans la miniaturisation du système de mesure, 5 moins complexe, moins coûteux, non intrusif et moins invasif que les systèmes de référence en hôpital, et dans sa capacité d'accéder à une mesure fiable, robuste et rapide. D'autres caractéristiques et avantages de la présente invention apparaîtront plus clairement à la lecture de la description suivante donnée à titre d'exemple illustratif et non limitatif et faite en 10 référence aux figures annexées. DESCRIPTIF DES DESSINS la figure 1 illustre un ensemble 6 de signaux bruts issus selon les 3 voies d'un 15 accéléromètre thoracique et d'un accéléromètre abdominal, ainsi qu'un signal de référence obtenu par une bande pléthysmographique abdominale et un signal de référence obtenu par une bande pléthysmographique abdominale, la figure 2 illustre de manière synchrone un signal PHI expérimental et un signal PRI de référence, obtenus selon les signaux de la figure 1, la figure 3 illustre la concordance entre les deux signaux PHI et PRI de la figure 2 selon un graphique de Bland-Altman, la figure 4 illustre de manière synchrone un signal dPHI, dérivé du signal PHI de la figure 2, un signal dPRI, dérivé du signal PRI de la figure 2, et un signal FLOW de référence, la figure 5 illustre de manière synchrone un signal sPHI et un signal sPRI, la figure 6 illustre le calcul d'un index de caractérisation de la périodicité d'un effort ventilatoire thoracique, la figure 7 illustre le calcul d'un index de caractérisation de la périodicité d'un autre effort ventilatoire thoracique, la figure 8A illustre un mode de réalisation d'un dispositif de mesure selon l'invention la figure 8B illustre un autre mode de réalisation d'un dispositif de mesure selon l'invention DESCRIPTION DETAILLEE Par concision, on entend par : « accéléromètre thoracique » ou « premier accéléromètre », un premier ensemble d'au moins un accéléromètre tridimensionnel positionné en position thoracique du sujet, 5 - « accéléromètre abdominal » ou « deuxième accéléromètre »: un deuxième ensemble d'au moins un accéléromètre tridimensionnel positionné en position abdominal du sujet, - « Individu », « patient » ou « sujet », toute personne physique, sain ou atteint d'une pathologie quelconque, 5 - « signaux » ou données, les valeurs mesurées selon chaque axe du premier et du deuxième accéléromètre, - « compartiment » le thorax ou l'abdomen. Dispositif de mesure 10 On prévoit un dispositif de mesure unique et compact comportant au moins deux accéléromètres configurés pour mesurer l'accélération selon trois axes. En particulier, on prévoit un premier ensemble d'au moins un premier accéléromètre, destiné 15 à être placé en position thoracique de l'individu, de préférence en position xiphoïde. On prévoit également un deuxième ensemble d'au moins un accéléromètre, destiné à être placé en position abdominale dudit individu. 20 Chaque ensemble d'au moins un accéléromètre peut comprendre une pluralité d'accéléromètres par exemple pour des raisons de redondance et robustesse. Par concision, le premier ensemble d'au moins un accéléromètre est appelé premier accéléromètre 110, et le deuxième ensemble d'au moins un accéléromètre est appelé deuxième 25 accéléromètre 120. En l'espèce, un seul accéléromètre par ensemble peut suffire, c'est à dire un unique premier accéléromètre placé en position thoracique et un unique deuxième accéléromètre en position abdominale. 30 Chaque accéléromètre émet un signal respectif comprenant un ensemble de données et appelé « signal initial », qui correspond à l'évolution des accélérations (en g) en fonction du temps. Par concision, on entend indistinctement « signal » et « données ». 35 De préférence, on prévoit que le dispositif de mesure comprend des moyens d'activation / désactivation des accéléromètres, par exemple sous forme d'interrupteur de commande d'activation, analogique ou numérique, éventuellement sans fil. 6 Dans une première variante, on prévoit que le premier accéléromètre et le deuxième accéléromètre sont solidaires d'un unique support 100 (figure 8A). 5 Dans une seconde variante, on prévoit que le premier accéléromètre et le deuxième accéléromètre sont solidaires d'un support respectif 200, 210 (figure 8B). Quelle que soit la variante, chaque support comprend une face supérieure et une face inférieure. La face supérieure supporte au moins un accéléromètre et la face inférieure, destinée à 10 être au contact de l'individu, est typiquement enduite d'adhésif (colle), de préférence repositionnable. L'avantage de l'aspect repositionnable est que l'individu peut retirer le dispositif de mesure au réveil, après une douche, etc. et réutiliser le même dispositif de mesure ensuite, et ce sur plusieurs jours consécutifs. Le dispositif est par exemple alimenté par une batterie. 15 Le premier accéléromètre et le deuxième accéléromètre peuvent être identiques. En l'espèce, le premier accéléromètre et le deuxième accéléromètre sont des accéléromètres trois axes standard. Les données issues de chaque capteur sont enregistrées dans une mémoire informatique. 20 Elles sont envoyées par voie filaire ou sans fil. Dans une première variante, la mémoire est en connexion informatique filaire avec le premier et le deuxième accéléromètre. La mémoire peut par exemple être disposée sur le support unique ou l'un quelconque des supports des accéléromètres. 25 De préférence, on prévoit que la mémoire est en connexion informatique avec un port de communication d'entrée/sortie (I/O), par exemple un port USB, ce qui permet d'exporter les données enregistrées vers un dispositif de traitement des données comprenant un logiciel de traitement des données, typiquement tout objet communicant, c'est-à-dire un dispositif 30 électronique comprenant des moyens de communication avec ou sans fil, un processeur et de préférence un écran d'affichage, par exemple un ordinateur personnel, un smartphone, une tablette tactile etc. On peut aussi prévoir que la mémoire est amovible, par exemple sous forme de support de 35 donnée de type carte mémoire à connexion USB, connectable au premier accéléromètre et au deuxième accéléromètre pour l'enregistrement des données issues de ceux-ci, et connectable ensuite à un ordinateur équipé d'un logiciel de filtrage permettant de filtrer lesdites données. 7 Dans une deuxième variante, la mémoire est déportée, c'est-à-dire sans connexion filaire entre ladite mémoire et les capteurs. Par exemple, la mémoire est déportée dans le dispositif de traitement des données. 5 Dans ce cas, on prévoit avantageusement que le dispositif de mesure comprend des moyens de communication avec ledit dispositif de traitement des données, que ce soit une communication filaire ou sans fil. 10 La mémoire peut donc être embarquée ou déportée, sur un serveur, en particulier d'un système informatique en nuage. Les accéléromètres sont montés avantageusement sur une carte électronique, comportant une batterie, un processeur (microcontrôleur) et une unité de stockage pouvant stocker les 15 données dans l'optique d'un traitement embarqué. Alternativement, chaque accéléromètre est monté sur une carte électronique indépendante, chacune disposant d'une pile d'alimentation, d'un processeur et d'une unité de stockage. On prévoit que les capteurs optionnels sont solidaires de la carte électronique du premier ensemble d'au moins un accéléromètre La synchronisation des signaux est faite par construction par le microcontrôleur. 20 De préférence, on prévoit également des moyens de traitement des signaux issus des capteurs. En l'espèce, les moyens de traitement des signaux comprennent des moyens de filtrage, 25 configurés pour filtrer le signal initial du premier accéléromètre, et configurés pour filtrer le signal initial du deuxième accéléromètre. De manière similaire à la mémoire, les moyens de filtrage peuvent être intégrés au dispositif de mesure, ou déportés sur le dispositif de traitement des données ou encore sur un serveur en 30 communication avec le dispositif de traitement des données. Dans ce cas, on prévoit avantageusement que le dispositif de mesure comprend des moyens de communication avec lesdits moyens de traitement des signaux, que ce soit une communication filaire ou sans fil. 35 8 Chaque accéléromètre tridimensionnel possède 3 voies (ou axes) orthonormées et appelées respectivement X, Y et Z, dont il conviendra de déterminer la position angulaire par rapport à la gravité. 5 Par convention, on entend par : - X1, Y1 et Z1 les 3 voies du premier accéléromètre, - X2, Y2 et Z2 les 3 voies du deuxième accéléromètre. A noter qu'une fois le premier accéléromètre et le deuxième accéléromètre positionnés sur 10 un individu, les axes X1, Y1 et Z1 ne sont pas nécessairement respectivement confondus avec les axes X2, Y2 et Z2. Pour chaque accéléromètre, les données échantillonnées sont enregistrées sous forme de matrice de départ à trois colonnes, une colonne pour l'axe X, une colonne pour l'axe Y et une 15 colonne pour l'axe Z; la matrice comprenant autant de lignes que de points d'échantillonnage, par exemple, L=60*T*F avec L le nombre total de lignes, T le temps d'échantillonnage en minutes et F la fréquence d'échantillonnage. On a donc une matrice de départ du compartiment thoracique et une matrice de départ du 20 compartiment abdominal. La respiration de l'individu génère des mouvements de sa cage thoracique et/ou de son abdomen. Le premier accéléromètre et le deuxième accéléromètre permettent, comme décrit ultérieurement, de reconstruire ces mouvements, notamment par calcul de leur vitesse angulaire 25 respective. Comme les axes XYZ sont orthonormés, la connaissance de la position de deux axes parmi X, Y et Z par rapport à la gravité permet de reconstruire l'autre. Et lorsqu'il existe une convention de pose déterminée de l'accéléromètre sur l'individu, cela permet de fixer le référentiel individu 30 /accéléromètre, alors dans ce cas la connaissance de la position d'un seul axe parmi X, Y et Z par rapport à la gravité permet de reconstruire les 2 autres. Quand un individu respire, le déplacement de la cage thoracique ou de l'abdomen est très faible. Par conséquent, la vitesse de déplacement, donc l'accélération est négligeable. Le signal 35 issu d'un accéléromètre n'est donc pas réellement lié au déplacement de la cage thoracique ou au déplacement de l'abdomen au premier ordre. 9 En revanche, quand un individu respire, les accéléromètres (thoracique ou abdominal) subissent une rotation par rapport à la gravité. Quand un accéléromètre pivote, le vecteur gravité pivote également dans le référentiel X1Y1Z1 du premier accéléromètre et dans le référentiel X2Y2Z2 du deuxième accéléromètre. 5 La présente invention vise à suivre ce pivotement, et à calculer l'évolution de l'orientation du vecteur gravité dans le référentiel X1Y1Z1 ou X2Y2Z2 au cours du temps. De préférence, on considère que l'individu est globalement statique pendant l'enregistrement 10 des données issues des accéléromètres, de sorte à obtenir une stationnarité des signaux, ce qui permet d'éviter que le mouvement de l'individu ne vienne bruiter les mesures des accéléromètres. Comme décrit ci-après, le principe de l'invention consiste à mesurer les variations d'inclinaison des accéléromètres par rapport à la gravité en fonction d'un axe de référence à 15 chaque changement de stationnarité du signal, puis à traiter ces signaux pour reconstruire la ventilation de l'individu. Prétraitement 20 On prévoit une étape de prétraitement consistant à sous-échantillonner les signaux issus de chaque voie des accéléromètres. De préférence, on prévoit que les signaux issus du premier accéléromètre et les signaux issus du deuxième accéléromètre sont échantillonnés à la même fréquence d'échantillonnage Fe. De préférence, on prévoit que la fréquence d'échantillonnage Fe est inférieure à une valeur seuil prédéterminée, de sorte à sous-échantillonner lesdits signaux. Par exemple, Fe = 256 Hz. De préférence, on prévoit également que la fréquence d'échantillonnage Fe est supérieure à une autre valeur seuil prédéterminée. Par exemple, on prévoit que la Fréquence de Nyquist = Fe/2 est supérieure à la fréquence de coupure haute multipliée par une valeur de référence, en l'espèce 10, de sorte à obtenir une fréquence de Nyquist au moins égale à 8 Hz. 25 30 35 10 Filtrage Après l'étape de prétraitement, on prévoit une étape de filtrage des signaux. 5 On prévoit de filtrer chaque voie des accéléromètres par un même filtre. De préférence, on prévoit un filtrage passe-bande. On peut prévoir une fréquence de coupure basse comprise entre 0,01 et 0,1 Hz et de fréquence de coupure haute comprise entre 0,6 et 1 Hz. En l'espèce on prévoit un filtre à réponse impulsionnelle réelle (FI R) passe bande de fréquence de coupure basse égale à 0,05 Hz et de fréquence de coupure haute égale à 0,8 Hz. 10 Reconstruction Après l'étape de filtrage, on prévoit une étape de reconstruction des efforts respiratoires, c'est à dire des mouvements à l'origine des signaux issus du premier accéléromètre (thoracique) 15 et des signaux issus du deuxième accéléromètre (abdominal), selon l'une au moins des deux variantes ci-dessous. De préférence, chacune des deux variantes ci-dessous est mise en oeuvre après chaque détection de changement de position de l'individu. 1e variante : analyse en composantes principales (ACP) On peut prévoir une étape d'analyse en composantes principales de chaque matrice de départ, pour chaque compartiment, et réduire ainsi la dimensionnalité de chaque matrice. En l'espèce, cela permet de passer d'une matrice L*3 à un vecteur L*1. Pour les données issues d'un accéléromètre donné, une analyse en composantes principales de ces données est une projection dans un nouvel espace, c'est à dire une 30 combinaison linéaire de chaque variable initiale, dans un nouvel espace qui décrit le maximum de variance dans le repère XYZ dudit accéléromètre. Ainsi, à chaque instant, l'ACP permet de déterminer l'axe X1, Y1 ou Z1 du premier accéléromètre et l'axe X2, Y2 ou Z2 du deuxième accéléromètre qui décrit le maximum de 35 variance (au moins 80%) et qui pivote par rapport à l'axe de la gravité. 20 25 11 Les étapes précédentes ont permis de filtrer (supprimer) les bruits d'origines physiologiques (par exemple l'activité cardiaque) et les artefacts divers (mouvements , bruits, etc...) et conserver l'activité respiratoire sur les signaux accélérométriques, autrement dit sur la mesure du pivotement du premier et du deuxième accéléromètre. Par conséquent, l'axe qui présente le maximum de 5 variabilité est celui qui porte la respiration. En l'espèce, pour chaque accéléromètre, les données échantillonnées sont enregistrées sous forme de matrice de départ à trois colonnes, une colonne pour l'axe X, une colonne pour l'axe Y et une colonne pour l'axe Z; la matrice comprenant autant de lignes que de points 10 d'échantillonnage, par exemple, L=60*T*F avec L le nombre total de lignes, T le temps d'échantillonnage en minutes et F la fréquence d'échantillonnage. L'ACP consiste à transformer la matrice de départ en un vecteur d'arrivée dont chaque coefficient est une combinaison linéaire respective de la ligne de la matrice de départ 15 correspondante, donc des 3 axes de la matrice de départ. En l'espèce, l'ACP consiste en la mesure des valeurs propres et des vecteurs propres de la matrice de variance-covariance de la matrice de départ. Le vecteur propre dont la valeur propre est la plus élevée est considérée comme la composante principale et servira de vecteur pour 20 projection dans un nouvel espace où chaque coordonnée est une combinaison linéaire de chaque axe (X, Y et Z de départ). On peut alors calculer le vecteur A_Respi résultant de la multiplication de la matrice de départ A_PCA par la composante principale CP : CF 25 APCA p, = Apc A* CP = [xi Yi 30 12 Chaque coordonnée du vecteur résultant A_Respi est donc, avec comme exemple à la coordonnée i : AResp- CP - 5 On obtient ainsi : un premier vecteur final, résultant d'une analyse en composantes principales sur une première matrice de départ correspondant aux données issues du premier accéléromètre, et un deuxième vecteur final, résultant d'une analyse en composantes principales sur une 10 deuxième matrice de départ correspondant aux données issues du deuxième accéléromètre. Chaque vecteur final est une reconstruction : le premier vecteur final est une reconstruction de l'effort respiratoire du compartiment thoracique et le deuxième vecteur final est une 15 reconstruction de l'effort respiratoire du compartiment abdominale. Toutefois, l'ACP est une méthode qui est sensible au bruit. Le filtrage selon l'invention précédemment décrit permet de réduire le bruit. Cependant, il peut être avantageux, en outre, de détecter certains événements générateurs de bruit, par exemple un changement de position de 20 l'individu. A cet effet, on peut prévoir de mesurer l'énergie du signal du premier et/ou du deuxième accéléromètre sur un ensemble de fenêtres temporelles glissantes, de préférence deux fenêtres temporelles adjacentes se recouvrant partiellement. 25 On prévoit alors une étape de seuillage consistant à comparer l'énergie du signal du premier et/ou du deuxième accéléromètre à une valeur seuil prédéterminée, pour chaque fenêtre temporelle, et si l'énergie du signal est supérieure à la valeur seuil prédéterminée, alors les données de cette fenêtre temporelle ne sont pas prises en compte dans l'ACP (i.e. dans la mesure 30 du vecteur propre servant de base pour la projection dans le nouvel espace). 35 Filtrage final On peut alors prévoir ensuite une étape de filtrage final du vecteur final. De préférence, le filtrage final est le même que le filtrage passe bande décrit précédemment. 13 2ème variante : méthode angulaire Il peut également être avantageux, en remplacement ou en complément de la méthode ACP, 5 de mettre en oeuvre une méthode dite angulaire. Dans cette variante, on vise à calculer l'angle de rotation de chaque accéléromètre par rapport à la gravité à chaque temps de mesure. 10 A cet effet, on prévoit de mesurer et de normaliser le vecteur accélération, noté a. En conditions normales, l'utilisateur est statique et les accélérations dues à la respiration sont négligeables. Par conséquent, en conditions normales, le vecteur accélération a est sensiblement égal à l'accélération due à la gravité g. Quand un accéléromètre pivote, le vecteur gravité pivote également dans le référentiel intrinsèque dudit accéléromètre. Un tel pivotement est du à la respiration et varie donc au cours du temps. Par conséquent, on peut considérer que le mouvement angulaire dudit accéléromètre a lieu selon un seul axe de rotation r. Il s'agit donc de mesurer le mouvement angulaire de l'axe de rotation r dudit accéléromètre. L'angle de rotation et et l'axe de rotation rt du vecteur accélération a entre deux mesures consécutives au temps t-1 et au temps t sont donnés respectivement par le produit scalaire et le 25 produit vectoriel des deux vecteurs suivants : cos - at-i) 61?-1 La respiration est un mouvement oscillatoire. Par conséquent, l'axe de rotation r du vecteur accélération a s'inverse quand la direction de rotation s'inverse. Il apparaît donc utile de normaliser cette direction de l'axe de rotation r dans un hémisphère prédéterminé, par comparaison à un axe 30 de référence rref. On peut ainsi définir : 15 20 14 De préférence, afin de diminuer l'influence du bruit et de favoriser l'influence des mesures proches du temps t, on prévoit de déterminer l'axe de rotation r normalisé moyenné sur une fenêtre de temps de durée VV tel que : " normalise 5 Avec H(n) une fenêtre de Hamming. On peut utiliser d'autres fenêtres, par exemple une fenêtre rectangulaire, une fenêtre triangulaire, une fenêtre de Hann, une fenêtre de Blackman, une fenêtre de Kaiser, etc. 10 De manière similaire, on prévoit de déterminer at l'accélération moyenne normalisée moyennée sur une fenêtre de temps de durée VV, de préférence la même fenêtre temporelle que pour , tel que : = non tal-,r ( .p f f On peut alors calculer l'angle de rotation dudit accéléromètre tel que : 1 On peut aussi calculer, si nécessaire, la vitesse angulaire wt telle que : d(4, (1t. Au minimum les valeurs de l'angle de rotation sont enregistrées sous forme de vecteur PRI virtuelle Ainsi, quelle que soit la première variante, la deuxième variante ou la combinaison des deux variantes mise en oeuvre, on obtient un premier vecteur final qui qualifie l'effort thoracique et un 25 deuxième vecteur final qui qualifie l'effort abdominal, ce qui permet d'obtenir une PRI virtuelle. at+i 15 final. 20 15 Grâce à ce premier vecteur final et ce deuxième vecteur final, on peut qualifier l'effort respiratoire global de l'individu, c'est à dire le volume relatif, ainsi que la dynamique temporelle des efforts de chaque compartiment, donc de qualifier leur cohérence temporelle. 5 La figure 1 illustre de manière synchrone : les données brutes AX_T, AY_T et AZ_T issues du premier accéléromètre selon les axes XYZ dudit premier accéléromètre ; les données brutes AX_A, AY_A et AZ_A issues du deuxième accéléromètre, selon les axes XYZ dudit deuxième accéléromètre ; 10 les données brutes PRI_T et PRI_A d'une bande de pléthysmographie positionnée respectivement en position thoracique et en position abdominale. La PRI thoracique (PRI_T) et la PRI abdominale (PRI_A) illustrées en figure 1 sont des variations de surface de section et servent de référence. 15 La somme de la PRI thoracique (PRI_T) et de la PRI abdominale (PRI_A) est un signal PRI qui donne une image du volume ventilé courant. La figure 2 illustre la superposition normalisée entre les valeurs -1 et +1 pour un signal PHI 20 expérimental et un signal PRI de référence. Le signal PRI de référence est la somme des données brutes PRI_T et PRI_A de la figure 1. Le signal PHI est calculé, en l'espèce selon la méthode angulaire, à partir des données brutes AX_T, AY_T et AZ_T issues du premier accéléromètre et des données brutes AX_A, AY_A 25 et AZ_A issues du deuxième accéléromètre. Il est clair que la figure 2 illustre la pertinence de la présente invention. Les deux signaux PHI et PRI sont quasiment toujours en phase. En l'espèce, la concordance entre les deux signaux PHI et PRI peut être établie, par exemple selon un graphique de Bland-Altman, voir figure 3, pour 30 lequel on obtient une valeur de r2 = 0.78. Reconstruction d'une image du volume Il existe une correspondance entre la vitesse angulaire d'un accéléromètre thoracique et 35 abdominal et la forme d'onde de l'effort ventilatoire globale respiratoire. 16 On peut donc prévoir une étape de calcul de la vitesse angulaire de chaque accéléromètre, ce qui permet une reconstruction d'une image du débit respiratoire, une reconstruction de la dynamique des mouvements respiratoires de l'individu. 5 A cet effet, on prévoit au moins l'une des combinaisons linéaires parmi : - une combinaison linéaire de l'ACP du premier accéléromètre (thoracique) et de l'ACP du deuxième accéléromètre (abdominal) ; et - une combinaison linéaire de la méthode angulaire appliquée au premier accéléromètre (thoracique) et de la méthode angulaire appliquée au deuxième accéléromètre 10 (abdominal). La figure 4 illustre de manière synchrone : - le signal dPHI qui est la dérivée par rapport au temps du signal PHI de la figure 2, - le signal dPRI qui est la dérivée par rapport au temps du signal PRI de la figure 2, et 15 - un signal FLOW de référence, qui est un signal de débit réel mesuré en l'espèce grâce à un masque respiratoire posé sur un individu. La figure 4 illustre également clairement la pertinence de la présente invention puisque la le calcul de la dérivé n'amène pas de bruit particulier. Les 3 signaux dPHI, dPRI et FLOW sont 20 globalement en phase. La présente invention permet de calculer une image du volume relatif, ainsi que sa variation en fonction du temps. En ce sens, elle est une PRI virtuelle. Le premier accéléromètre et le deuxième accéléromètre fonctionnent comme un réseau de capteurs. Leurs mesures respectives sont synchronisées et correspondent respectivement à 25 l'accélération du compartiment thoracique et du compartiment abdominal. Il est alors possible de calculer les vitesses angulaires respectives. Ensuite, on peut comparer les variations de la vitesse angulaire du premier accéléromètre et les variations de la vitesse angulaire du deuxième accéléromètre sur une même base temporelle, et détecter par 30 exemple si le premier accéléromètre et le deuxième accéléromètre sont en déphasage ou en opposition de phase et calculer par exemple un index de désynchronisation abdomino-thoracique. De manière similaire, la figure 5 illustre de manière synchrone un signal sPHI, qui est la somme des signaux PHI en l'espèce de la figure 1, et un signal sPRI, qui est la somme des 35 signaux PRI en l'espèce de la figure 1. 17 Sur la figure 5, il est clair qu'une modification du rythme respiratoire (zone PB) est détectée par le procédé selon l'invention. La présente invention permet de reconstruire une « image volume » des compartiments 5 respiratoires et de suivre des évènements d'intérêt respiratoire. Elle présente une robustesse algorithmique et physiologique. Comme exposé précédemment la première variante (ACP) et la deuxième variante (angulaire) sont combinables. Par exemple, on peut prévoir dans une première phase de mettre en oeuvre la première variante (ACP) sur une première fenêtre temporelle des mesures de l'un au moins des accéléromètres parmi le premier accéléromètre et le deuxième accéléromètre, puis de mettre en oeuvre la deuxième variante (angulaire) sur une sous-fenêtre temporelle de la première fenêtre. Qualification de l'effort respiratoire global Grâce à la présente invention, on peut, à l'aide d'un seul couple d'accéléromètres, qualifier le mode de fonctionnement ventilatoire d'un individu. Dans une première variante, on prévoit de calculer une combinaison linéaire du premier vecteur final et du deuxième vecteur final. Par exemple, on prévoit des coefficients unité dans la combinaison linéaire. 25 Dans une deuxième variante, on prévoit de calculer au moins l'un des descripteurs parmi : une corrélation croisée normalisée sur au moins l'un des vecteurs finaux afin de caractériser la périodicité de l'effort respiratoire pour caractériser la périodicité (ou non périodicité) de l'un au moins des signaux parmi le signal issu du premier accéléromètre et le signal issu du deuxième accéléromètre,; 30 o et par exemple comparer la périodicité (ou non périodicité) de l'un des deux signaux issu du premier accéléromètre et du deuxième accéléromètre par rapport à l'autre deux signaux issu du premier accéléromètre et du deuxième accéléromètre, la différence entre le premier vecteur final et le deuxième vecteur final, pour calculer au 35 moins l'un parmi : o le déphasage éventuel le signal issu du premier accéléromètre et le signal issu du deuxième accéléromètre, 10 15 20 18 o le couplage entre le signal issu du premier accéléromètre et le signal issu du deuxième accéléromètre par intercorrélation ou mutuelle information croisée, o un index de désynchronisation thoraco-abdominal, o et par exemple qualifier ainsi la cohérence du compartiment thoracique et du 5 compartiment abdominal, - une quantification de la contribution de chaque compartiment dans la ventilation globale. Instabilités 10 En outre, on peut prévoir une étape de détection d'instabilités, ou événements anormaux, par exemple non stationnaires, qui peuvent par exemple correspondre à un changement de position de l'individu et/ou de mouvement corporels, ce qui peut être le cas y compris la nuit. Ainsi chaque signal issu du premier accéléromètre et du deuxième accéléromètre comprend 15 des fenêtres temporelles de stabilité et éventuellement un ensemble d'au moins une fenêtre temporelle d'instabilité. De préférence, on prévoit de mettre en oeuvre les calculs décrits ci-avant dans une fenêtre de stabilité. 20 A cet effet, on peut prévoir de comparer la valeur du signal de l'un au moins des accéléromètres à une valeur seuil prédéterminée, et pour chaque fenêtre temporelle pour laquelle la valeur du signal est supérieure à la valeur seuil prédéterminée, alors les données de cette fenêtre temporelle ne sont pas prises en compte dans les calculs ; ainsi les calculs ne sont 25 effectués que sur des phases stationnaires pendant lesquelles l'individu est immobile, ce qui améliore le rapport signal / bruit. On peut alors lancer une nouvelle séquence d'étapes de prétraitement, filtrage, calculs de vecteurs finaux par ACP / méthode angulaire et qualification de l'effort respiratoire global. 30 Par exemple, la figure 6 et la figure 7 illustrent de manière synchrone : - un signal FLOW de référence, qui est un signal de débit réel mesuré en l'espèce grâce à un masque respiratoire posé sur un individu, - un signal PHI_T, qui est le signal PHI calculé en l'espèce selon la méthode angulaire, à 35 partir des données brutes AX_T, AY_T et AZ_T issues du premier accéléromètre, - un signal CROSS, qui est la corrélation croisée normalisée PHI_T, 19 un index INDEX_T, qui est la caractérisation de la périodicité de l'effort ventilatoire thoracique reconstruit PHI_T, l'index INDEX_T étant un signal binaire pour lequel la valeur 1 correspond à un signal PHI_T périodique, et la valeur 2 à un signal PHI_T non périodique, obtenu en l'espèce par comparaison de la valeur du signal CROSS à une 5 valeur seuil prédéterminée. Sur la figure 6 et sur la figure 7, le signal PHI_T présente un défaut de périodicité, ce défaut étant cerclé en pointillés sur les figures 6 et 7 (observable également sur la référence FLOVV). 10 Comme illustré par le signal INDEX_T, ce défaut de périodicité est clairement détecté par le procédé selon l'invention (front montant du signal INDEX_T sur la figure 6 et la figure 7). Un retour à la normale, c'est à dire un signal PHI_T redevenant périodique est également détecté grâce au procédé selon l'invention, comme illustré par le front descendant du signal 15 INDEX_T REVENDICATIONS1. Procédé de traitement de signal, comprenant des étapes consistant à: - enregistrer préalablement dans une mémoire un premier signal initial issu d'un premier ensemble d'au moins un accéléromètre (110) tridimensionnel positionné en position thoracique d'un individu, caractérisé en ce qu'il comprend en outre : - enregistrer préalablement dans une mémoire un deuxième signal initial issu d'un deuxième ensemble d'au moins un accéléromètre (120) tridimensionnel positionné en position abdominale de l'individu, le deuxième signal étant synchronisé avec le premier signal, traiter les données du premier signal initial pour calculer un premier vecteur final, représentatif des efforts thoraciques subis par le premier ensemble d'au moins un accéléromètre tridimensionnel, traiter les données du deuxième signal initial pour calculer un deuxième vecteur final, représentatif des efforts abdominaux subis par le deuxième ensemble d'au moins un accéléromètre tridimensionnel. 2. Procédé selon la revendication 1, dans lequel l'étape de traitement des données du premier signal initial, respectivement du deuxième signal initial, comprend une étape de détection d'une fenêtre d'instabilité temporelle dudit premier signal initial, respectivement deuxième signal initial et, en cas de détection d'une fenêtre d'instabilité temporelle, le procédé comprenant alors en outre une étape d'inhibition temporaire du calcul dudit premier vecteur final, respectivement dudit deuxième vecteur final, de préférence jusqu'à retrouver une fenêtre de stabilité temporelle. 3. Procédé selon l'une quelconque des revendications précédentes, dans lequel les données du premier signal initial, respectivement du deuxième signal, sont enregistrées sous forme d'une première matrice, respectivement d'une deuxième matrice ; Et dans lequel au moins une des étapes parmi : traiter les données du premier signal initial et traiter les données du deuxième signal initial, comprend une étape d'analyse en composantes principales de la première matrice, respectivement deuxième matrice, pour obtenir ledit premier vecteur final, respectivement ledit deuxième vecteur final ; le procédé comprenant optionnellement en outre une étape de filtrage du premier vecteur final. 21 4. Procédé selon l'une quelconque des revendications précédentes, comprenant en outre une étape consistant à calculer, à chaque temps de mesure, l'angle de rotation par rapport à la gravité dudit premier ensemble d'au moins un accéléromètre tridimensionnel, respectivement du deuxième ensemble d'au moins un accéléromètre tridimensionnel, pour 5 obtenir ledit premier vecteur final, respectivement ledit deuxième vecteur final. 5. Procédé selon rune quelconque des revendications précédentes, comprenant une étape de filtrage passe-bande de l'un au moins parmi ledit premier signal et ledit deuxième signal, de préférence avec fréquence de coupure basse égale à 0,05 Hz et une fréquence de 10 coupure haute égale à 0,8 Hz; et optionnellement avec un filtre à réponse impulsionnelle réelle. 6. Procédé selon l'une quelconque des revendications précédentes, comprenant une étape consistant à qualifier l'effort respiratoire global de l'individu. 15 7. Procédé selon la revendication 6, comprenant une combinaison linéaire du premier vecteur final et du deuxième vecteur final. 8. Procédé selon la revendication 6 ou 7, comprenant une étape consistant à 20 calculer, pour l'un au moins parmi le premier ensemble (110) et le deuxième ensemble (120) d'au moins un accéléromètre tridimensionnel, une corrélation croisée normalisée, pour caractériser le caractère périodique de l'effort respiratoire. 9. Procédé selon l'une quelconque des revendications 6 à 8, comprenant une étape 25 consistant à calculer au moins l'un des descripteurs parmi : le déphasage éventuel entre le signal issu du premier ensemble d'au moins un accéléromètre et le signal issu du deuxième ensemble d'au moins un accéléromètre, le couplage entre le signal issu du premier ensemble d'au moins un accéléromètre 30 et le signal issu du deuxième ensemble d'au moins un accéléromètre par intercorrélation ou mutuelle information croisée, et un index de désynchronisation thoraco-abdominal. 10. Procédé selon l'une quelconque des revendications précédentes, comprenant en 35 outre une étape de détection de changement de position de l'individu.