CN103284703A - Aortic pulse wave transfer time measuring method based on upper extremity artery information - Google Patents
Aortic pulse wave transfer time measuring method based on upper extremity artery information Download PDFInfo
- Publication number
- CN103284703A CN103284703A CN2013102415279A CN201310241527A CN103284703A CN 103284703 A CN103284703 A CN 103284703A CN 2013102415279 A CN2013102415279 A CN 2013102415279A CN 201310241527 A CN201310241527 A CN 201310241527A CN 103284703 A CN103284703 A CN 103284703A
- Authority
- CN
- China
- Prior art keywords
- aortic
- pulse wave
- artery
- upper extremity
- 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
Images
Landscapes
- Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)
Abstract
The invention provides an aortic pulse wave transfer time measuring method based on upper extremity artery information. The method is characterized by comprising the following steps: (1) collecting the upper extremity artery information composed of a radial artery pulse wave, a finger tip volume pulse wave and a brachial artery blood pressure value; (2) by means of a generalized transfer function model, converting the radial artery pulse wave into an aortic pulse wave which is calibrated by the brachial artery blood pressure value, and obtaining an aortic pulse pressure waveform p(t); (3) by means of a blood flow transfer function model, converting a finger artery blood wave into an aortic blood flow waveform q(t); (4) moving the staring point of the p(t) and the starting point of the q(t) behind the same point, by means of an impedance formula, dividing the p(t) into a forward wave Pf(t) and a reflective wave Pb(t); calculating aortic pulse wave transfer time according to normalization correlation coefficients of the Pf(t) and the Pb(t). By means of collecting of the upper extremity artery information, the aortic pulse wave transfer time can be obtained, testing procedures are simple, test difficulty is low, the upper extremity artery information of a human body is stable, signal quality is good, and testing precision is high.
Description
Technical field
The present invention relates to the cardiovascular medicine field in the biomedicine, be specifically related to a kind of aortic pulse wave passing time measuring method based on artery of upper extremity information.
Background technology
Cardiovascular disease is the murderous main reason of global range." Chinese cardiovascular diseases reports 2010 " shows: China's cardiovascular patient has reached 2.3 hundred million at present, and among per 10 adults 2 people being arranged is cardiovascular diseasess; China has 3,000,000 people to die from the cardiovascular diseases every year approximately, and it is the cardiovascular diseases that every dead 3 philtrums just have 1 people.In JIUYUE, 2011 World Health Organization (WHO) report claims: the whole world had 1,730 ten thousand people to die from the cardiovascular diseases approximately in 2008, accounted for 30% of total death; As not controlled, the year two thousand thirty will have 2,360 ten thousand people to die from the cardiovascular diseases.The cardiovascular diseases has become one of great public health problem of 21st century threat human health.
The essence of cardiovascular disease is the artery structure that causes of arteriosclerosis and the change of spring function.The change of arterial elasticity function is early than the change of structure, and the change of structure is irreversible often.Therefore, the key of reduction cardiovascular disease incidence rate and mortality rate is arteriosclerotic inspection morning, discovery early, intervention early.
Aortic pulse wave speed aoPWV measurement is important noinvasive arterial stiffness assay method, is widely used in therapeutic evaluation and the medication guide of arteriosclerosis earlier detection, intervention treatment.Wherein, aortic pulse wave passing time aoPTT is the key that aoPWV measures.The method of aoPTT detection commonly used comprises two-point method and single-point method at present.The key of measuring aoPTT based on " two-point method " is to identify the starting point of two-way (neck-femoral artery) pulse wave in the same cardiac cycle, calculates their time difference for assessment of aoPTT.But this measuring method may need the patient to expose inguinal region, operates relative complex, and is difficult to obtain stable signal for cervical region or the more patient of stomach fat.The key of measuring aoPTT based on " single-point method " is that accurately identification characterizes the characteristic point of forward-wave and echo starting point or peak value, and the time difference of calculating them is used for estimating aoPTT.But, when experimenter's second contraction crest is not obvious, be difficult to accurate recognition feature point from time domain, influence the certainty of measurement of aoPTT.Therefore, simple to operate, the aoPTT measuring method is demanded urgently breaking through accurately.
Summary of the invention
The present invention is for avoiding above-mentioned existing in prior technology weak point, a kind of aortic pulse wave passing time measuring method based on artery of upper extremity information being provided, in the hope of effectively improving measuring accuracy, simplifying testing process.
The present invention adopts following technical scheme for the technical solution problem:
The characteristics that the present invention is based on the aortic pulse wave passing time measuring method of artery of upper extremity information are: described artery of upper extremity information is made up of radial artery pulse wave, finger volume pulse ripple and brachial arterial pressure value;
Described aortic pulse wave passing time measuring method based on artery of upper extremity information is to carry out as follows:
A, gather radial artery pulse wave at the human body wrist, gather the finger volume pulse ripple at finger tip, and measure the brachial arterial pressure value;
B, utilizing the generalized transfer function model, is input with the Laplace transformation of radial artery pulse wave, obtains the pull-type conversion of aortic pulse wave, and then obtains aortic pulse wave; With the brachial arterial pressure value aortic pulse wave is demarcated, the aortic pulse wave that acquisition has the blood pressure characteristic, described aortic pulse wave with blood pressure characteristic is aortic pressure waveform p (t);
C, described finger volume pulse ripple is carried out independent component analysis, obtain to refer to the artery blood flow waveform; Utilizing the blood flow transfer function model, is input with the Laplace transformation that refers to the artery blood flow wave mode, obtains the pull-type conversion of aortic flow waveform, and then obtains aortic flow waveform q (t);
D, described aortic pressure waveform p (t) and described aortic flow waveform q (t) are adopted Fourier analysis calculated characteristics impedance Z
c
E, the starting point of described aortic pressure waveform p (t) and described aortic flow waveform q (t) moved to same point after, utilize formula of impedance that aortic pressure waveform p (t) is decomposed into forward-wave Pf (t) and echo Pb (t); Described formula of impedance is suc as formula shown in (1) and the formula (2):
Pf(t)=（p(t)+Z
c×q(t)）÷2 （1）
Pb(t)=（p(t)-Z
c×q(t)）÷2 （2）；
The normalized crosscorrelation coefficient of f, calculating forward-wave Pf (t) and echo Pb (t), the maximum time corresponding point of described normalized crosscorrelation coefficient is TR2, then aortic pulse wave passing time aoPTT is TR2/2.
The characteristics that the present invention is based on the aortic pulse wave passing time measuring method of artery of upper extremity information also are: described generalized transfer function model is: the Laplace transformation RP (S) of pull-type conversion AP (the S)/radial artery pulse wave of generalized transfer function GTF (S)=aortic pulse wave.
Described blood flow transfer function model is: the pull-type conversion AB (S) of blood flow transfer function BTF (S)=aortic flow waveform/the refer to pull-type conversion FB (S) of artery blood flow waveform.
Compared with the prior art, beneficial effect of the present invention is embodied in:
(1) the present invention only need gather the passing time that artery of upper extremity information just can be obtained aortic pulse wave, and testing process is simple, difficulty of test is low, and human upper limb tremulous pulse information stability, and signal quality is good, the measuring accuracy height;
(2) compare and traditional " two-point method ", the present invention has transferred to artery of upper extremity (radial artery, refer to tremulous pulse) with collection position from large artery trunks, has overcome the shortcoming that is difficult to obtain stable signal for cervical region or the more patient of stomach fat effectively;
(3) the present invention is incorporated into the fields of measurement of aoPTT with the impedance analysis technology, also has good robustness even put at wave character when not obvious, has well solved the not enough problem of certainty of measurement that " single-point method " exists, and has wider application prospect.
Description of drawings
Fig. 1 obtains the general flow chart of aoPTT method for the present invention;
Fig. 2 obtains the sketch map of aortic pressure Reeb for the present invention by the radial pulse waveform;
Fig. 3 obtains the sketch map of aortic flow waveform by the finger volume pulse ripple for the present invention;
Fig. 4 obtains the sketch map of forward-wave and echo for the present invention utilizes impedance analysis;
Fig. 5 obtains the sketch map of aoPTT for the present invention utilizes the normalized crosscorrelation coefficient of forward-wave and echo.
The specific embodiment
Embodiment 1
Present embodiment based on the characteristics of the aortic pulse wave passing time measuring method of artery of upper extremity information is: artery of upper extremity information is made up of radial artery pulse wave, finger volume pulse ripple, brachial arterial pressure value;
As shown in Figure 1, present embodiment is to carry out as follows based on the aortic pulse wave passing time measuring method of artery of upper extremity information:
1, gathers radial artery pulse wave at the human body wrist, gather the finger volume pulse ripple at finger tip, and measure the brachial arterial pressure value;
2, as shown in Figure 2, utilizing the generalized transfer function model, is input with the Laplace transformation of radial artery pulse wave, obtains the pull-type conversion of aortic pulse wave, and then obtains aortic pulse wave; With the brachial arterial pressure value aortic pulse wave is demarcated, the aortic pulse wave that acquisition has the blood pressure characteristic, described aortic pulse wave with blood pressure characteristic is aortic pressure waveform p (t);
To the periphery transmission, owing to be subjected to the effect of factors such as blood vessel wall and Peripheral resistance, features such as the amplitude of pulse wave, form can change the pulse wave that heart contraction produces along arterial tree.Multinomial clinical research shows, has specific transitive relation between radial artery pulse wave and the aortic pulse wave.The generalized transfer function model that utilizes mathematical method to set up can be used for obtaining aortic pulse wave based on radial artery pulse wave.
The generalized transfer function model is: the Laplace transformation RP (S) of pull-type conversion AP (the S)/radial artery pulse wave of generalized transfer function GTF (S)=aortic pulse wave.
The generalized transfer function model can adopt prior art, also can set up by following method: the radial artery pulse wave and the aortic pulse wave that at first obtain the experimenter, these two kinds of waveforms are obtained the Laplace transformation RP (S) of radial artery pulse wave and the pull-type conversion AP (S) of aortic pulse wave after the pull-type conversion of process respectively, and then obtain generalized transfer function GTF (S).This is single experimenter's generalized transfer function computational methods, and a plurality of experimenters' generalized transfer function is averaged after for the calculating of carrying out single experimenter respectively.The generalized transfer function model has general applicability.
3, as shown in Figure 3, described finger volume pulse ripple is carried out independent component analysis, obtain to refer to the artery blood flow waveform; Utilizing the blood flow transfer function model, is input with the Laplace transformation that refers to the artery blood flow ripple, obtains the pull-type conversion of aortic flow waveform, and then obtains aortic flow waveform q (t);
The finger volume pulse ripple comprises abundant hemodynamics information, by the independent component analysis technology, can extract from the finger volume pulse ripple and refers to the artery blood flow waveform.To the peripheral arterial transmittance process, owing to be subjected to the effect of factors such as blood vessel viscoelasticity, features such as its amplitude, form can change blood flow equally from aorta.The same with the pulse wave transfer law, refer to that also there are a kind of transitive relation in artery blood flow waveform and aortic flow waveform, set up the blood flow transfer function model by this relation, can be used for obtaining the aortic flow waveform.
The blood flow transfer function model is: the pull-type conversion AB (S) of blood flow transfer function BTF (S)=aortic flow waveform/the refer to pull-type conversion FB (S) of artery blood flow waveform.
The blood flow transfer function model can adopt prior art, also can set up by following method: at first obtain experimenter's aortic flow waveform and refer to the artery blood flow waveform, these two kinds of waveforms are obtained the pull-type conversion (AB (S)) of aortic flow waveform and the pull-type conversion (FB (S)) of finger artery blood flow waveform after the pull-type conversion of process respectively, and then obtain blood flow transfer function BTF (S).This is the computational methods of single experimenter's blood flow transfer function model, and a plurality of experimenters' blood flow transfer function model is averaged after for the calculating of carrying out single experimenter respectively.The blood flow transfer function model has general applicability.
4, described aortic pressure waveform p (t) and described aortic flow waveform q (t) are adopted Fourier analysis calculated characteristics impedance Z
c
5, as shown in Figure 4, the starting point of described aortic pressure waveform p (t) and described aortic flow waveform q (t) moved to same point after, utilize formula of impedance that aortic pressure waveform p (t) is decomposed into forward-wave Pf (t) and echo Pb (t); Described formula of impedance is suc as formula shown in (1) and the formula (2):
Pf(t)=（p(t)+Z
c×q(t)）÷2 （1）
Pb(t)=（p(t)-Z
c×q(t)）÷2 （2）；
6, as shown in Figure 5, calculate the normalized crosscorrelation coefficient of forward-wave Pf (t) and echo Pb (t), the maximum time corresponding point of described normalized crosscorrelation coefficient is TR2, and then aortic pulse wave passing time aoPTT is TR2/2.
Claims (3)
1. based on the aortic pulse wave passing time measuring method of artery of upper extremity information, it is characterized in that: described artery of upper extremity information is made up of radial artery pulse wave, finger volume pulse ripple and brachial arterial pressure value;
Described aortic pulse wave passing time measuring method based on artery of upper extremity information is to carry out as follows:
A, gather radial artery pulse wave at the human body wrist, gather the finger volume pulse ripple at finger tip, and measure the brachial arterial pressure value;
B, utilizing the generalized transfer function model, is input with the Laplace transformation of described radial artery pulse wave, obtains the pull-type conversion of aortic pulse wave, and then obtains aortic pulse wave; With described brachial arterial pressure value aortic pulse wave is demarcated, the aortic pulse wave that acquisition has the blood pressure characteristic, described aortic pulse wave with blood pressure characteristic is aortic pressure waveform p (t);
C, described finger volume pulse ripple is carried out independent component analysis, obtain to refer to the artery blood flow waveform; Utilizing the blood flow transfer function model, is input with the Laplace transformation that refers to the artery blood flow wave mode, obtains the pull-type conversion of aortic flow waveform, and then obtains aortic flow waveform q (t);
D, described aortic pressure waveform p (t) and described aortic flow waveform q (t) are adopted Fourier analysis calculated characteristics impedance Z
c
E, the starting point of described aortic pressure waveform p (t) and described aortic flow waveform q (t) moved to same point after, utilize formula of impedance that aortic pressure waveform p (t) is decomposed into forward-wave Pf (t) and echo Pb (t); Described formula of impedance is suc as formula shown in (1) and the formula (2):
Pf(t)=（p(t)+Z
c×q(t)）÷2 （1）
Pb(t)=（p(t)-Z
c×q(t)）÷2 （2）；
The normalized crosscorrelation coefficient of f, calculating forward-wave Pf (t) and echo Pb (t), the maximum time corresponding point of described normalized crosscorrelation coefficient is TR2, then aortic pulse wave passing time aoPTT is TR2/2.
2. the aortic pulse wave passing time measuring method based on artery of upper extremity information shown in the claim 1, it is characterized in that: described generalized transfer function model is: the Laplace transformation RP (S) of pull-type conversion AP (the S)/radial artery pulse wave of generalized transfer function GTF (S)=aortic pulse wave.
3. the aortic pulse wave passing time measuring method based on artery of upper extremity information shown in the claim 1, it is characterized in that: described blood flow transfer function model is: the pull-type conversion AB (S) of blood flow transfer function BTF (S)=aortic flow waveform/the refer to pull-type conversion FB (S) of artery blood flow waveform.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201310241527.9A CN103284703B (en) | 2013-06-18 | 2013-06-18 | Aortic pulse wave transfer time measuring method based on upper extremity artery information |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201310241527.9A CN103284703B (en) | 2013-06-18 | 2013-06-18 | Aortic pulse wave transfer time measuring method based on upper extremity artery information |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN103284703A true CN103284703A (en) | 2013-09-11 |
| CN103284703B CN103284703B (en) | 2015-05-13 |
Family
ID=49086645
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201310241527.9A Active CN103284703B (en) | 2013-06-18 | 2013-06-18 | Aortic pulse wave transfer time measuring method based on upper extremity artery information |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN103284703B (en) |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104983412A (en) * | 2015-05-28 | 2015-10-21 | 中国科学院合肥物质科学研究院 | Central pulse systole average normalization blood flow waveform model and method for obtaining aorta pulse wave transmission time based on same |
| CN108236459A (en) * | 2016-12-26 | 2018-07-03 | 深圳先进技术研究院 | A kind of central aortic blood pressure measuring method and system |
| CN108464837A (en) * | 2018-05-04 | 2018-08-31 | 重庆大学 | A kind of blood viscosity optics diagnosis and treatment integrated instrument and method towards domestic medicine |
| CN109602401A (en) * | 2019-01-28 | 2019-04-12 | 徐州市心血管病研究所 | A kind of microvascular blood flow kinetic parameter analyzer and analysis method |
| WO2022110020A1 (en) * | 2020-11-26 | 2022-06-02 | 苏州润迈德医疗科技有限公司 | Method for obtaining aortic pressure curve of interest, and storage medium |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108937896B (en) * | 2017-05-23 | 2021-08-03 | 深圳先进技术研究院 | Central arterial blood pressure measuring method and equipment |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2003079586A (en) * | 2001-07-02 | 2003-03-18 | Nippon Colin Co Ltd | Arteriosclerosis evaluation apparatus |
| JP2006102250A (en) * | 2004-10-06 | 2006-04-20 | Terumo Corp | Circulation index measuring device, circulation index measuring method, control program, and computer-readable storage medium |
| CN101084828A (en) * | 2006-06-05 | 2007-12-12 | 中国科学院力学研究所 | Device and method for measuring pulse waving speed |
| CN101317756A (en) * | 2007-06-06 | 2008-12-10 | 索尼株式会社 | Apparatus for obtaining pulse wave velocity information and method thereof |
| CN101810474A (en) * | 2010-01-22 | 2010-08-25 | 徐亮禹 | Device for realizing no-invasive measurement of blood pressure and pulse wave velocity |
-
2013
- 2013-06-18 CN CN201310241527.9A patent/CN103284703B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2003079586A (en) * | 2001-07-02 | 2003-03-18 | Nippon Colin Co Ltd | Arteriosclerosis evaluation apparatus |
| JP2006102250A (en) * | 2004-10-06 | 2006-04-20 | Terumo Corp | Circulation index measuring device, circulation index measuring method, control program, and computer-readable storage medium |
| CN101084828A (en) * | 2006-06-05 | 2007-12-12 | 中国科学院力学研究所 | Device and method for measuring pulse waving speed |
| CN101317756A (en) * | 2007-06-06 | 2008-12-10 | 索尼株式会社 | Apparatus for obtaining pulse wave velocity information and method thereof |
| CN101810474A (en) * | 2010-01-22 | 2010-08-25 | 徐亮禹 | Device for realizing no-invasive measurement of blood pressure and pulse wave velocity |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104983412A (en) * | 2015-05-28 | 2015-10-21 | 中国科学院合肥物质科学研究院 | Central pulse systole average normalization blood flow waveform model and method for obtaining aorta pulse wave transmission time based on same |
| CN108236459A (en) * | 2016-12-26 | 2018-07-03 | 深圳先进技术研究院 | A kind of central aortic blood pressure measuring method and system |
| CN108464837A (en) * | 2018-05-04 | 2018-08-31 | 重庆大学 | A kind of blood viscosity optics diagnosis and treatment integrated instrument and method towards domestic medicine |
| CN109602401A (en) * | 2019-01-28 | 2019-04-12 | 徐州市心血管病研究所 | A kind of microvascular blood flow kinetic parameter analyzer and analysis method |
| CN109602401B (en) * | 2019-01-28 | 2021-08-27 | 徐州市心血管病研究所 | Microvascular hemodynamic parameter analyzer and analysis method |
| WO2022110020A1 (en) * | 2020-11-26 | 2022-06-02 | 苏州润迈德医疗科技有限公司 | Method for obtaining aortic pressure curve of interest, and storage medium |
Also Published As
| Publication number | Publication date |
|---|---|
| CN103284703B (en) | 2015-05-13 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN103284703B (en) | Aortic pulse wave transfer time measuring method based on upper extremity artery information | |
| CN106037694B (en) | A kind of continuous blood pressure measurer based on pulse wave | |
| CN103385702B (en) | A kind of non-invasive blood pressure continuous detection apparatus and method | |
| CN101732040B (en) | Non-invasive multipath pulse wave detection device, system and analytical system | |
| CN105726000B (en) | A kind of device of the cardiovascular functional parameter based on blood pressure of four limbs pulse | |
| CN103054562B (en) | Cardiovascular function detection method based on multi-channel pulse wave form analysis and device thereof | |
| CN104188637B (en) | A kind of aortic pulse wave conduction time acquisition methods based on Waveform Matching method | |
| CN109044302A (en) | Measure device, electronic equipment and the computer readable storage medium of blood pressure | |
| CN104757959B (en) | Pulse wave transmission velocity detecting method and system based on image foldover | |
| CN104068841B (en) | A kind of measuring method and device measuring Indices of Systolic Time parameter | |
| CN101919704B (en) | Heart sound signal positioning and segmenting method | |
| CN204121004U (en) | A kind of pulse diagnostic information acquisition device in conjunction with multiple pressure-type sensor | |
| WO2020182010A1 (en) | Method for determining pulse transmission time, and arteriosclerosis detection device and system | |
| CN103876723A (en) | Method for obtaining blood pressure value by calculating pulse wave conduction time through non-invasive radial artery waves | |
| CN108354597A (en) | A kind of rapid blood pressure computational methods based on the extraction of optimal wave | |
| CN103610452B (en) | Non-contact magnetic induction type pulse detection method | |
| CN104983412B (en) | A kind of central aortic systolic phase average normalized blood flow waveform model and the method that aortic pulse wave passing time is obtained based on it | |
| CN109009004A (en) | A kind of physical examinations method based on Chinese medicine pulse analysis | |
| CN102334986A (en) | Method for identifying dicrotic notch point in arterial tension signal | |
| CN105167759A (en) | Human pulse wave velocity measuring method and system based on intelligent mobile phone | |
| Hu et al. | Feature extraction and choice in PCG based on Hilbert Transfer | |
| CN105708441A (en) | Wearable finger stall and electrocardiogram and pulse condition acquisition terminal | |
| CN104027097A (en) | Vascular function noninvasive detecting method and device | |
| Thakker et al. | Outlier pulse detection and feature extraction for wrist pulse analysis | |
| Almeida et al. | Hemodynamic features extraction from a new arterial pressure waveform probe |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| C14 | Grant of patent or utility model | ||
| GR01 | Patent grant |