FR2826854A1

FR2826854A1 - Blood flow monitor estimates velocity and pressure from maximum derivative timing

Info

Publication number: FR2826854A1
Application number: FR0108891A
Authority: FR
Inventors: Roland Asmar
Original assignee: ARTECH MEDICAL
Current assignee: ARTECH MEDICAL
Priority date: 2001-07-04
Filing date: 2001-07-04
Publication date: 2003-01-10
Anticipated expiration: 2021-07-04
Also published as: FR2826854B1

Abstract

A blood flow monitor measures the time interval between passage of a derivative maximum in the blood pressure at the carotid (M1) and brachial (M2) arteries and calculates (10) the central pressure from the velocities.

Description

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é. The present invention relates to a method for determining the time shift between the instants of passage of the same pulse wave at two different measurement points of an arterial network of a living being. It also relates to a method for determining the speed of propagation of the pulse wave between these two measurement points, from the determined time offset.

Elle concerne enfin un procédé d'estimation de la pression artérielle pulsée centrale à partir de vitesses de propagation de l'onde de pouls déterminées. Finally, it relates to a method for estimating the central pulsed arterial pressure from determined pulse wave propagation speeds.

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émorale, au niveau de l'aine. Ce dernier capteur ou un capteur supplémentaire peut également être disposé sur l'artère brachiale ou radiale ou un autre site artériel à l'extrémité distale de l'avant-bras. An installation for measuring the speed of propagation of the pulse wave between two points of the body is currently known. This implements a method of determining the time shift between the instants of passage of the same pulse wave at two measurement points distinct from the arterial network of a human being. These measurement points are arranged for example on the primary carotid artery at the neck and on the femoral artery, at the level of the groin. The latter or additional sensor may also be disposed on the brachial or radial artery or other arterial site at the distal end of the forearm.

Pour déterminer la vitesse de propagation de l'onde de pouls, la distance séparant les deux points de mesure est mesurée manuellement à l'aide d'un mètre à ruban ou déterminé par tout autre procédé. To determine the speed of propagation of the pulse wave, the distance between the two measuring points is measured manually using a tape measure or determined by any other method.

A partir de la distance séparant les deux points et du décalage temporel 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 effectué 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. From the distance separating the two points from the determined time offset, the propagation speed is estimated by simply dividing the distance by the time offset, taking into account, however, the distance separating the core from one of the measurement points. or a possible path taken in one opposite direction by one of the two waves relative to the other, this distance to be deduced from the distance separating the two measuring points.

Pour la détermination du décalage temporel entre les instants de passage 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 In order to determine the time difference between the instants of passage of the same pulse wave at two distinct points of the arterial network, the known method consists in recording the evolution over time of the arterial pressure, at each point of measurement, then to deduce, by analysis

des deux enregistrements le décalage temporel entre les instants de passage de la même onde de pouls. of the two recordings, the time difference between the instants of passage of the same pulse wave.

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. For this purpose, the curves obtained from each record are analyzed. These curves present, during the passage of each pulse wave a characteristic peak representative of a heartbeat.

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, ou toute autre portion autour du pic et on détermine le décalage temporel entre ces deux maxima. Ce décalage peut être obtenu manuellement par mesure directe sur des courbes imprimées représentatives des deux enregistrements ou par traitement informatique des deux enregistrements. According to the current method, in order to determine the time shift, the peak pressure maxima corresponding to the same pulse wave, or any other portion around the peak, are considered and the time difference between these two maxima is determined. This offset can be obtained manually by direct measurement on printed curves representative of the two recordings or by computer processing of the two recordings.

Toutefois, la détermination des maxima des pics est relativement déli- cate et variable selon plusieurs paramètres cliniques et biologiques. Notamment, cette portion est affectée par les ondes de réflexion. 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. However, the determination of peak maxima is relatively sensitive and variable according to several clinical and biological parameters. In particular, this portion is affected by the reflection waves. Thus, the accuracy of determining the instants for calculating the time offset is low, which affects the overall quality of the measurement.

L'invention a pour but de fournir un procédé de précision accrue, notamment 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. The aim of the invention is to provide a method of increased precision, in particular by improving the determination of instants enabling the calculation of the time shift between the passages of the same pulse wave at two points distant from the arterial network of a living being.

A cet effet, l'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ésentatif 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 : To this end, the subject of the invention is a method for determining the time difference between the instants of passage of the same pulse wave at two distinct measurement points of an arterial network of a living being comprising the steps of: a) placing a sensor of a parameter representative of the pulse wave at each measurement point of the arterial network; b) record the evolution over time of the parameter representative of the pulse wave at each measurement point; c) deriving, from the two recordings, said time shift, characterized in that the step of deducing the time offset comprises the steps of:

c1) déterminer un premier instant où une dérivée d'un ordre prédéterminé par rapport au temps du premier enregistrement est maximale ; c2) sélectionner un intervalle de temps de référence contenant ledit premier instant ; c3) évaluer le décalage temporel entre ledit intervalle de temps de référence et un intervalle de temps mobile de même durée pour lequel une distance entre la plage de valeurs du premier enregistrement correspondant audit intervalle de temps de référence et la plage de valeurs du deuxième enregistrement correspondant audit intervalle de temps mobile est minimale, la distance étant définie par la somme des valeurs absolues des écarts point à point entre les plages des valeurs des premier et deuxième enregistrements ramenés à une même origine temporelle ; et c4) considérer le décalage temporel à déterminer comme égal au décalage temporel évalué. c1) determining a first time when a derivative of a predetermined order with respect to the time of the first record is maximum; c2) selecting a reference time interval containing said first instant; c3) evaluating the time difference between said reference time interval and a movable time interval of the same duration for which a distance between the value range of the first record corresponding to said reference time interval and the value range of the corresponding second record. said mobile time interval is minimal, the distance being defined by the sum of the absolute values of the point-to-point differences between the ranges of the values of the first and second records brought back to the same time origin; and c4) consider the time offset to be determined as equal to the evaluated time offset.

