CN110141196A - Peripheral arterial blood vessel elasticity evaluation method and system based on double triangle flow models - Google Patents
Peripheral arterial blood vessel elasticity evaluation method and system based on double triangle flow models Download PDFInfo
- Publication number
- CN110141196A CN110141196A CN201910509662.4A CN201910509662A CN110141196A CN 110141196 A CN110141196 A CN 110141196A CN 201910509662 A CN201910509662 A CN 201910509662A CN 110141196 A CN110141196 A CN 110141196A
- Authority
- CN
- China
- Prior art keywords
- wave
- pulse wave
- point
- peripheral arterial
- waveform
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/02—Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
- A61B5/02007—Evaluating blood vessel condition, e.g. elasticity, compliance
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7203—Signal processing specially adapted for physiological signals or for diagnostic purposes for noise prevention, reduction or removal
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7235—Details of waveform analysis
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7271—Specific aspects of physiological measurement analysis
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Heart & Thoracic Surgery (AREA)
- Public Health (AREA)
- Physiology (AREA)
- Biophysics (AREA)
- Pathology (AREA)
- Signal Processing (AREA)
- Biomedical Technology (AREA)
- Veterinary Medicine (AREA)
- Medical Informatics (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Physics & Mathematics (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Psychiatry (AREA)
- Vascular Medicine (AREA)
- Cardiology (AREA)
- Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)
Abstract
The present invention provides a kind of peripheral arterial blood vessel elasticity evaluation method and system based on double triangle flow models, is related to medical treatment and healthy electronic technology field.The present invention passes through acquisition peripheral arterial pulse wave first, and to pulse wave removal Hz noise, baseline drift these noises.Characteristic point is extracted from filtered signal, is fitted flow velocity waveforms using double triangle flow models according to characteristic point position, pulse wave is carried out on the basis of double triangle flow velocity waveforms and decomposes to obtain forward-wave and retonation wave.The retonation wave Augmentation index RM and reflection index RI of reflection peripheral arterial blood vessel elasticity are finally calculated according to the feature point value of the forward-wave of extraction and retonation wave.The present invention carries out the accurate decomposition of peripheral arterial pulse wave according to fitting blood flow waveform, relevant parameter is calculated with result after decomposition, the characteristics of comprehensive forward-wave and retonation wave, it realizes and the vascular stiffness degree assessment of peripheral arterial and systema cariovasculare functional is evaluated, so that more convenient to cardiovascular function detection.
Description
Technical field
The present invention relates to medical treatment and healthy electronic technology field more particularly to a kind of periphery based on double triangle flow models
Arteries photoelastic evaluation method and system.
Background technique
Pulse wave is able to reflect periodic vascular pressure variation, can analyze diastole by peripheral arterial pulse waveform
Phase index and contraction phase index, so that it is determined that cardiovascular risk index and blood vessel elasticity indexes.
The patent of invention of Patent No. 201010558798.3 extracts the characteristic point of pulse wave by wavelet transformation, realizes
Vascular stiffness is detected by pulse wave velocity, but characteristic point used is directly extracted from pulse wave, cannot achieve pair
The accurate assessment of vascular stiffness;The patent of invention of Patent No. 201510955386.6 proposes a kind of by left and right finger tip
Blood volume pulsation wave acquisition, realize function of vascular endothelium index calculate, using only arm arteria brachialis block mode surveyed
Amount, is not easy to measure;The patent of invention of Patent No. 201811213852.3 proposes a kind of human body based on support vector machines
Pulse detection method measures pulse wave by portable wrist strap arteries and veins device and extracts feature, according to existing model to pulse wave into
Row classification assessment hardened condition, but the invention is only extracted time and frequency domain characteristics, and the feature for not carrying out pulse Wave Decomposition deeply mentions
It takes;Pulse wave is to be propagated in human body along arterial tree, the patent of invention of Patent No. 201710369140.X, by obtaining
The waveform for taking peripheral arterial to change over time obtains propagation time of the blood flow from central aortic to peripheral arterial, realizes central aortic
The measurement of blood pressure, but this method lacks the step of carrying out decomposition acquisition forward-wave to peripheral arterial, it is difficult to realize the propagation time
It is accurate to calculate.
Summary of the invention
It is a kind of based on double triangle blood flows the technical problem to be solved by the present invention is in view of the above shortcomings of the prior art, provide
The peripheral arterial blood vessel elasticity evaluation method and system of model are realized to pulse wave decomposition, and then evaluate peripheral arterial blood vessel
Elasticity.
To solve the above-mentioned problems, on the one hand, the technical solution adopted by the present invention is that: one kind is based on double triangle flow models
Peripheral arterial blood vessel elasticity evaluation method, comprising the following steps:
Step 1, measurement peripheral arterial at pressure waveform, obtain pulse wave signal, and it is removed Hz noise,
These noises of baseline drift;The peripheral arterial is radial artery, arteria brachialis, arteria carotis or femoral artery;
Step 2 is fitted flow velocity waveforms using double triangle flow models, and is carried out on the basis of the waveform to pulse wave
It decomposes, obtains forward-wave and retonation wave waveform, method particularly includes:
Step 2.1 extracts characteristic point in the pulse wave signal of acquisition;The characteristic point includes pulse wave at peripheral arterial
Starting point, reflection point, incisura point, dicrotic wave peak point and pulse wave end point;
Pulse wave starting point is determined by iteration tangential method at the peripheral arterial, method particularly includes:
(1) first derivative for calculating peripheral arterial rising edge, finds rising edge gradient maxima point;
(2) baseline of pulse wave is extracted;
(3) tangent line on pulse wave about gradient maxima point is done, the intersection point of tangent line Yu pulse wave baseline, intersection point position are found
Set the initial time as pulse wave;
(4) it using the intersection point of pulse wave baseline and its gradient maxima point tangent line as identification point, is projected to pulse wave, gained
Subpoint is as pulse wave initial point;
The reflection point is that diastasis heart starts to penetrate the mark of blood, is determined according to the quality of pulse wave, such as
Dicrotic wave waveform in fruit pulse wave is obvious, then calculates the second dervative of measured pulse wave, take first of second dervative
Otherwise reflection point of the minimum point as pulse wave at peripheral arterial will be used as reflection point at the 25% of ejection time;
The incisura point is the mark terminated the cardiac ejection phase, by second maximum of points of the second dervative of pulse wave
It determines;
The dicrotic wave peak point is the peak point of peripheral arterial end diastolic velocity, according to the quality of pulse wave into
Row determines, if the dicrotic wave waveform in pulse wave is obvious, dicrotic wave peak point is determined using difference threshold algorithm, with dicrotic pulse
Crest value and time of occurrence given threshold, otherwise according to the position of incisura point, finding first peak point backward is dicrotic wave
Peak point position;
The determination method of the pulse wave end point is identical as the tangential method of pulse wave starting point, pulse wave next cycle
Starting point be upper a cycle pulse wave end point;
The maximum value or minimum value of the pulse Propagation guide function, are all obtained using curve-fitting method;The curve is quasi-
Conjunction method is fitted pulse wave segmentation using lower order polynomial expressions the following steps are included: the first step;Second step, it is multinomial to being fitted
Formula is differentiated;Third step is extracted characteristic point from the derived function of polynomial fitting, and is corresponded in pulse wave;Curve matching side
Method is for accurately extracting used second dervative minimum value, maximum value, the maximum value of first derivative;
Step 2.2 constructs peripheral arterial flow velocity waveforms using double triangle flow models;It is extracted according to step 2.1
Five characteristic points construct double triangle blood flow waveforms, and double triangle blood flow waves are by two triangulars at embodying actual flow of blood
The feature of velocity wave form;
For first triangle using first minimum value of the second dervative of pulse wave as vertex position, this is peripheral arterial
Reflection point;The bottom edge of first triangle is determined by starting point to the incisura point position of pulse wave;The incisura point is arteries and veins
Fight wave second order derived function second maximum value where position;Second vertex of a triangle is located at peripheral arterial
At dicrotic pulse crest value, i.e. the last one peak value of peripheral arterial;The bottom edge of second triangle is by incisura point position and arteries and veins
Fight wave end point position determine;
Step 2.3 carries out pulse wave decomposition on the basis of the flow velocity waveforms of double triangle blood flow waveforms building, obtains
Forward-wave and retonation wave waveform, shown in following formula:
Wherein, PmIt (t) is the pulse pressure wave of actual measurement, Q (t) is the Hemodynamic environment using double triangle flow models fittings
It writes music line, PfIt (t) is forward-wave, PbIt (t) is retonation wave, t represents time, ZcFor the characteristic impedance of peripheral arterial, calculation method
Are as follows:
Zc=(Pr-Pd)/H
Wherein, PrIt is pulse wave reflection point pressure value, PdIt is end diastolic pressure force value, H is the amplitude of flow velocity waveforms;
Step 2.4, after being decomposed to pulse wave, a starting point in forward-wave is acquired using tangential method, if preceding
The starting point of Xiang Bozhong is different from the starting point in pulse wave, then the initial point position of pulse wave is replaced with pulse wave starting point
Position is averaged with forward-wave initial point position, and return step 2.2, until the starting extracted from pulse wave and forward-wave
Point is in the same position;
Step 3 is removed baseline drift processing to the obtained forward-wave of decomposition and retonation wave waveform, and to treated
Forward-wave and retonation wave extract characteristic point;The characteristic point includes: preceding to crest value, backward crest value, dicrotic pulse wave height;
Step 4 calculates reflection peripheral arterial blood vessel elasticity according to the forward-wave of extraction and the feature point value of retonation wave
Retonation wave Augmentation index RM and reflection index RI;
The calculation formula of the RM and RI is as follows:
RI=Pb/(Pf+Pb)
Wherein, PfTo be preceding to crest value, PbTo be rear to crest value;
Step 5, according to the retonation wave Augmentation index RM that is calculated and reflection index RI to periphery arteries elasticity into
Row analysis.
On the other hand, the present invention also provides a kind of, and the peripheral arterial blood vessel elasticity based on double triangle flow models evaluates system
System, including pulse wave acquisition unit, transmission unit, waveform decomposition unit and vascular function assessment unit;
The pulse wave acquisition unit obtains pulse wave, and pass through transmission for measuring pressure waveform at peripheral arterial
Collected pulse wave is transferred to waveform decomposition unit and vascular function assessment unit by unit;The waveform decomposition unit will be adopted
The pulse wave that collection obtains is decomposed, and forward-wave and retonation wave waveform is obtained by calculation;The vascular function assessment unit root
Peripheral arterial vascular function is assessed according to the Parameter analysis of forward-wave and retonation wave, and by the assessment result and pulse of vascular function
Wave and decomposition result simultaneous display;
The pulse wave acquisition unit, including velcro cuff and pressure sensor;The velcro cuff is used for pressure
The precise measurement of peripheral arterial pulse wave is realized in the fixation of sensor;The pressure sensor is used for precise measurement peripheral arterial
Locate the pressure change of skin surface, and pressure value is sent by transmission unit to waveform decomposition unit;
The waveform decomposition unit, including signal conditioning module, blood flow velocity fitting module and waveform decomposing module;It is described
Signal conditioning module, for the pulse wave signal of measurement to be filtered and calibrated and is passed to blood flow velocity fitting module;It is described
Blood flow velocity fitting module is fitted periphery artery flow velocity waveform using double triangle flow models;The waveform decomposes
The result that blood flow velocity is fitted is used for waveform and decomposed by module;
The vascular function assessment unit includes blood vessel elasticity evaluation module and oscillography module;The blood vessel elasticity assesses mould
Block is used for the calculating of blood vessel elasticity parameter;The blood vessel elasticity parameter includes reflection index RI, retonation wave Augmentation index RM;It is described
Oscillography module involves backward for displaying in real time the peripheral arterial pulse waveform that measurement obtains and the forward direction that decomposition obtains
Wave, and comprehensive analysis parameter calculated result, show periphery arteries flexibility analysis and its functional assessment.
Preferably, the pressure sensor of the measurement pulse wave, uses the varistor of cantilever beam structure;The cantilever beam
Structure, including probe and varistor two parts;The probe is fixed on peripheral arterial, varistor position by velcro wrist strap
In probe rear.
The beneficial effects of adopting the technical scheme are that the peripheral arterial blood vessel based on double triangle flow models
Photoelastic evaluation method and system, constructed system are mainly the peripheral arterial pulse Wave Decomposition side according to double triangle flow models
Method and build, for peripheral arterial, the characteristics of for its systole phase and end diastolic velocity waveform, use first triangle
Systole phase flow velocity waveforms are fitted, are fitted end diastolic velocity waveform using second triangle, which can be accurate
Reflection actual flow of blood velocity wave form whole features.It is accurately decomposed according to the pulse wave that fitting blood flow waveform carries out radial artery,
The characteristics of being calculated relevant parameter with result after decomposition, integrated forward-wave and retonation wave, realizes the vascular sclerosis journey to peripheral arterial
Degree and its functional evaluation.
Detailed description of the invention
Fig. 1 is the peripheral arterial blood vessel elasticity evaluation method provided in an embodiment of the present invention based on double triangle flow models
Flow chart;
Fig. 2 is the peripheral arterial blood vessel elasticity evaluation system provided in an embodiment of the present invention based on double triangle flow models
Structural block diagram;
Fig. 3 is the radial artery pulse wave waveform that acquisition provided in an embodiment of the present invention obtains;
Fig. 4 be it is provided in an embodiment of the present invention by double triangle flow models construct peripheral arterial flow velocity waveforms with
The comparison of actual flow of blood velocity wave form;
Fig. 5 is pulse wave decomposition result figure provided in an embodiment of the present invention.
Specific embodiment
With reference to the accompanying drawings and examples, specific embodiments of the present invention will be described in further detail.Implement below
Example is not intended to limit the scope of the invention for illustrating the present invention.
In the present embodiment, the peripheral arterial blood vessel elasticity evaluation method based on double triangle flow models, as shown in Figure 1, packet
Include following steps:
Pressure waveform at step 1, measurement peripheral arterial (such as radial artery, arteria brachialis, arteria carotis or femoral artery), obtains arteries and veins
It fights wave signal, and is removed Hz noise, baseline drift these noises to it;
Step 2 is fitted flow velocity waveforms using double triangle flow models, and is carried out on the basis of the waveform to pulse wave
It decomposes, obtains forward-wave and retonation wave waveform, method particularly includes:
Step 2.1 extracts characteristic point in the pulse wave signal of acquisition;The characteristic point includes pulse wave at peripheral arterial
Starting point, reflection point, incisura point, dicrotic wave peak point and pulse wave end point;
The starting point of pulse wave is determined by iteration tangential method at the peripheral arterial, method particularly includes:
(1) first derivative for calculating peripheral arterial rising edge, finds rising edge gradient maxima point;
(2) baseline of pulse wave is extracted;
(3) tangent line on pulse wave about gradient maxima point is done, the intersection point of tangent line Yu pulse wave baseline, intersection point position are found
Set the initial time as pulse wave;
(4) it using the intersection point of pulse wave baseline and its gradient maxima point tangent line as identification point, is projected to pulse wave, gained
Subpoint is as pulse wave initial point;
The reflection point is that diastasis heart starts to penetrate the mark of blood, is determined according to the quality of pulse wave, such as
Dicrotic wave waveform in fruit pulse wave is obvious, then calculates the second dervative of measured pulse wave, take first of second dervative
Otherwise reflection point of the minimum point as pulse wave at peripheral arterial will be used as reflection point at the 25% of ejection time;
The incisura point is the mark terminated the cardiac ejection phase, by second maximum of points of the second dervative of pulse wave
It determines;
The dicrotic wave peak point is the peak point of peripheral arterial end diastolic velocity, according to the quality of pulse wave into
Row determines, if the dicrotic wave waveform in pulse wave is obvious, dicrotic wave peak point is determined using difference threshold algorithm, with dicrotic pulse
Crest value and time of occurrence given threshold, otherwise according to the position of incisura point, finding first peak point backward is dicrotic wave
Peak point position;
The determination method of the pulse wave end point is identical as the tangential method of pulse wave starting point, pulse wave next cycle
Starting point be upper a cycle pulse wave end point;
The maximum value or minimum value of the pulse Propagation guide function, are all obtained using curve-fitting method;The curve is quasi-
Conjunction method is fitted pulse wave segmentation using lower order polynomial expressions the following steps are included: the first step;Second step, it is multinomial to being fitted
Formula is differentiated;Third step is extracted characteristic point from the derived function of polynomial fitting, and is corresponded in pulse wave;Curve matching side
Method is for accurately extracting used second dervative minimum value, maximum value, the maximum value of first derivative;
Step 2.2 constructs peripheral arterial flow velocity waveforms using double triangle flow models;It is extracted according to step 2.1
Five characteristic points construct double triangle blood flow waveforms, and double triangle blood flow waves are by two triangulars at embodying actual flow of blood
The feature of velocity wave form;
For first triangle using first minimum value of the second dervative of pulse wave as vertex position, this is peripheral arterial
Reflection point;The bottom edge of first triangle is determined by starting point to the incisura point position of pulse wave;The incisura point is arteries and veins
Fight wave second order derived function second maximum value where position;Second vertex of a triangle is located at peripheral arterial
At dicrotic pulse crest value, i.e. the last one peak value of peripheral arterial;The bottom edge of second triangle is by incisura point position and arteries and veins
Fight wave end point position determine;
Step 2.3 carries out pulse wave decomposition on the basis of the flow velocity waveforms of double triangle blood flow waveforms building, obtains
Forward-wave and retonation wave waveform, shown in following formula:
Wherein, PmIt (t) is the pulse pressure wave of actual measurement, Q (t) is the Hemodynamic environment using double triangle flow models fittings
It writes music line, PfIt (t) is forward-wave, PbIt (t) is retonation wave, t represents time, ZcFor the characteristic impedance of peripheral arterial, calculation method
Are as follows:
Zc=(Pr-Pd)/H
Wherein, PrIt is pulse wave reflection point pressure value, PdIt is end diastolic pressure force value, H is the amplitude of flow velocity waveforms;
Step 2.4, after being decomposed to pulse wave, a starting point in forward-wave is acquired using tangential method, if preceding
The starting point of Xiang Bozhong is different from the starting point in pulse wave, then the initial point position of pulse wave is replaced with pulse wave starting point
Position is averaged with forward-wave initial point position, and return step 2.2, until the starting extracted from pulse wave and forward-wave
Point is in the same position;
Step 3 is removed baseline drift processing to the obtained forward-wave of decomposition and retonation wave, and to treated forward direction
Wave and retonation wave extract characteristic point;The characteristic point includes: preceding to crest value, backward crest value and dicrotic pulse wave height;
Step 4 calculates reflection peripheral arterial blood vessel elasticity according to the forward-wave of extraction and the feature point value of retonation wave
Retonation wave Augmentation index RM and reflection index RI;
The calculation formula of the RM and RI is as follows:
RI=Pb/(Pf+Pb)
Wherein, PfTo be preceding to crest value, PbTo be rear to crest value;
Step 5, according to the retonation wave Augmentation index RM that is calculated and reflection index RI to periphery arteries elasticity into
Row analysis.
Peripheral arterial blood vessel elasticity evaluation system based on double triangle flow models, as shown in Fig. 2, including that pulse wave acquires
Unit, transmission unit, waveform decomposition unit and vascular function assessment unit;
The pulse wave acquisition unit obtains pulse wave, and single by transmission for measuring peripheral arterial pressure wave shape
Collected pulse wave is transferred to waveform decomposition unit and vascular function assessment unit by member;The waveform decomposition unit will acquire
The pulse wave of acquisition is decomposed, and forward-wave and retonation wave is obtained by calculation;The vascular function assessment unit according to before to
The Parameter analysis of wave and retonation wave assesses cardiovascular function, and the assessment result of vascular function and pulse wave and decomposition result is same
Step display;
The pulse wave acquisition unit, including velcro cuff and pressure sensor;The velcro cuff is used for pressure
The precise measurement of peripheral arterial pulse wave is realized in the fixation of sensor;The pressure sensor is used for precise measurement peripheral arterial
Locate the pressure change of skin surface, and pressure value is sent by transmission unit to waveform decomposition unit;Pressure sensor uses outstanding
The varistor of arm girder construction, being capable of the contacted skin surface pressure change of accurate detection probe;The cantilever beam structure,
Including probe and varistor two parts;The probe is fixed on peripheral arterial by velcro wrist strap, and varistor, which is located at, to be visited
Head rear, due to unique beam type structure, can it is sensitive, accurately return to pulse waveform.
The waveform decomposition unit, including signal conditioning module, blood flow velocity fitting module and waveform decomposing module;It is described
Signal conditioning module, for the pulse wave signal of measurement to be filtered and calibrated and is passed to blood flow velocity fitting module;It is described
Blood flow velocity fitting module is using double fittings of the triangle flow model to flow velocity waveforms at peripheral arterial, by measurement
Characteristic point is extracted on pressure waveform, and flow velocity waveforms can be made with more accurate approximation;The waveform decomposing module, will
The result of blood flow velocity fitting is decomposed for waveform, and the method has been proven that its feasibility through Westerhof et al., for outer
All arteries, the present invention decompose pulse wave to trigonometry and are improved, and the approximation of flow velocity waveforms is done using two triangles,
It realizes with accurately decomposing;
The vascular function assessment unit includes blood vessel elasticity evaluation module and oscillography module;The blood vessel elasticity assesses mould
Block is used for the calculating of blood vessel elasticity parameter;The blood vessel elasticity parameter includes reflection index RI, radial artery retonation wave Augmentation index
RM;The oscillography module involves backward for displaying in real time the peripheral arterial waveform that measurement obtains and the forward direction that decomposition obtains
Wave, and comprehensive analysis parameter calculated result, show periphery arteries flexibility analysis and its functional assessment.
In the present embodiment, the sample frequency of pressure sensor is 1000 hertz, acquires the pulse wave at radial artery, waveform is such as
Shown in Fig. 3, collected signal is pre-processed, the noises such as removal Hz noise, baseline drift, pretreated clean letter
It number is sent by transmission module to waveform decomposing module;The flow velocity waveforms and actual flow of blood constructed by double triangle flow models
The comparison of velocity wave form is as shown in figure 4, pulse wave decomposition result is as shown in Figure 5;As can be seen from Figure 4 and Figure 5, the method for the present invention
Whole features of actual flow of blood velocity wave form can accurately be reflected by passing through double triangle flow models with system, according to fitting blood flow
Waveform realizes the accurate decomposition to radial artery pulse wave.
Finally, it should be noted that the above embodiments are merely illustrative of the technical solutions of the present invention, rather than its limitations;Although
Present invention has been described in detail with reference to the aforementioned embodiments, those skilled in the art should understand that: it still may be used
To modify to technical solution documented by previous embodiment, or some or all of the technical features are equal
Replacement;And these are modified or replaceed, model defined by the claims in the present invention that it does not separate the essence of the corresponding technical solution
It encloses.
Claims (6)
1. a kind of peripheral arterial blood vessel elasticity evaluation method based on double triangle flow models, it is characterised in that: including following step
It is rapid:
Pressure waveform at step 1, measurement peripheral arterial, obtains pulse wave signal, and Hz noise, baseline are removed to it
Drift about these noises;The peripheral arterial is radial artery, arteria brachialis, arteria carotis or femoral artery;
Step 2 is fitted flow velocity waveforms using double triangle flow models, and is divided on the basis of the waveform pulse wave
Solution obtains forward-wave and retonation wave waveform, method particularly includes:
Step 2.1 extracts characteristic point in the pulse wave signal of acquisition;The characteristic point includes of pulse wave at peripheral arterial
Initial point, reflection point, incisura point, dicrotic wave peak point and pulse wave end point;
Step 2.2 constructs peripheral arterial flow velocity waveforms using double triangle flow models;Five extracted according to step 2.1
Characteristic point constructs double triangle blood flow waveforms, and the triangular wave is by two triangulars at embodying actual flow of blood velocity wave form
Feature;
For first triangle using first minimum value of the second dervative of pulse wave as vertex position, this is peripheral arterial reflection
Point;The bottom edge of first triangle is determined by starting point to the incisura point position of pulse wave;The incisura point is pulse wave
Second order derived function second maximum value where position;Second vertex of a triangle is located at the dicrotic pulse of peripheral arterial
At crest value, i.e. the last one peak value of peripheral arterial;The bottom edge of second triangle is by incisura point position and pulse wave
End point position determines;
Step 2.3, by double triangle blood flow waveforms building flow velocity waveforms on the basis of carry out pulse wave decomposition, obtain before to
Wave and retonation wave waveform;
Step 2.4, after being decomposed to pulse wave, a starting point in forward-wave is acquired using tangential method, if forward-wave
In starting point it is different from the starting point in pulse wave, then the initial point position of pulse wave is replaced with into pulse wave initial point position
With being averaged for forward-wave initial point position, and return step 2.2, until the starting point extracted from pulse wave and forward-wave exists
The same position;
Step 3 is removed baseline drift processing to the obtained forward-wave of decomposition and retonation wave, and to treated forward-wave and
Retonation wave extracts characteristic point;The characteristic point includes: preceding to crest value, backward crest value, dicrotic pulse wave height;
Step 4, the retonation wave that anti-peripheral arterial blood vessel elasticity is calculated according to the forward-wave of extraction and the feature point value of retonation wave
Augmentation index RM and reflection index RI;
Step 5 divides periphery arteries elasticity according to the retonation wave Augmentation index RM and reflection index RI that are calculated
Analysis.
2. the peripheral arterial blood vessel elasticity evaluation method according to claim 1 based on double triangle flow models, feature
Be: the starting point of pulse wave is determined by iteration tangential method at peripheral arterial described in step 2.1, method particularly includes:
(1) first derivative for calculating peripheral arterial rising edge, finds rising edge gradient maxima point;
(2) baseline of pulse wave is extracted;
(3) tangent line on pulse wave about gradient maxima point is done, the intersection point of tangent line Yu pulse wave baseline is found, intersection position is
For the initial time of pulse wave;
(4) it using the intersection point of pulse wave baseline and its gradient maxima point tangent line as identification point, is projected to pulse wave, gained projection
Point is used as pulse wave initial point;
The reflection point is that diastasis heart starts to penetrate the mark of blood, is determined according to the quality of pulse wave, if arteries and veins
The dicrotic wave waveform fought in waveform is obvious, then calculates the second dervative of measured pulse wave, take first minimum of second dervative
Reflection point of the value point as pulse wave at peripheral arterial, otherwise, by the 25% of ejection time as reflection point;
The incisura point is the mark terminated the cardiac ejection phase, and second maximum of points by the second dervative of pulse wave is true
It is fixed;
The dicrotic wave peak point is the peak point of peripheral arterial end diastolic velocity, is carried out according to the quality of pulse wave true
It is fixed, if the dicrotic wave waveform in pulse wave is obvious, dicrotic wave peak point is determined using difference threshold algorithm, with dicrotic pulse wave crest
Value and time of occurrence given threshold, otherwise according to the position of incisura point, finding first peak point backward is dicrotic pulse crest value
Point position;
The determination method of the pulse wave end point is identical as the tangential method of pulse wave starting point, and pulse wave next cycle rises
Initial point is the end point of upper a cycle pulse wave;
The maximum value or minimum value of the pulse Propagation guide function, are all obtained using curve-fitting method;The curve matching side
Method is fitted pulse wave segmentation using lower order polynomial expressions the following steps are included: the first step;Second step seeks polynomial fitting
Derivative;Third step is extracted characteristic point from the derived function of polynomial fitting, and is corresponded in pulse wave;Curve-fitting method is used
Second dervative minimum value, maximum value, the maximum value of first derivative used in accurate extract.
3. the peripheral arterial blood vessel elasticity evaluation method according to claim 1 based on double triangle flow models, feature
It is: shown in the forward-wave and the following formula of retonation wave that step 2.3 obtains:
Wherein, PmIt (t) is the pulse pressure wave of actual measurement, Q (t) is to be write music using the Hemodynamic environment of double triangle flow models fittings
Line, PfIt (t) is forward-wave, PbIt (t) is retonation wave, t represents time, ZcFor the characteristic impedance of peripheral arterial, calculation method are as follows:
Zc=(Pr-Pd)/H
Wherein, PrIt is pulse wave reflection point pressure value, PdIt is end diastolic pressure force value, H is the amplitude of flow velocity waveforms.
4. the peripheral arterial blood vessel elasticity evaluation method according to claim 3 based on double triangle flow models, feature
Be: the calculation formula of retonation wave Augmentation index RM and reflection index RI described in step 4 are as follows:
RI=Pb/(Pf+Pb)
Wherein, PfTo be preceding to crest value, PbTo be rear to crest value.
5. a kind of peripheral arterial blood vessel elasticity evaluation system based on double triangle flow models is based on base described in claim 1
It is evaluated in the peripheral arterial blood vessel elasticity evaluation method of double triangle flow models, it is characterised in that: acquired including pulse wave
Unit, transmission unit, waveform decomposition unit and cardiovascular function assessment unit;
The pulse wave acquisition unit obtains pulse wave, and pass through transmission unit for measuring pressure waveform at peripheral arterial
Collected pulse wave is transferred to waveform decomposition unit and cardiovascular function assessment unit;The waveform decomposition unit will acquire
The pulse wave of acquisition is decomposed, and forward-wave and retonation wave is obtained by calculation;The cardiovascular function assessment unit is before
Cardiovascular function is assessed to the Parameter analysis of wave and retonation wave, and by the assessment result of vascular function and pulse wave and decomposes knot
Fruit simultaneous display;
The pulse wave acquisition unit, including velcro cuff and pressure sensor;The velcro cuff is used for pressure sensing
The precise measurement of peripheral arterial pulse wave is realized in the fixation of device;The pressure sensor, for skin at precise measurement peripheral arterial
The pressure change on skin surface, and pressure value is sent by transmission unit to waveform decomposition unit;
The waveform decomposition unit, including signal conditioning module, blood flow velocity fitting module and waveform decomposing module;The signal
Conditioning module, for the pulse wave signal of measurement to be filtered and calibrated and is passed to blood flow velocity fitting module;The blood flow
Velocity fitting module is fitted flow velocity waveforms at peripheral arterial using double triangle flow models;The waveform decomposes mould
The result that blood flow velocity is fitted is used for waveform and decomposed by block;
The cardiovascular function assessment unit includes blood vessel elasticity evaluation module and oscillography module;The painstaking effort tubular elastic assesses mould
Block is used for the calculating of blood vessel elasticity parameter;The blood vessel elasticity parameter includes reflection index RI, radial artery retonation wave Augmentation index
RM;The oscillography module involves backward for displaying in real time the peripheral arterial waveform that measurement obtains and the forward direction that decomposition obtains
Wave, and comprehensive analysis parameter calculated result show periphery arteries flexibility analysis and cardiovascular function assessment.
6. the peripheral arterial blood vessel elasticity evaluation system according to claim 5 based on double triangle flow models, feature
Be: the pressure sensor of the measurement pulse wave uses the varistor of cantilever beam structure;The cantilever beam structure, including
Probe and varistor two parts;The probe is fixed on peripheral arterial by velcro wrist strap, after varistor is located at probe
Side.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910509662.4A CN110141196B (en) | 2019-06-13 | 2019-06-13 | Peripheral arterial vessel elasticity evaluation system based on double-triangle blood flow model |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910509662.4A CN110141196B (en) | 2019-06-13 | 2019-06-13 | Peripheral arterial vessel elasticity evaluation system based on double-triangle blood flow model |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110141196A true CN110141196A (en) | 2019-08-20 |
CN110141196B CN110141196B (en) | 2021-06-01 |
Family
ID=67591217
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910509662.4A Active CN110141196B (en) | 2019-06-13 | 2019-06-13 | Peripheral arterial vessel elasticity evaluation system based on double-triangle blood flow model |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110141196B (en) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111508077A (en) * | 2020-04-29 | 2020-08-07 | 中国人民解放军总医院 | Method, device, equipment and medium for intelligently evaluating heart function by using double-triangle model |
CN112274120A (en) * | 2020-10-28 | 2021-01-29 | 河北工业大学 | Noninvasive arteriosclerosis detection method and device based on one-way pulse wave |
CN112826468A (en) * | 2019-11-05 | 2021-05-25 | 深圳市大富智慧健康科技有限公司 | Blood pressure detection device, blood pressure detection system and blood pressure monitoring method |
CN112932423A (en) * | 2021-01-25 | 2021-06-11 | 中山大学附属第八医院(深圳福田) | Cardiovascular and cerebrovascular disease prediction method, system and equipment based on external counterpulsation intervention |
CN113010969A (en) * | 2021-05-11 | 2021-06-22 | 中国空气动力研究与发展中心计算空气动力研究所 | Calculation and correction method for maximum value of pressure wave of high-speed train entering tunnel |
CN113951846A (en) * | 2021-12-17 | 2022-01-21 | 北京麦邦光电仪器有限公司 | Pulse wave signal processing method and device and readable storage medium |
CN114190902A (en) * | 2021-11-22 | 2022-03-18 | 湖南云医链生物科技有限公司 | Body health index characterization method, device, equipment and storage medium |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6527725B1 (en) * | 2001-01-25 | 2003-03-04 | Colin Corporation | Blood pressure estimating apparatus |
WO2007062456A1 (en) * | 2005-12-01 | 2007-06-07 | Atcor Medical Pty Ltd | A method of estimating pulse wave velocity |
AU2006319737A1 (en) * | 2005-12-01 | 2007-06-07 | Atcor Medical Pty Ltd | A method of estimating pulse wave velocity |
CN102894982A (en) * | 2012-09-28 | 2013-01-30 | 北京工业大学 | Non-invasive detecting method for blood viscosity based on pulse wave |
CN103690152A (en) * | 2014-01-06 | 2014-04-02 | 山东大学 | Arterial elasticity evaluating device based on pulse analysis |
CN106691406A (en) * | 2017-01-05 | 2017-05-24 | 大连理工大学 | Detection method of vascular elasticity and blood pressure based on single probe photoplethysmography pulse wave |
US20180263506A1 (en) * | 2017-03-13 | 2018-09-20 | Rutgers, The State University Of New Jersey | Methods and apparatus for express estimation of the arterial elastic property in a subject |
-
2019
- 2019-06-13 CN CN201910509662.4A patent/CN110141196B/en active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6527725B1 (en) * | 2001-01-25 | 2003-03-04 | Colin Corporation | Blood pressure estimating apparatus |
WO2007062456A1 (en) * | 2005-12-01 | 2007-06-07 | Atcor Medical Pty Ltd | A method of estimating pulse wave velocity |
AU2006319737A1 (en) * | 2005-12-01 | 2007-06-07 | Atcor Medical Pty Ltd | A method of estimating pulse wave velocity |
CN102894982A (en) * | 2012-09-28 | 2013-01-30 | 北京工业大学 | Non-invasive detecting method for blood viscosity based on pulse wave |
CN103690152A (en) * | 2014-01-06 | 2014-04-02 | 山东大学 | Arterial elasticity evaluating device based on pulse analysis |
CN106691406A (en) * | 2017-01-05 | 2017-05-24 | 大连理工大学 | Detection method of vascular elasticity and blood pressure based on single probe photoplethysmography pulse wave |
US20180263506A1 (en) * | 2017-03-13 | 2018-09-20 | Rutgers, The State University Of New Jersey | Methods and apparatus for express estimation of the arterial elastic property in a subject |
Non-Patent Citations (2)
Title |
---|
AHMED QASEM AND ALBERTO AVOLIO: "Determination of Aortic Pulse Wave Velocity From Waveform Decomposition of the Central Aortic Pressure Pulse", 《HYPERTENSION》 * |
叶骏: "基于单个颈动脉压力波形获取主动脉脉搏波速度的可行性研究", 《北京生物医学工程》 * |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112826468A (en) * | 2019-11-05 | 2021-05-25 | 深圳市大富智慧健康科技有限公司 | Blood pressure detection device, blood pressure detection system and blood pressure monitoring method |
CN112826468B (en) * | 2019-11-05 | 2023-12-01 | 深圳市大富智慧健康科技有限公司 | Blood pressure detection device, blood pressure detection system, and blood pressure monitoring method |
CN111508077A (en) * | 2020-04-29 | 2020-08-07 | 中国人民解放军总医院 | Method, device, equipment and medium for intelligently evaluating heart function by using double-triangle model |
CN112274120A (en) * | 2020-10-28 | 2021-01-29 | 河北工业大学 | Noninvasive arteriosclerosis detection method and device based on one-way pulse wave |
CN112274120B (en) * | 2020-10-28 | 2023-07-25 | 河北工业大学 | Noninvasive arteriosclerosis detection method and device based on one-way pulse wave |
CN112932423A (en) * | 2021-01-25 | 2021-06-11 | 中山大学附属第八医院(深圳福田) | Cardiovascular and cerebrovascular disease prediction method, system and equipment based on external counterpulsation intervention |
CN112932423B (en) * | 2021-01-25 | 2024-05-28 | 中山大学附属第八医院(深圳福田) | Method, system and equipment for predicting cardiovascular and cerebrovascular diseases based on external counterpulsation intervention |
CN113010969A (en) * | 2021-05-11 | 2021-06-22 | 中国空气动力研究与发展中心计算空气动力研究所 | Calculation and correction method for maximum value of pressure wave of high-speed train entering tunnel |
CN114190902A (en) * | 2021-11-22 | 2022-03-18 | 湖南云医链生物科技有限公司 | Body health index characterization method, device, equipment and storage medium |
CN113951846A (en) * | 2021-12-17 | 2022-01-21 | 北京麦邦光电仪器有限公司 | Pulse wave signal processing method and device and readable storage medium |
Also Published As
Publication number | Publication date |
---|---|
CN110141196B (en) | 2021-06-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110141196A (en) | Peripheral arterial blood vessel elasticity evaluation method and system based on double triangle flow models | |
CN106037694B (en) | A kind of continuous blood pressure measurer based on pulse wave | |
EP3157416B1 (en) | System for cuff-less blood pressure (bp) measurement of a subject | |
Kachuee et al. | Cuffless blood pressure estimation algorithms for continuous health-care monitoring | |
Kachuee et al. | Cuff-less high-accuracy calibration-free blood pressure estimation using pulse transit time | |
Forouzanfar et al. | Oscillometric blood pressure estimation: past, present, and future | |
McLaughlin et al. | Piezoelectric sensor determination of arterial pulse wave velocity | |
CN103385702B (en) | A kind of non-invasive blood pressure continuous detection apparatus and method | |
Trachet et al. | Numerical validation of a new method to assess aortic pulse wave velocity from a single recording of a brachial artery waveform with an occluding cuff | |
US6413223B1 (en) | Cuffless continuous blood pressure monitor | |
Sola et al. | Chest pulse-wave velocity: A novel approach to assess arterial stiffness | |
CN103479343B (en) | Central aortic pressure detection system and method based on oscillating sphygmomanometer signals | |
CN101156771A (en) | Method and apparatus for improving vascellum hardness measurement precision base on pulse wave frequency spectrum analysis | |
CN111493855B (en) | System and method for non-invasive measurement of individualized cardiac output | |
CN104042200B (en) | A kind of non-invasive monitoring devices and methods therefor of arteriotony by shooting | |
Bertaglia et al. | Computational hemodynamics in arteries with the one-dimensional augmented fluid-structure interaction system: viscoelastic parameters estimation and comparison with in-vivo data | |
EP2499967A1 (en) | Pulse wave propagation speed measurement device and pulse wave propagation speed measurement program | |
CN104068841B (en) | A kind of measuring method and device measuring Indices of Systolic Time parameter | |
CN113907727B (en) | Beat-by-beat blood pressure measurement system and method based on photoplethysmography | |
Xu et al. | Quantitative analyses of pulse images in Traditional Chinese Medicine | |
Aguado-Sierra et al. | Arterial pulse wave velocity in coronary arteries | |
CN116392091A (en) | Continuous arterial blood pressure waveform detection system and device based on single-channel finger tip PPG | |
Nabeel et al. | Deep learning for blood pressure estimation: an approach using local measure of arterial dual diameter waveforms | |
CN101332084B (en) | Area barycenter displacement analytical method of pulse wave | |
JP2023510943A (en) | System and method for pulse transit time measurement with optical data |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |