FR 2794961 A1 20001222 FR 9907632 A 19990616 La présente invention concerne un procédé de détermination du dé- calage temporel entre les instants de passage d'une même onde de pouls en deux points de mesure distincts d'un réseau artériel d'un être vivant. Elle concerne en outre un procédé de détermination de la vitesse de propagation de l'onde de pouls entre ces deux points de mesure, à partir du décalage temporel déterminé. Elle concerne enfin un procédé d'estimation de la pression artérielle pulsée aortique à partir de vitesses de propagation de l'onde de pouls dé- terminées. Une installation de mesure de la vitesse de propagation de l'onde de pouls entre deux points du corps est actuellement connue. Celle-ci met en oeuvre un procédé de détermination du décalage temporel entre les instants de passage d'une même onde de pouls en deux points de mesure distincts du réseau artériel d'un être humain. Ces points de mesure sont disposés par exemple sur l'artère carotide primitive au niveau du cou et sur l'artère fémo- rale, au niveau de l'aine. Ce dernier capteur ou un capteur supplémentaire peut également être disposé sur l'artère brachiale à l'extrémité distale de l'avant-bras. Pour déterminer la vitesse de propagation de l'onde de pouls, la dis- tance séparant les deux points de mesure est mesurée manuellement à l'aide d'un mètre à ruban. A partir de la distance séparant les deux points et du décalage tempo- rel déterminé, la vitesse de propagation est estimée par simple division de la distance par le décalage temporel, en tenant compte toutefois de la distance séparant le coeur de l'un des points de mesure ou d'un éventuel trajet effec- tué en une direction opposée par l'une des deux ondes par rapport à l'autre, cette distance devant être déduite de la distance séparant les deux points de mesure. Pour la détermination du décalage temporel entre les instants de pas- sage de la même onde de pouls en deux points distincts du réseau artériel, le procédé connu consiste à enregistrer l'évolution au cours du temps de la pression artérielle, à chaque point de mesure, puis à déduire, par analyse des deux enregistrements le décalage temporel entre les instants de pas- sage de la même onde de pouls. A cet effet, les courbes obtenues à partir de chaque enregistrement sont analysées. Ces courbes présentent, lors du passage de chaque onde de pouls un pic caractéristique représentatif d'une pulsation cardiaque. Selon le procédé actuel, afin de déterminer le décalage temporel, on considère les maxima de pression de pics correspondant à une même onde de pouls, et on détermine le décalage temporel entre ces deux maxima. Ce décalage peut être obtenu manuellement par mesure directe sur des cour- bes imprimées représentatives des deux enregistrements ou par traitement informatique des deux enregistrements. Toutefois, la détermination des maxima des pics est relativement déli- cate et variable selon plusieurs paramètres cliniques et biologiques. Ainsi, la précision de détermination des instants permettant de calculer le décalage temporel est faible, ce qui nuit à la qualité d'ensemble de la mesure. L'invention a pour but de fournir un procédé de précision accrue, no- tamment en améliorant la détermination des instants permettant le calcul du décalage temporel entre les passages d'une même onde de pouls en deux points distants du réseau artériel d'un être vivant. A cet effet, I'invention a pour objet un procédé de détermination du décalage temporel entre les instants de passage d'une même onde de pouls en deux points de mesure distincts d'un réseau artériel d'un être vivant, comportant les étapes de: a) placer un capteur d'un paramètre représentatif de l'onde de pouls à chaque point de mesure du réseau artériel; b) enregistrer l'évolution au cours du temps du paramètre représenta- tif de l'onde de pouls à chaque point de mesure; c) déduire des deux enregistrements ledit décalage temporel, caractérisé en ce que l'étape de déduction du décalage temporel comporte les étapes de: cl) déterminer, pour un premier enregistrement, un premier instant d'apparition de l'onde de pouls en identifiant le pied de l'onde de pouls, c2) déterminer, pour le deuxième enregistrement, un deuxième ins- tant d'apparition de la même onde de pouls en identifiant, dans le deuxième enregistrement, I'apparition de la partie correspondante de l'onde de pouls, c3) calculer le décalage temporel à partir des instants et déterminés sur les premier et deuxième enregistrements. Suivant des modes particuliers de mise en oeuvre, le procédé com- porte l'une ou plusieurs des caractéristiques suivantes: - le premier instant d'apparition du pied de l'onde de pouls est déter- miné comme l'instant o la dérivée de la fonction enregistrée est maximale; - le deuxième instant d'apparition de la même onde de pouls sur le deuxième enregistrement est obtenu par recherche, dans le deuxième enre- gistrement, de l'instant o un coefficient de corrélation est maximal entre une plage du premier enregistrement définie autour du premier instant et une plage correspondante du deuxième enregistrement; - ledit coefficient de corrélation est un coefficient de corrélation li- néaire; et -ledit paramètre représentatif de l'onde de pouls est la pression arté- rielle. L'invention a en outre pour objet un procédé de détermination de la vitesse de propagation d'une onde de pouls entre deux points de mesure distants d'un réseau artériel d'un être vivant, comportant les étapes de: A) mesurer la distance séparant les deux points de mesure sur l'être vivant, B) déterminer le décalage temporel entre les instants de passage d'une même onde de pouls aux deux points de mesure; C) calculer la vitesse de propagation à partir de la distance séparant les deux points de mesure et du décalage temporel, caractérisé en ce que le décalage temporel entre les instants de passage de l'onde de pouls aux deux points de mesure est déterminé par un procédé selon l'une quelconque des revendications précédentes. L'invention a enfin pour objet un procédé d'estimation de la pression pulsée aortique, caractérisé en ce qu'il comporte les étapes de: a) déterminer la vitesse de propagation aortique d'une onde de pouls entre un point de l'artère carotide et un point de l'artère fémorale par mise en oeuvre d'un procédé tel que décrit ci-dessus, b) déterminer la vitesse de propagation brachiale d'une onde de pouls entre un point de l'artère carotide et un point de l'artère brachiale, par mise en oeuvre d'un procédé tel que décrit ci-dessus, c) mesurer la pression pulsée brachiale; d) estimer la pression aortique à partir des vitesses de propagation brachiale et aortique et de la pression pulsée brachiale. L'invention sera mieux comprise à la lecture de la description qui va suivre, donnée uniquement à titre d'exemple et faite en se référant aux des- sins, sur lesquels: - La figure 1 est une vue schématique de l'installation selon l'inven- tion, montrant les points d'implantation des capteurs sur un être humain; - la figure 2 est un organigramme explicitant les étapes du procédé selon l'invention; et - La figure 3 est une vue d'un enregistrement au cours du temps des mesures de pressions effectuées aux trois points de mesure du réseau arté- riel d'un être humain. L'installation représentée sur la figure 1 est adaptée pour le calcul, d'une part, de la pression pulsée aortique d'un être humain et, d'autre part, de la vitesse de propagation d'une onde de pouls entre deux points de me- sure distants du réseau artériel de l'être humain. A cet effet, elle met en ceuvre un procédé de détermination du déca- lage temporel entre les instants de passage d'une même onde de pouls en deux points de mesure distincts du réseau artériel de l'être humain. L'installation selon l'invention, représentée sur la figure 1, comporte une unité centrale de traitement d'informations 10, telle qu'un micro- ordinateur. Elle comporte de plus une carte d'acquisition 12 assurant la liai- son de trois capteurs de pression artérielle 14, 16, 18 à l'unité centrale de traitement d'informations 10. Elle comporte en outre une pédale pour arrêter l'acquisition et démarrer le traitement des données recueillies. Le premier capteur de pression 14 est destiné à être placé sur le cou du patient, au-dessus de l'artère carotide primitive au point M1. Le second capteur 16 est adapté pour être placé sur l'aine du patient au-dessus de l'artère fémorale au point M2. Le troisième capteur 18 est adapté pour être placé sur l'extrémité distale de l'avant-bras du patient, au-dessus de l'artère radiale au point M3. L'emplacement des capteurs peut être modifié selon le souhait de l'investigateur. Ainsi, il pourrait placer le capteur 14 en regard de l'artère bra- chiale et le capteur 18 en regard de l'artère radiale pour l'étude du segment artériel: brachial-radial. Ces trois capteurs sont tenus sur le corps du patient en regard de l'artère correspondante par un bracelet élastique ou tout autre système de- fixation adapté. L'unité centrale de traitement d'informations 10 est munie d'une base de temps constituée par exemple par son horloge interne. Elle met en ceuvre un programme adapté dont l'organigramme est représenté sur la figure 2. La première étape, notée 100, consiste à entrer dans l'unité centrale , la distance séparant les capteurs 14 et 16 d'une part, et 14 et 18 d'autre part. Ces distances sont notées dl2 et d13 respectivement. Elles sont mesurées par exemple à l'aide d'un mètre à ruban mainte- nu manuellement le long du corps du patient. A l'étape 102, l'unité centrale de traitement d'informations 10 procède, sur un intervalle de temps déterminé, à l'enregistrement simultané de l'évo- lution de la pression artérielle aux points de mesure Ml, M2, M3, à partir des capteurs 14, 16, 18 et de la carte d'acquisition 12. A cet effet, I'unité centrale enregistre, par exemple sur un disque dur, trois séries de mesures dis- crètes notées (an), (bn), (Cn) effectuées à des instants d'acquisition espacés d'un intervalle constant propre à la carte d'acquisition 12. Après acquisition des enregistrements (an), (bn) et (Cn), chaque enre- gistrement est lissé à l'étape 104 afin de filtrer les hautes fréquences typi- ques des mesures physique. A cet effet, l'algorithme suivant est utilisé pour chaque enregistrement. Il est présenté ici pour l'enregistrement (an). for i: = 2 à N-1 do moyenne:= (a[i-1]+a[i]+a[i+l])/3 dl = abs(a[i-1]moyenne); d2:= abs(a[i]-moyenne); d3:= abs(a[i+l]-moyenne); if dl construction de courbes telles qu'illustrées sur la figure 3. Ces courbes, no- tées respectivement A, B, C, représentent l'évolution de l'onde de pression au cours du temps aux points Ml, M2, M3. On constate sur ces courbes que l'onde de pouls passe à chacun des points de mesure avec un décalage temporel. A l'étape 106, I'unité centrale 10 détermine, à partir de l'analyse des enregistrements (an) et (bn) d'une part et (a,) et (c,) d'autre part, les décala- ges temporels entre les instants de passage d'une même onde de pouls aux points de mesure Ml et M2, d'une part, et Ml et M3 d'autre part. Ces décala- ges temporels sont notés t12, t13, respectivement. Le même algorithme est mis en ceuvre pour le calcul du décalage t12 et pour le calcul du décalage t13. Aussi, seul le procédé de mesure du déca- lage temporel t12 sera décrit en détail. Pour la détermination du décalage temporel t12, l'unité centrale 10 détermine d'abord pour le premier enregistrement (an) I'instant t1 d'appari- tion de l'onde de pouls en identifiant le pied de l'onde de pouls. Elle déter- mine ensuite pour le deuxième enregistrement (bn), I'instant t2 d'apparition de la même onde de pouls en identifiant, dans le deuxième enregistrement (bn), I'apparition de la partie correspondante du pied de l'onde de pouls. Elle calcule enfin le décalage temporel à partir des instants t1 et t2 déterminés sur les premier et deuxième enregistrements. L'instant ti d'apparition de l'onde de pouls sur le premier enregistre- ment (an) est déterminé avantageusement comme l'instant o la dérivée de la fonction enregistrée est maximale. L'instant t2 d'apparition de l'onde de pouls correspondante au point M2 est obtenu par recherche dans le deuxième enregistrement (bn) de l'instant o le coefficient de corrélation linéaire est maximal entre une plage du pre- mier enregistrement définie autour de l'instant t1 et une plage similaire cor- respondante du deuxième enregistrement. Cela permet de déterminer l'ins- tant o la dérivée de la fonction enregistrée est maximale. Plus précisément, la détermination du décalage temporel t12 est réali- sée par mise en oeuvre des étapes suivantes. On calcule d'abord en utilisant un schéma de différences finies d'ordre 1 la valeur maximale en valeur absolue de la "dérivée" de la première courbe établie à partir de l'enregistrement (an), i.e. la valeur C définie par C = {nelN tel que Vie[0,N-1] J ai+1 - ai éléments successifs est maximale. On définit ensuite une fenêtre de travail [mo, Mo] centrée en C en po- sant mo = C - Valeur Min, Mo = C + Valeur Max; On cherche alors la valeur Ke[O, N - Mo] qui maximise le coefficient de corrélation linéaire entre les courbes discrètes a = (an)m0o -------- -- -=max(--------------- c(a)a(bK) c(a)a(bk) On détermine enfin l'instant t2 correspondant à l'instant o a été ef- fectuée la mesure d'indice K dans le deuxième enregistrement. Le décalage t12 est alors calculé par la relation t12 = t2 - t,. A l'étape 108, l'unité centrale 10 calcule, les vitesses v12 et V13 de propagation de l'onde de pouls, d'une part, entre les points M1 et M2 et d'au- tre part, entre les points Ml et M3. La vitesse v12 correspond à la vitesse de propagation aortique encore désignée par vitesse de propagation tronculaire. La vitesse V13 correspond à la vitesse de propagation brachiale. Le calcul des vitesses est effectué à partir des distances d12 et d,3 mesurées, ainsi que des décalages temporels t12 et t13 déterminés à l'étape 106. L'onde de pouls étant issue du coeur, il convient de corriger les dis- tances d12 et d13 par la soustraction de la distance parcourue dans un sens opposé par l'une des deux ondes qui parcourt le trajet artériel dans un sens opposé. Ainsi, à l'étape 108, les vitesses de propagation de l'onde de pouls sont calculées par les relations suivantes d12 - d d13 d' V12 =----------- et v13= ------------ t12 t13 o d et d' représentent les longueurs du trajet artériel parcouru par l'onde dans des sens opposés après rencontre d'une bifurcation conduisant celle-ci vers chacun des points de mesure. Par exemple, la distance d mesurée en- tre la carotide et l'aorte est soustraite de la distance mesurée entre la caro- tide et la fémorale puisque l'onde entre l'aorte et la carotide se propage dans un sens opposé (ascendant) à celui entre l'aorte et la fémorale (descen- dant). Après calcul, I'unité centrale 10 met à disposition du praticien, no- tamment par affichage à l'écran, les valeurs v12 et v13. Aux étapes 110 et 112, I'unité centrale procède à une estimation de la pression pulsée aortique. Pour l'estimation de la pression pulsée aortique, le praticien possède d'abord manuellement à la mesure de la pression pulsée brachiale à l'aide d'un brassard gonflable et d'un stéthoscope appliqué de manière classique autour du bras du patient. La pression pulsée brachiale, encore appelée pression différentielle brachiale, est obtenue par différence de la pression systolique et de la pression diastolique mesurées par le praticien. Cette pression pulsée brachiale pourrait être aussi déterminée par des appareils automatiques de mesure de la pression artérielle (appareil électronique). A l'étape 110, la pression pulsée brachiale ainsi obtenue est entrée dans l'unité centrale de traitement d'informations 10. A l'étape 112, I'unité de traitement d'informations calcule la pression pulsée aortique estimée, notée Pa, à partir de la relation Pa = PbF Va/Vb o: Va = vitesse de propagation aortique Vb = vitesse de propagation brachiale, et Pb = pression pulsée brachiale. Les vitesses va et vb sont prises égales respectivement aux vitesses v12 et v13 calculées à l'étape 108. La pression pulsée aortique ainsi calculée donne une estimation sa- tisfaisante de la valeur réelle de la pression pulsée aortique. Cette dernière est difficilement mesurable puisque l'aorte est inaccessible depuis l'extérieur du corps du patient. Pour comprendre que l'estimation de la pression pulsée aortique est satisfaisante, il est possible de considérer le modèle d'une artère selon le- quel, en l'absence de réflexion, le rapport de l'amplitude de la pression proximale Pp à l'amplitude de la pression distale Pd aux deux extrémités d'une artère est proportionnel à la racine carrée des caractéristiques d'im- pédance de l'artère, aux deux extrémités d'une artère, d'o la relation Pp/Pd = ZZ o: Zp est l'impédance proximale de l'artère, et Zd est l'impédance distale de l'artère. En outre, d'après la formule de "Water Hammer", selon laquelle la vitesse d'une onde pulsée est égale au produit de l'impédance caractéristi- que d'une artère par la densité du sang, I'expression précédente devient: Pp/Pd = Vp/'Vd o: Vp est la vitesse de propagation proximale de l'onde de pouls dans le sys- tème artériel, et Vd est la vitesse de propagation distale de l'onde de pouls dans le système artériel. Ainsi, on comprend que du fait de la non uniformité de l'élasticité arté- rielle mesurée sur un être vivant, il est possible, à partir de la pression pul- sée brachiale, de déterminer une estimation satisfaisante de la pression pul- sée aortique. On constate qu'avec le procédé mis en oeuvre pour la détermination des décalages temporels t12 et t13, la précision des mesures est accrue. En effet, I'utilisation du pied du pic traduisant l'onde de pouls sur les courbes obtenues à partir des enregistrements effectués aux différents points de me- sure permet, en étant associée à une analyse de corrélation entre les diffé- rentes courbes, une estimation précise du décalage temporel. Dans le procédé décrit ici, la détection du passage d'une onde de pouls en un point donné du réseau artériel est réalisée par suivi de l'évolu-- tion de la pression en ce point. Toutefois, tout autre paramètre représentatif de l'onde de pouls peut être utilisé, tel que le diamètre de l'artère ou la vi- tesse du flux sanguin. Un capteur adapté est alors utilisé en chaque point de mesure à la place du capteur de pression. REVENDICATIONS 1.- Procédé de détermination du décalage temporel entre les instants de passage d'une même onde de pouls en deux points de mesure distincts (M1, M2) d'un réseau artériel d'un être vivant, comportant les étapes de: a) placer un capteur d'un paramètre représentatif de l'onde de pouls (14,16) à chaque point de mesure du réseau artériel; b) enregistrer l'évolution au cours du temps du paramètre représenta- tif de l'onde de pouls à chaque point de mesure; c) déduire des deux enregistrements ledit décalage temporel, caractérisé en ce que l'étape de déduction du décalage temporel comporte les étapes de: cl) déterminer, pour un premier enregistrement, un premier instant- (t,) d'apparition de l'onde de pouls en identifiant le pied de l'onde de pouls, c2) déterminer, pour le deuxième enregistrement, un deuxième ins- tant (t2) d'apparition de la même onde de pouls en identifiant, dans le deuxième enregistrement, I'apparition de la partie correspondante de l'onde de pouls, c3) calculer le décalage temporel à partir des instants t1 et t2 détermi- nés sur les premier et deuxième enregistrements. 2.- Procédé selon la revendication 1, caractérisé en ce que le premier instant (t,) d'apparition du pied de l'onde de pouls est déterminé comme l'instant o la dérivée de la fonction enregistrée est maximale. 3.- Procédé selon la revendication 1 ou 2, caractérisé en ce que le deuxième instant (t2) d'apparition de la même onde de pouls sur le deuxième enregistrement est obtenu par recherche, dans le deuxième enregistrement, de l'instant o un coefficient de corrélation est maximal entre une plage du premier enregistrement définie autour du premier instant (tl) et une plage correspondante du deuxième enregistrement. 4.- Procédé selon la revendication 3, caractérisé en ce que ledit coef- ficient de corrélation est un coefficient de corrélation linéaire. 5.- Procédé selon l'une quelconque des revendications précédentes, caractérisé en ce que ledit paramètre représentatif de l'onde de pouls est la pression artérielle. 6.- Procédé de détermination de la vitesse de propagation d'une onde de pouls entre deux points de mesure distants d'un réseau artériel d'un être vivant, comportant les étapes de: A) mesurer la distance séparant les deux points de mesure sur l'être vivant, B) déterminer le décalage temporel entre les instants de passage d'une même onde de pouls aux deux points de mesure; C) calculer la vitesse de propagation à partir de la distance séparant les deux points de mesure et du décalage temporel, caractérisé en ce que le décalage temporel entre les instants de passage de l'onde de pouls aux deux points de mesure est déterminé par un procédé selon l'une quelconque des revendications précédentes. 7.- Procédé d'estimation de la pression pulsée aortique, caractérisé en ce qu'il comporte les étapes de: a) déterminer la vitesse de propagation aortique d'une onde de pouls entre un point (Mi) de l'artère carotide et un point (M2) de l'artère fémorale par mise en oeuvre d'un procédé selon la revendication 6, b) déterminer la vitesse de propagation brachiale d'une onde de pouls entre un point (Mi) de l'artère carotide et un point (M3) de l'artère brachiale, par mise en oeuvre d'un procédé selon la revendication 6, c) mesurer la pression pulsée brachiale; d) estimer la pression aortique à partir des vitesses de propagation brachiale et aortique et de la pression pulsée brachiale. 8.- Procédé selon la revendication 7, caractérisé en ce que la pres- sion pulsée aortique Pa est déterminée par la relation: Pa = Pb 'I7/b o: Va est la vitesse de propagation aortique Vb est la vitesse de propagation brachiale, et Pb est la pression pulsée brachiale.