Suivant des modes particuliers de réalisation, le procédé comporte l'une ou plusieurs des caractéristiques suivantes : - ledit premier instant est déterminé comme l'instant où la dérivée première par rapport au temps du premier enregistrement est maximale ; - ledit intervalle de temps de référence est un intervalle s'achevant sensiblement audit premier instant ; - ledit premier instant est déterminé comme l'instant où la dérivée seconde par rapport au temps du premier enregistrement est maximale ; - ledit intervalle de temps de référence est tel que ledit premier instant se situe sensiblement à mi-durée dudit intervalle ; - ladite dérivée d'ordre prédéterminé est évaluée par application successive d'une méthode des différences finies d'ordre 1 ; - ledit paramètre représentatif de l'onde de pouls est la pression artérielle. According to particular embodiments, the method comprises one or more of the following characteristics: said first instant is determined as the time when the first derivative with respect to the time of the first recording is maximum; said reference time interval is an interval ending substantially at said first instant; said first instant is determined as the moment when the second derivative with respect to the time of the first recording is maximum; said reference time interval is such that said first instant is substantially mid-duration of said interval; said predetermined order derivative is evaluated by successive application of a finite difference method of order 1; said representative parameter of the pulse wave is the arterial pressure.

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, The invention further relates to a method for determining the speed of propagation of a pulse wave between two measurement points remote from an arterial network of a living being, comprising the steps of:
A) measuring the distance separating the two measuring points on the living being,

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é tel que défini ci-dessus. B) determining the time difference between the instants of passage of the same pulse wave at the two measurement points;
C) calculating the propagation velocity from the distance separating the two measuring points and the time offset, characterized in that the time difference between the instants of passage of the pulse wave at the two measurement points is determined by a process as defined above.

L'invention a enfin pour objet un procédé d'estimation de la pression pulsée centrale, caractérisé en ce qu'il comporte les étapes de : a) déterminer la vitesse de propagation centrale 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 dans le membre supérieur ; d) estimer la pression centrale à partir des vitesses de propagation brachiale et centrale et de la pression pulsée dans le membre supérieur. Finally, the subject of the invention is a method for estimating the central pulsating pressure, characterized in that it comprises the steps of: a) determining the central propagation speed of a pulse wave between a point of the artery carotid artery and a point of the femoral artery using a method as described above, b) determining the brachial velocity of propagation of a pulse wave between a point of the carotid artery and a point of the brachial artery, by carrying out a method as described above, c) measuring the pulsed pressure in the upper limb; d) estimate the central pressure from the brachial and central propagation velocities and the pulsed pressure in the upper limb.

L'invention concerne enfin une installation de diagnostic caractérisée en ce qu'elle comporte : - une base de données de valeurs de référence, - des moyens de détermination par mise en oeuvre d'un procédé tel que défini ci-dessus de l'une des grandeurs comprises dans le groupe consistant en le décalage temporel entre les instants de passage d'une même onde de pouls en deux points de mesure distants d'un réseau artériel d'un être vivant, 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 et la pression pulsée centrale d'un être vivant ; et - des moyens de comparaison de ladite grandeur déterminée avec au moins une valeur de référence. Finally, the invention relates to a diagnostic installation characterized in that it comprises: a database of reference values; means for determining by implementing a method as defined above of one of the following: magnitudes included in the group consisting of the temporal shift between the instants of passage of the same pulse wave at two measurement points distant from an arterial network of a living being, the speed of propagation of a pulse wave between two measuring points distant from an arterial network of a living being and the central pulsating pressure of a living being; and means for comparing said determined quantity with at least one reference value.

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 dessins, sur lesquels : The invention will be better understood on reading the description which follows, given solely by way of example and with reference to the drawings, in which:

- la figure 1 est une vue schématique de l'installation selon l'invention, montrant les points d'implantation des capteurs sur un être humain ; - la figure 2 est un organigramme explicitant les étapes d'un premier mode de mise en oeuvre du procédé selon l'invention ; - 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 ; - la figure 4 est une vue schématique agrandie d'un enregistrement au cours du temps de deux mesures de pression illustrant la mise en oeuvre de l'algorithme illustré sur la figure 2 ; - la figure 5 est un organigramme explicitant les étapes d'un second mode de mise en oeuvre du procédé selon l'invention ; et - la figure 6 est une vue schématique agrandie d'un enregistrement au cours du temps de deux mesures de pression effectuées en deux points illustrant la mise en oeuvre du procédé illustré sur la figure 5. FIG. 1 is a schematic view of the installation according to the invention, showing the locations of the sensors on a human being; FIG. 2 is a flow diagram explaining the steps of a first mode of implementation of the method according to the invention; FIG. 3 is a view of a recording over time of pressure measurements taken at the three measuring points of the arterial network of a human being; FIG. 4 is an enlarged schematic view of a recording over time of two pressure measurements illustrating the implementation of the algorithm illustrated in FIG. 2; FIG. 5 is a flowchart explaining the steps of a second embodiment of the method according to the invention; and FIG. 6 is an enlarged schematic view of a recording over time of two pressure measurements made at two points illustrating the implementation of the method illustrated in FIG. 5.

L'installation représentée sur la figure 1 est adaptée pour le calcul, d'une part, de la pression pulsée centrale d'un être humain et, d'autre part, de la vitesse de propagation d'une onde de pouls entre deux points de mesure distants du réseau artériel de l'être humain. The installation shown in FIG. 1 is suitable for calculating, on the one hand, the central pulsating pressure of a human being and, on the other hand, the speed of propagation of a pulse wave between two points. remote measurement of the arterial network of the human being.

A cet effet, elle 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 de l'être humain. For this purpose, it implements a method of determining the time shift between the instants of passage of the same pulse wave at two measurement points distinct from the arterial network of the human being.

L'installation selon l'invention, représentée sur la figure 1, comporte une unité centrale de traitement d'informations 10, telle qu'un microordinateur. Elle comporte de plus une carte d'acquisition 12 assurant la liaison à l'unité centrale de traitement d'informations 10 de plusieurs capteurs, notamment trois, pour la mesure par exemple de la pression artérielle 14, 16,18. Elle comporte en outre une pédale pour arrêter l'acquisition et démarrer le traitement des données recueillies. The installation according to the invention, represented in FIG. 1, comprises a central information processing unit 10, such as a microcomputer. It furthermore comprises an acquisition card 12 connecting the central information processing unit 10 with several sensors, in particular three sensors, for measuring, for example, arterial pressure 14, 16, 18. It also includes a pedal to stop the acquisition and start processing the data collected.

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. The first pressure sensor 14 is intended to be placed on the patient's neck, above the primary carotid artery at 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. The second sensor 16 is adapted to be placed on the groin of the patient above the femoral artery at 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. The third sensor 18 is adapted to be placed on the distal end of the patient's forearm, above the radial artery at the M3 point.

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 brachiale au point M4 et le capteur 18 en regard de l'artère radiale pour l'étude du segment artériel : brachial-radial ou encore fémoro-tibiale. The location of the sensors can be changed according to the wishes of the investigator. Thus, it could place the sensor 14 opposite the brachial artery at point M4 and the sensor 18 facing the radial artery for the study of the arterial segment: brachial-radial or femoro-tibial.

Les capteurs peuvent également être placés à tout autre emplacement tel qu'au point M5 au dessus de l'artère tibiale ou au point M6 au dessus de l'artère pédieuse. The sensors may also be placed at any other location such as at point M5 above the tibial artery or at point M6 above the pedicle artery.

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é. These three sensors are held on the body of the patient facing the corresponding artery by an elastic strap or other suitable fastening system.

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 oeuvre un programme adapté dont l'organigramme d'un premier exemple de réalisation est illustré sur la figure 2. The central information processing unit 10 is provided with a time base constituted for example by its internal clock. It implements an adapted program whose organization chart of a first exemplary embodiment is illustrated in FIG. 2.

La première étape, notée 100, consiste à entrer dans l'unité centrale 10, la distance séparant les capteurs 14 et 16 d'une part, et 14 et 18 d'autre part. Ces distances sont notées dis et d13 respectivement. The first step, denoted 100, consists in entering the central unit 10, the distance separating the sensors 14 and 16 on the one hand, and 14 and 18 on the other hand. These distances are denoted dis and d13 respectively.

Elles sont mesurées par exemple à l'aide d'un mètre à ruban maintenu manuellement le long du corps du patient. Elle peut être mesurée par tout autre procédé. They are measured for example using a tape measure held manually along the body of the patient. It can be measured by any other method.

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'évolution de la pression artérielle aux points de mesure Mi, M2, M3, à partir des capteurs 14,16, 18 et de la carte d'acquisition 12. A cet effet, l'unité centrale 10 enregistre, par exemple sur un disque dur, trois séries de mesures discrètes notées (an), (bn), (cn) effectuées à des instants d'acquisition espacés d'un intervalle d'échantillonnage constant propre à la carte d'acquisition 12. In step 102, the central information processing unit 10 proceeds, on a determined time interval, to the simultaneous recording of the evolution of the arterial pressure at the measurement points Mi, M2, M3, from sensors 14, 16, 18 and the acquisition card 12. For this purpose, the central unit 10 records, for example on a hard disk, three series of discrete measurements noted (an), (bn), (cn ) performed at acquisition times spaced a constant sampling interval specific to the acquisition card 12.

La fréquence d'échantillonnage est par exemple d'environ 800 Hz pour chaque série de mesure. La résolution est par exemple de 0,025 %, chaque valeur de mesure étant codée sur 12 bits. Chaque enregistrement est supposé contenir N valeurs. The sampling frequency is for example about 800 Hz for each measurement series. The resolution is, for example, 0.025%, each measurement value being coded on 12 bits. Each record is assumed to contain N values.

Après acquisition des enregistrements (an), (bn) et (cn), chaque enregistrement est lissé à l'étape 104 afin de filtrer les hautes fréquences typiques 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 : =OàN-1do moyenne : = (a [i-1] +a [i] +a [i+1])/3 ; d1 : = abs (a [i-1]-moyenne) ; d2 : = abs (a [i]-moyenne) ; d3 : = abs (a[i+1]-moyenne) ; ifd1 < d2thena [i] : = a[i-1] ; fi ; if d3 < d2 then a [i] : = a[i+1] ; fi ; od ;
Après filtrage, les séries de mesures (an), (bn), (Cn) permettent la construction de courbes telles qu'illustrées sur la figure 3. Ces courbes, notées respectivement A, B, C, représentent l'évolution de l'onde de pression au cours du temps aux points Mi, M2, M3. After acquiring records (an), (bn) and (cn), each record is smoothed in step 104 to filter the high frequencies typical of the physical measurements. For this purpose, the following algorithm is used to

each record. It is presented here for recording (year).

for i: = O-N-1do mean: = (a [i-1] + a [i] + a [i + 1]) / 3; d1: = abs (a [i-1] -middle); d2: = abs (a [i] -middle); d3: = abs (a [i + 1] -middle); ifd1 <d2thena [i]: = a [i-1]; fi; if d3 <d2 then a [i]: = a [i + 1]; fi; od;
After filtering, the series of measurements (an), (bn), (Cn) allow the construction of curves as illustrated in Figure 3. These curves, denoted respectively A, B, C, represent the evolution of the pressure wave over time at points Mi, M2, M3.

On constate sur ces courbes que l'onde de pouls passe à chacun des points de mesure avec un décalage temporel. These curves show that the pulse wave passes at each of the measurement points with a time shift.

A partir de l'étape 106, t'unité centrale 10 détermine, à partir de l'analyse des enregistrements (an) et (bn), d'une part, et (an) et (en), d'autre part, les décalages temporels entre les instants de passage d'une même onde de pouls aux points de mesure Mi et M2, d'une part, et M1 et M3 d'autre part. From step 106, the central unit 10 determines, from the analysis of the records (an) and (bn), on the one hand, and (an) and (en), on the other hand, the time offs between the instants of passage of the same pulse wave at the measurement points Mi and M2, on the one hand, and M1 and M3, on the other hand.

Ces décalages temporels sont notés tis, tis, respectivement. These time offsets are denoted tis, tis, respectively.

Le même algorithme est mis en oeuvre pour le calcul du décalage t12 et pour le calcul du décalage tig. The same algorithm is implemented for calculating the shift t12 and for calculating the shift tig.

Pour la détermination du décalage temporel t12, t'unité centrale 10 détermine d'abord pour le premier enregistrement (an) une plage de valeurs d'apposition de l'onde de pouls. Elle calcule enfin le décalage temporel entre la plage de valeurs identifiée dans le premier enregistrement et la plage de valeurs correspondante dans le deuxième enregistrement. Ce calcul se fait en recherchant le minimum d'une distance d définie entre des plages de valeurs des premier et deuxième enregistrements. For the determination of the time shift t12, the central unit 10 first determines for the first record (an) a range of values of apposition of the pulse wave. Finally, it calculates the time difference between the range of values identified in the first record and the corresponding range of values in the second record. This calculation is done by searching for the minimum of a distance d defined between ranges of values of the first and second records.

A l'étape 106, t'unité centrale 10 détermine d'abord, pour chacun des enregistrements (an), (bn) et (cn), la dérivée première par rapport au temps. In step 106, the central unit 10 first determines, for each of the records (an), (bn) and (cn), the first derivative with respect to time.

Celle-ci est notée respectivement (a'n), (b'n) et (c'n). This is noted respectively (a'n), (b'n) and (c'n).

Les valeurs de la dérivée première par rapport au temps sont déterminées par application d'une méthode des différences finies d'ordre 1. The values of the first derivative with respect to time are determined by applying a finite difference method of order 1.

Ainsi, dans la mesure où les intervalles d'échantillonnage sont égaux, on a par exemple :

Thus, insofar as the sampling intervals are equal, for example:

A l'étape 110, J'unité centrale 10 détermine l'indice de la valeur maximale de chacune des dérivées premières des enregistrements (an), (bn) et (cn). Les indices des valeurs pour lesquelles les dérivées premières (a'n), (b'n) et (c'n) sont maximales sont notés respectivement p, q et r pour les enregistrements (an), (bn) et (cn).

In step 110, the central unit 10 determines the index of the maximum value of each of the first derivatives of the records (an), (bn) and (cn). The indices of the values for which the first derivatives (a'n), (b'n) and (c'n) are maximal are respectively denoted p, q and r for the records (an), (bn) and (cn) .

Ces indices sont donc définis par les relations suivantes :

p = {ne [0, N-1] tel que ViE[O, N-1] a') a'n} q = {ne [0, N-1] tel que V, < = [0, N-1] b', < b'n} r = {ne [0, N-1] tel que V, e [0, N-1] c', < c'n}

A l'étape 112, le nombre de périodes de décalage entre deux plages des deux enregistrements (an) et (bn) correspondant au passage d'une même onde de pouls est déterminé. These indices are therefore defined by the following relations:

p = {ne [0, N-1] such that ViE [O, N-1] a ') a'n} q = {ne [0, N-1] such that V, <= [0, N- 1] b ', <b'n} r = {ne [0, N-1] such that V, e [0, N-1] c', <c'n}

In step 112, the number of offset periods between two ranges of the two records (an) and (bn) corresponding to the passage of a same pulse wave is determined.

A cet effet, un premier intervalle de temps de référence, noté IR, est sélectionné. Cet intervalle de temps est choisi de sorte que l'instant noté tp où la dérivée (a'n) est maximale soit contenu dans cet intervalle de temps de référence IR. For this purpose, a first reference time interval, denoted IR, is selected. This time interval is chosen so that the moment noted tp where the derivative (a'n) is maximum is contained in this reference time interval IR.

Avantageusement, l'intervalle de temps de référence IR est choisi de telle sorte que sa durée soit inférieure à la durée séparant deux battements cardiaques. En outre, l'intervalle de temps de référence est tel que l'instant tp où la dérivée première (a'n) est maximale constitue la borne supérieure de cet intervalle de temps, c'est-à-dire que cet intervalle de temps s'achève sensiblement à cet instant tp. Advantageously, the IR reference time interval is chosen such that its duration is less than the duration between two heart beats. In addition, the reference time interval is such that the time tp where the first derivative (a'n) is maximum constitutes the upper bound of this time interval, that is to say that this time interval ends substantially at this instant tp.

Dans la pratique, cet intervalle de temps de référence IR a une durée égale à un nombre entier de périodes d'échantillonnage, cette période étant notée T. Le nombre de périodes est compris entre 20 et 40 de préférence. In practice, this IR reference time interval has a duration equal to an integer number of sampling periods, this period being noted T. The number of periods is preferably between 20 and 40.

Dans l'exemple considéré, la durée de l'intervalle de référence IR est égale à 30 T, c'est-à-dire environ 37,5 millisecondes. In the example considered, the duration of the reference interval IR is equal to 30 T, that is to say approximately 37.5 milliseconds.

La sélection de l'intervalle de temps de référence IR permet de déterminer une plage de valeurs du premier enregistrement (an), cette plage étant constituée des valeurs du premier enregistrement (an) pour chacun des instants d'échantillonnage compris dans l'intervalle de référence IR. The selection of the IR reference time interval makes it possible to determine a range of values of the first record (an), this range consisting of the values of the first record (an) for each of the sampling times in the range of IR reference.

L'instant tp et l'intervalle de référence IR sont illustrés sur la figure 4. The instant tp and the reference interval IR are illustrated in FIG. 4.

Afin de déterminer le décalage temporel entre les enregistrement (an) et (bn), un ensemble de plages de valeurs du deuxième enregistrement (bn) est défini. Les plages de valeurs du deuxième enregistrement constituant cet ensemble correspondent à des intervalles de temps mobiles IM ayant la même durée que l'intervalle de temps de référence IR et décalés de celui-ci d'un nombre entier de périodes d'échantillonnage T. In order to determine the time offset between the records (an) and (bn), a set of ranges of values of the second record (bn) is defined. The ranges of values of the second record constituting this set correspond to mobile time slots IM having the same duration as the reference time interval IR and shifted therefrom by an integer number of sampling periods T.

En particulier, les intervalles mobiles considérés IM sont constitués de tous les intervalles de temps décalés d'un nombre fini de périodes d'échantillonnage T depuis l'intervalle de référence IR, les durées de décalage étant comprises entre 0 et une valeur maximale égale à q-p + 50, p et q étant les indices où les dérivées premières des premier et deuxième enregistrements (an) et (bn) sont maximales. L'instant tq est illustré sur la figure 4 et correspond à l'instant où la dérivée première (b'n) du deuxième enregistrement est maximale. In particular, the mobile intervals considered IM consist of all the time intervals shifted by a finite number of sampling periods T from the reference interval IR, the offset times being between 0 and a maximum value equal to qp + 50, where p and q are the indices where the first derivatives of the first and second records (an) and (bn) are maximal. The instant tq is illustrated in FIG. 4 and corresponds to the moment when the first derivative (b'n) of the second record is maximum.

Dans chacune des plages de valeurs correspondant aux intervalles mobiles lM, la distance d entre la plage du premier enregistrement (an) correspondant à l'intervalle de référence IR et chacune des plages de va- leurs du deuxième enregistrement (bn) correspondant aux intervalles mobi- les IM est calculée. In each of the ranges of values corresponding to the moving intervals lM, the distance d between the range of the first record (an) corresponding to the reference interval IR and each of the ranges of values of the second record (bn) corresponding to the intervals mobi - IM is calculated.

Cette distance d est définie par la somme des valeurs absolues des écarts point à point entre les plages de valeurs des premier et deuxième enregistrements ramenés à une même origine temporelle. This distance d is defined by the sum of the absolute values of the point-to-point differences between the ranges of values of the first and second records brought back to the same time origin.

Plus précisément, cette distance d est définie pour deux enregistrements (Xn), (yn) par :

More precisely, this distance d is defined for two records (Xn), (yn) by:

On détermine alors, à l'étape 112, indice i appartenant à [0, q-p + 50] pour lequel la distance entre les plages de valeurs correspondantes des

premier et deuxième enregistrements (an), (bn) est minimale.

En particulier, on recherche l'indice i tel que L 1 an - bn+, 1 ne [p-30, p] soit minimale. Then, in step 112, index i belonging to [0, qp + 50] is determined for which the distance between the corresponding ranges of values of

first and second records (year), (bn) is minimal.

In particular, we search the index i such that L 1 year - bn +, 1 not [p-30, p] is minimal.

Un calcul analogue est effectué pour déterminer un indice j correspondant au décalage pour lequel la distance entre des plages correspondantes des premier et troisième enregistrements (an), (cn) est minimale. A similar calculation is made to determine an index j corresponding to the offset for which the distance between corresponding ranges of the first and third records (an), (cn) is minimal.

A l'étape 114, les décalages temporels t12 et t13 sont calculés par multiplication des indices i et j trouvés précédemment par la durée de la période d'échantillonnage T.

In step 114, the time offsets t12 and t13 are calculated by multiplying the indices i and j found previously by the duration of the sampling period T.

Ainsi, on a : t12 = ixT, et t13 = jxT. Thus, we have: t12 = ixT, and t13 = jxT.

A l'étape 130, 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'autre part, entre les points M1 et M3. In step 130, the central unit 10 calculates the velocity V12 and V13 for propagation of the pulse wave, on the one hand, between the points M1 and M2 and, on the other hand, between the points M1 and M3. .

La vitesse V12 correspond à la vitesse de propagation centrale encore désignée par vitesse de propagation tronculaire. La vitesse V13 correspond à la vitesse de propagation brachiale. The velocity V12 corresponds to the central propagation velocity further designated by truncular propagation velocity. The speed V13 corresponds to the brachial propagation velocity.

Le calcul des vitesses est effectué à partir des distances d12 et d13 mesurées, ainsi que des décalages temporels t12 et t13 déterminés à l'étape 112. The velocities are calculated from the distances d12 and d13 measured, as well as time offsets t12 and t13 determined in step 112.

L'onde de pouls étant issue du coeur, il convient de corriger les distances 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é. The pulse wave coming from the heart, it is appropriate to correct the distances d12 and d13 by subtracting the distance traveled in an opposite direction by one of the two waves that travels the arterial path in an opposite direction.

Ainsi, à l'étape 130, les vitesses de propagation de l'onde de pouls sont calculées par les relations suivantes :

di2-d dl3-d Vl2=----etVi3=---tl2 tl3

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 entre 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 (descendant). Thus, in step 130, the propagation speeds of the pulse wave are calculated by the following relationships:

di2-d dl3-d Vl2 = ---- andVi3 = --- tl2 tl3

where d and d represent the lengths of the arterial path traveled by the wave in opposite directions after encountering a bifurcation leading it to each of the measurement points. For example, the distance d measured between the carotid and the aorta is subtracted from the distance measured between the carotid and the femoral since the wave between the aorta and the carotid is propagated in an opposite (ascending) direction to that between the aorta and the femoral (descending).

Après calcul, l'unité centrale 10 met à disposition du praticien, notamment par affichage à l'écran, les valeurs V12 et Vis. After calculation, the central unit 10 makes available to the practitioner, in particular by display on the screen, the values V12 and Vis.

Aux étapes 132 et 134, l'unité centrale procède à une estimation de la pression pulsée centrale. In steps 132 and 134, the central unit makes an estimate of the central pulsating pressure.

Pour l'estimation de la pression pulsée centrale, 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). For the estimation of the central pulsating pressure, the practitioner first manually measures the pressure of the brachial pressure using an inflatable cuff and a stethoscope conventionally applied around the patient's arm. Brachial pulse pressure, also called brachial pressure differential, is obtained by the difference between the systolic pressure and the diastolic pressure measured by the practitioner. This brachial pressure could also be determined by automatic blood pressure measuring devices (electronic device).

A l'étape 132, la pression pulsée brachiale ainsi obtenue est entrée dans l'unité centrale de traitement d'informations 10. In step 132, the brachial pulse pressure thus obtained is inputted to the central information processing unit 10.

A l'étape 134, l'unité de traitement d'informations calcule la pression pulsée centrale estimée, notée Pa, à partir de la relation

où : Va = vitesse de propagation centrale Vb = vitesse de propagation brachiale, et Pb = pression pulsée brachiale. In step 134, the information processing unit calculates the estimated central pulsating pressure, denoted Pa, from the relation

where: Va = central propagation velocity Vb = brachial propagation velocity, and Pb = brachial pulsed pressure.

La pression puisée brachiale peut-être remplacée par toute mesure de pression dans le membre supérieur telle que la pression radiale. The brachial pressure can be replaced by any pressure measurement in the upper limb such as radial pressure.

Les vitesses Va et Vb sont prises égales respectivement aux vitesses V12 et V13 calculées à l'étape 130. The speeds Va and Vb are taken equal to the speeds V12 and V13 respectively calculated in step 130.

La pression pulsée centrale ainsi calculée donne une estimation satisfaisante de la valeur réelle de la pression pulsée centrale. Cette dernière est difficilement mesurable puisque l'aorte est inaccessible depuis l'extérieur du corps du patient. The central pulsed pressure thus calculated gives a satisfactory estimate of the actual value of the central pulsating pressure. The latter is difficult to measure since the aorta is inaccessible from outside the body of the patient.

Pour comprendre que l'estimation de la pression pulsée centrale est satisfaisante, il est possible de considérer le modèle d'une artère selon lequel, 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'impédance de l'artère, aux deux extrémités d'une artère, d'où la relation :

où : Zp est l'impédance proximale de l'artère, et Zd est l'impédance distale de l'artère. To understand that the estimate of the central pulsating pressure is satisfactory, it is possible to consider the model of an artery according to which, in the absence of reflection, the ratio of the amplitude of the proximal pressure Pp to the amplitude Pd distal pressure at both ends of an artery is proportional to the square root of the impedance characteristics of the artery, at both ends of an artery, hence the relation:

where: Zp is the proximal impedance of the artery, and Zd is the distal impedance of the artery.

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éristique d'une artère par la densité du sang, l'expression précédente devient :

où : Vp est la vitesse de propagation proximale de l'onde de pouls dans le système artériel, et Vd est la vitesse de propagation distale de l'onde de pouls dans le système artériel. Moreover, according to the formula of "Water Hammer", according to which the speed of a pulsed wave is equal to the product of the characteristic impedance of an artery by the density of the blood, the preceding expression becomes:

where: Vp is the velocity of proximal propagation of the pulse wave in the arterial system, and Vd is the rate of distal propagation of the pulse wave in the arterial system.

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 puisée brachiale, de déterminer une estimation satisfaisante de la pression pulsée centrale. Thus, it is understood that because of the nonuniformity of the arterial elasticity measured on a living being, it is possible, from the brachial pulse pressure, to determine a satisfactory estimate of the central pulsating pressure.

A l'étape 136, t'un des décalages calculés à l'étape 114, l'une des vitesses de propagation de l'onde de pouls calculées à l'étape 130 où la pression pulsée centrale calculée à l'étape 134 est comparé à une valeur de référence standard contenue dans une base de données. Cette base de données comporte des valeurs de référence obtenues chez des sujets indemnes de maladie cardio-vasculaire. In step 136, one of the offsets calculated in step 114, one of the pulse wave propagation velocities calculated in step 130 where the central pulsating pressure calculated in step 134 is compared to a standard reference value contained in a database. This database contains reference values obtained in subjects free from cardiovascular disease.

Les valeurs de références peuvent également être obtenues par application de formules mathématiques à partir de données de référence mémorisées. Reference values can also be obtained by applying mathematical formulas from stored reference data.

En variante, un index de rigidité artérielle est déterminé à partir d'une vitesse de propagation de l'onde de pouls déterminée à l'étape 130 et la comparaison est effectuée à l'étape 136 entre cet index de rigidité artérielle calculé et un ensemble d'index de rigidité artérielle de référence. Alternatively, an arterial stiffness index is determined from a pulse wave velocity determined in step 130 and the comparison is made at step 136 between this computed arterial stiffness index and a set index of arterial stiffness of reference.

La valeur de la vitesse de l'onde de pouls et la rigidité artérielle, notamment en étant comparée à des valeurs de référence peuvent permettre à un praticien de formuler un pronostic cardio-vasculaire, ces données pouvant être utilisées pour le dépistage précoce des patients à haut risque cardio-vasculaire. The value of pulse wave velocity and arterial stiffness, especially when compared to reference values, may allow a practitioner to formulate a cardiovascular prognosis, which data can be used for the early detection of patients with high cardiovascular risk.

Le procédé mis en oeuvre pour la détermination des décalages tem- porels t12 et tis, permet de déterminer le passage d'une même partie caractéristique de l'onde de pouls sur les courbes obtenues à partir des enregistrements effectués aux différents points de mesure. La détermination des maxima des dérivées premières des enregistrements en cause et le calcul des minima de distances définies entre les enregistrements sont très simples à mettre en oeuvre. Ces calculs reposent en effet seulement sur des opérations mathématiques élémentaires. The method implemented for the determination of the temporal shifts t12 and tis makes it possible to determine the passage of the same characteristic part of the pulse wave on the curves obtained from the recordings made at the different measurement points. The determination of the maximums of the first derivatives of the recordings in question and the calculation of the minimum distances defined between the recordings are very simple to implement. These calculations are based only on elementary mathematical operations.

Ainsi, le procédé peut être mis en oeuvre sur des moyens de calcul de vitesse et de capacité réduite, permettant ainsi d'obtenir une installation de faible coût ayant une précision satisfaisante. Thus, the method can be implemented on speed calculation means and reduced capacity, thus obtaining a low cost installation having a satisfactory accuracy.

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'évolution 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- In the method described here, the detection of the passage of a pulse wave at a given point of the arterial network is performed by monitoring the evolution of the pressure at this point. However, any other parameter representative of the pulse wave can be used, such as the diameter of the artery or the vi-

tesse du flux sanguin. Un capteur adapté est alors utilisé en chaque point de mesure à la place du capteur de pression. blood flow. A suitable sensor is then used at each measuring point in place of the pressure sensor.

Un autre mode de mise en oeuvre du procédé selon l'invention est illustré sur les figures 5 et 6. Another embodiment of the method according to the invention is illustrated in FIGS. 5 and 6.

Dans cette variante, les étapes 100 à 106 sont identiques aux étapes ayant le même numéro dans l'algorithme illustré sur la figure 2. In this variant, the steps 100 to 106 are identical to the steps having the same number in the algorithm illustrated in FIG.

Toutefois, à l'issue de l'étape 106, une étape 208 est mise en oeuvre pour le calcul de la dérivée seconde par rapport au temps de chacun des enregistrements (an), (bn) et (Cn). However, at the end of step 106, a step 208 is implemented for calculating the second derivative with respect to the time of each of the records (an), (bn) and (Cn).

A l'étape 208, les dérivées secondes par rapport au temps de chaque enregistrement (an), (bn) et (cn) sont calculées à partir des dérivées premières (a'n), (b'n) et (c'n) par application d'une méthode de différences finies d'ordre 1. Ces dérivées secondes sont notées respectivement (a"n), (b"n) et (C"n). In step 208, the second derivatives with respect to the time of each record (an), (bn) and (cn) are calculated from the first derivatives (a'n), (b'n) and (c'n) ) by applying a method of finite differences of order 1. These second derivatives are denoted respectively (a "n), (b" n) and (C "n).

A cet effet, les valeurs des dérivées secondes par rapport au temps sont données par les relations :

For this purpose, the values of the second derivatives with respect to time are given by the relations:

A l'étape 210, l'unité centrale 10 détermine l'indice de la valeur maximale de chacune des dérivées secondes des enregistrements (an), (bn) et (cn). Les indices des valeurs pour lesquelles les dérivées secondes (a"n), (b"n), (ciln) sont maximales sont notés respectivement s, t et u pour les enregistrements (an), (bn) et (cn).

In step 210, the CPU 10 determines the index of the maximum value of each of the second derivatives of the records (an), (bn) and (cn). The indices of the values for which the second derivatives (a "n), (b" n), (ciln) are maximal are denoted respectively s, t and u for the records (an), (bn) and (cn).

Ceux-ci sont définis par les relations suivantes :

s= {ne [0, N-1] tetqueV, e [0, N-1] a", < a"n} t = {n < = [0, N-1] tel que Vj&num; [0, N-1] b", < b"n} u = {nr= [O, N-1] tel que Vie [O, N-1] C"i : C"nl.

These are defined by the following relationships:

s = {ne [0, N-1] tetqueV, e [0, N-1] a ", <a" n} t = {n <= [0, N-1] such that Vj &num; [0, N-1] b ", <b"n> u = {nr = [O, N-1] such that Vie [O, N-1] C "i: C" nl.

A l'étape 212, le nombre de périodes de décalage entre deux plages correspondantes de deux enregistrements (an) et (bn) correspondant au passage d'une même onde de pouls est déterminé. In step 212, the number of offset periods between two corresponding ranges of two records (an) and (bn) corresponding to the passage of a same pulse wave is determined.

A cet effet, un premier intervalle de temps de référence, noté KR, est sélectionné. Cet intervalle de temps est choisi de sorte que l'instant noté ts For this purpose, a first reference time interval, denoted KR, is selected. This time interval is chosen so that the instant noted ts

où la dérivée seconde (a"n) est maximale soit contenu dans cet intervalle de temps de référence. where the second derivative (a "n) is maximum is contained in this reference time interval.

Avantageusement, l'intervalle de temps de référence KR est choisi de telle sorte que sa durée soit inférieure à la durée séparant deux battements cardiaques et qu'il soit de préférence sensiblement égal à un quart de cette durée. En outre, l'intervalle de temps de référence KR est tel que l'instant ts où la dérivée seconde (a"n) est maximale se situe sensiblement à mi-durée de cet intervalle de temps. Advantageously, the reference time interval KR is chosen such that its duration is less than the duration separating two heart beats and that it is preferably substantially equal to a quarter of this duration. In addition, the reference time interval KR is such that the time ts where the second derivative (a "n) is maximum is substantially mid-duration of this time interval.

Dans la pratique, cet intervalle de temps de référence KR a une durée égale à un nombre entier de périodes d'échantillonnage, cette période étant notée T. In practice, this reference time interval KR has a duration equal to an integer number of sampling periods, this period being noted T.

Dans l'exemple considéré, la durée de l'intervalle de référence KR est égale à 30 T, c'est-à-dire environ 37,5 millisecondes. L'intervalle de référence KR s'étend donc de ts-15T à ts + 15T. In the example considered, the duration of the reference interval KR is equal to 30 T, that is to say approximately 37.5 milliseconds. The reference interval KR therefore ranges from ts-15T to ts + 15T.

La sélection de l'intervalle de temps de référence KR permet de déterminer une plage de valeurs du premier enregistrement (an), cette plage étant constituée des valeurs du premier enregistrement (an) pour chacun des instants d'échantillonnage compris dans l'intervalle de référence KR. The selection of the reference time interval KR makes it possible to determine a range of values of the first record (an), this range consisting of the values of the first record (an) for each of the sampling times included in the interval of reference KR.

L'instant ts et l'intervalle de référence KR sont illustrés sur la figure 6. The instant ts and the reference interval KR are illustrated in FIG.

Afin de déterminer le décalage temporel entre les enregistrement (an) et (bn), un ensemble de plages de valeurs du deuxième enregistrement (bn) est défini. Les plages de valeurs du deuxième enregistrement constituant cet ensemble correspondent à des intervalles de temps mobiles KM ayant la même durée que l'intervalle de temps de référence KR et décalés de celui-ci d'un nombre entier de périodes d'échantillonnage T. In order to determine the time offset between the records (an) and (bn), a set of ranges of values of the second record (bn) is defined. The ranges of values of the second record constituting this set correspond to mobile time slots KM having the same duration as the reference time interval KR and shifted therefrom by an integer number of sampling periods T.

En particulier, les intervalles mobiles considérés KM sont constitués de tous les intervalles de temps décalés d'un nombre fini de périodes d'échantillonnage T depuis l'intervalle de référence kR, les durées de décalage étant comprises entre 0 et une valeur maximale égale à t-s + 50, s et t étant les indices où les dérivées secondes des premier et deuxième enregistrements (an) et (bn) sont maximales. L'instant ts est illustré sur la figure 6 et correspond à l'instant où la dérivée seconde (b"n) du deuxième enregistrement est maximale. In particular, the mobile intervals considered KM consist of all the time intervals shifted by a finite number of sampling periods T from the reference interval kR, the offset times being between 0 and a maximum value equal to ts + 50, where s and t are the indices where the second derivatives of the first and second records (an) and (bn) are maximum. The instant ts is illustrated in FIG. 6 and corresponds to the moment when the second derivative (b "n) of the second record is maximum.

Dans chacune des plages de valeurs correspondant aux intervalles mobiles KM, la distance entre la plage du premier enregistrement (an) correspondant à l'intervalle de référence KR et chacune des plages de valeurs du deuxième enregistrement (bn) correspondant aux intervalles mobile IM est calculée. In each of the ranges of values corresponding to the mobile intervals KM, the distance between the range of the first record (an) corresponding to the reference interval KR and each of the ranges of values of the second record (bn) corresponding to the mobile intervals IM is calculated. .

Cette distance notée d est définie par la somme des valeurs absolues des écarts point à point entre les plages de valeurs des premier et deuxième enregistrements ramenés à une même origine temporelle. This distance denoted d is defined by the sum of the absolute values of the point-to-point differences between the ranges of values of the first and second records brought back to the same time origin.

More precisely, this distance d is defined for two records (Xn), (yn) by:

On détermine alors, à l'étape 112, l'indice k appartenant à [0, t-s + 50] pour lequel la distance entre les plages de valeurs correspondantes des

premier et deuxième enregistrements (an), (bn) est minimale.

En particulier, on recherche l'indice k tel que L 1 an - bn+, 1 ne [s- 15, s+ 15] soit minimale. Then, in step 112, the index k belonging to [0, ts + 50] is determined, for which the distance between the corresponding value ranges of

first and second records (year), (bn) is minimal.

In particular, we look for the index k such that L 1 year - bn +, 1 ne [s-15, s + 15] is minimal.

Un calcul analogue est effectué pour déterminer un indice 1 correspondant au décalage pour lequel la distance entre des plages correspondantes des premier et troisième enregistrements (an), (cn) est minimale. A similar calculation is performed to determine an index 1 corresponding to the offset for which the distance between corresponding ranges of the first and third records (an), (cn) is minimal.

A l'étape 214, les décalages temporels t12 et t13 sont calculés par multiplication des indices k et 1 trouvés précédemment par la durée de la période d'échantillonnage T.

In step 214, the time offsets t12 and t13 are calculated by multiplying the indices k and 1 found previously by the duration of the sampling period T.

Ainsi, on a :

t12 = kxT, et t13 = IxT. Thus, we have:

t12 = kxT, and t13 = IxT.

Les étapes suivantes du procédé, à partir de l'étape 130 sont identiques aux étapes de l'algorithme illustré sur la figure 2. The following steps of the method, starting from step 130, are identical to the steps of the algorithm illustrated in FIG.

Claims

CLAIMS 1. - Method for determining the time shift (t12) between the instants of passage of the same pulse wave at two different measurement points (Mi, M2) of an arterial network of a living being, comprising the steps to: a) placing a sensor of a parameter representative of the pulse wave (14,16) at each measurement point (Mi, M2) of the arterial network; b) record the evolution over time of the parameter representative of the pulse wave at each measurement point (Mi, M2); c) deducing, from the two records ((an), (bn)), said time offset (t12), characterized in that the time-lag deduction step (t12) comprises the steps of: c1) determining a first moment (tp) where a derivative of a predetermined order with respect to the time of the first record (an) is maximum; c2) selecting a reference time interval (IR) containing said first instant (tp); c3) evaluating the time difference between said reference time interval (IR) and a moving time interval (lM) of the same duration for which a distance (d) between the range of values of the first record (an) corresponding to said interval of time reference time (IR) and the range of values of the second record corresponding to said moving time interval (1M) is minimal, the distance being defined by the sum of the absolute values of the point-to-point differences between the ranges of the values of the first and second recordings ((an), (bn)) brought back to the same time origin; and c4) consider the time offset (t12) to be determined as equal to the evaluated time offset.

2. - Method according to claim 1, characterized in that said first instant (tp) is determined as the time when the first derivative with respect to time ((a'n)) of the first record ((an)) is maximum .

3. - Method according to claim 2, characterized in that said reference time interval (IR) is an interval ending substantially at said first instant (tp).

4. The method as claimed in claim 1, characterized in that said first instant (tp) is determined as the time when the second derivative with respect to the time ((a "n)) of the first record ((an)) is maximum .

substantially mid-duration of said interval (IR).

5. - Method according to claim 4, characterized in that said reference time interval (IR) is such that said first instant (tp) is

6. - Method according to any one of the preceding claims, characterized in that said derivative of predetermined order is evaluated by successive application of a finite difference method of order 1.

7. A method according to any one of the preceding claims, characterized in that said representative parameter of the pulse wave is the arterial pressure.

8. A method for determining the propagation velocity (V12) of a pulse wave between two distant measurement points (M1, M2) of an arterial network of a living being, comprising the steps of:

A) measuring the distance (dis) separating the two measurement points (min, M2) on the living being,

B) determining the time shift (tel2) between the instants of passage of the same pulse wave at the two measurement points (Mi, M2);

C) calculating the propagation velocity from the distance (dis) separating the two measurement points and the time shift (tis), characterized in that the temporal shift (tis) between the instants of passage of the pulse wave at the two measuring points (me, M2) is determined by a method according to any one of the preceding claims.

9. A method for estimating central pulsating pressure, characterized in that it comprises the steps of: a) determining the central propagation velocity of a pulse wave between a point (Mi) of the carotid artery and a point (M2) of the femoral artery by carrying out a method according to claim 8,

b) determining the brachial propagation velocity of a pulse wave between a point (Mi) of the carotid artery and a point (M3) of the brachial artery, by carrying out a method according to claim 8, c) measuring the pulse pressure in the upper limb; d) estimate the central pressure from the brachial and central propagation velocities and the pulsed pressure in the upper limb.

where: Va is the central propagation velocity V is the brachial propagation velocity, and Pb is the pulsed pressure in the upper limb.

10. Process according to claim 9, characterized in that the central pulsating pressure Pa is determined by the relation:

11. - Diagnostic installation characterized in that it comprises: - a database of reference values, - means for determining by implementing a method according to any one of the preceding claims of one of the magnitudes included in the group consisting of the temporal shift between the instants of passage of the same pulse wave at two measurement points distant from an arterial network of a living being, the speed of propagation of a pulse wave between two measuring points distant from an arterial network of a living being and the central pulsating pressure of a living being; and means for comparing said determined quantity with at least one reference value.