CN111513689A - Method for identifying aorta opening characteristic points in pulse waves - Google Patents
Method for identifying aorta opening characteristic points in pulse waves Download PDFInfo
- Publication number
- CN111513689A CN111513689A CN202010463628.0A CN202010463628A CN111513689A CN 111513689 A CN111513689 A CN 111513689A CN 202010463628 A CN202010463628 A CN 202010463628A CN 111513689 A CN111513689 A CN 111513689A
- Authority
- CN
- China
- Prior art keywords
- point
- signal
- order
- value
- pulse
- 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.)
- Pending
Links
- 210000000709 aorta Anatomy 0.000 title claims abstract description 20
- 238000000034 method Methods 0.000 title claims abstract description 17
- 238000005070 sampling Methods 0.000 claims abstract description 7
- 238000001914 filtration Methods 0.000 claims abstract description 6
- 238000004364 calculation method Methods 0.000 claims abstract description 5
- 230000000737 periodic effect Effects 0.000 claims description 10
- 230000009286 beneficial effect Effects 0.000 abstract description 5
- 238000000605 extraction Methods 0.000 abstract description 4
- 238000003745 diagnosis Methods 0.000 description 7
- 210000001367 artery Anatomy 0.000 description 4
- 230000001174 ascending effect Effects 0.000 description 3
- 239000003814 drug Substances 0.000 description 3
- 238000012360 testing method Methods 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 239000008280 blood Substances 0.000 description 2
- 210000004369 blood Anatomy 0.000 description 2
- 201000010099 disease Diseases 0.000 description 2
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 2
- 230000002526 effect on cardiovascular system Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000000717 retained effect Effects 0.000 description 2
- 230000002861 ventricular Effects 0.000 description 2
- 206010003210 Arteriosclerosis Diseases 0.000 description 1
- 208000024172 Cardiovascular disease Diseases 0.000 description 1
- 206010020772 Hypertension Diseases 0.000 description 1
- 210000001765 aortic valve Anatomy 0.000 description 1
- 208000011775 arteriosclerosis disease Diseases 0.000 description 1
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 description 1
- 230000036772 blood pressure Effects 0.000 description 1
- 210000005242 cardiac chamber Anatomy 0.000 description 1
- 230000000747 cardiac effect Effects 0.000 description 1
- 230000008602 contraction Effects 0.000 description 1
- 208000029078 coronary artery disease Diseases 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 230000002706 hydrostatic effect Effects 0.000 description 1
- 210000005240 left ventricle Anatomy 0.000 description 1
- 208000010125 myocardial infarction Diseases 0.000 description 1
- 229910052760 oxygen Inorganic materials 0.000 description 1
- 239000001301 oxygen Substances 0.000 description 1
Images
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
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/48—Other medical applications
- A61B5/4854—Diagnosis based on concepts of traditional oriental medicine
-
- 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
-
- 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
- A61B5/725—Details of waveform analysis using specific filters therefor, e.g. Kalman or adaptive filters
Abstract
The invention discloses a method for identifying aorta opening characteristic points in pulse waves, which comprises the following steps of carrying out low-pass filtering on first-order backward difference of pulse wave signals, calculating second-order difference signals, calculating a main wave peak point and a period of the pulse waves by using the filtered first-order and second-order difference signals, and further identifying aorta opening characteristic points, wherein the identification method comprises the following steps of: searching a right zero point corresponding to the maximum value of the first-order differential data in a period interval 10% before each main wave point; in the signal with lower sampling rate, the position of the minimum value is searched at the left and right positions of the identified right zero point to be used as the aorta opening characteristic point. The innovative identification step of the aorta opening point effectively reduces the complexity of calculation, thereby reducing the requirement on processing equipment and being more beneficial to the extraction of the characteristic point.
Description
Technical Field
The invention relates to a signal processing technology and medical signal processing, in particular to a method for identifying an aorta opening characteristic point in a pulse wave.
Background
The periodic contractions and relaxations of the heart form a pulse wave. Blood in the heart chambers of the human body is injected into the aorta and propagates in the arteries of the human body in the form of waves, namely, pulse waves. The pulse wave has rich cardiovascular information of human body, and can reflect heart rate, cardiovascular state and the like. At present, the pulse wave can be used for diagnosis and treatment of cardiovascular diseases, including determination of coronary heart disease, hypertension, myocardial infarction and the like, and also can be used for psychological state evaluation, body state diagnosis and the like, including noninvasive blood oxygen saturation test, fatigue driving test and the like. The extraction of parameters in the pulse waves is the basis of a series of applications of the pulse waves, and the identification of the characteristic points of the pulse waves is taken as one of key technologies, so that the method has a key significance in noninvasive blood pressure testing and arteriosclerosis.
The pulse wave generally includes 6 parts, ascending branch, descending branch, main wave, tidal wave, descending central isthmus and dicrotic wave, and its physiological significance corresponds to the following: ascending branches correspond to the rapid ejection period of ventricles, descending branches correspond to the period from the later period of ventricular ejection to the beginning of the next cardiac cycle, the main wave is the maximum value of a pulse wave and represents the maximum values of pressure and volume in an artery, and a tidal wave is before the main wave and a dicrotic wave and represents the conditions of stopping ejection of the left ventricle, expanding the artery, reducing the pressure and reversely reflecting the waves; the tidal wave is the aortic hydrostatic emptying time and is the boundary point of systole and diastole; the dicrotic wave is the aortic valve closure, aortic elastic recoil wave. Among these, the start point of the ascending branch is the start point of the pulse wave as the start of the rapid ventricular ejection period, and the identification thereof is significant.
Disclosure of Invention
The invention aims to provide a method for identifying an aorta opening characteristic point in a pulse wave.
The technical scheme adopted by the invention is as follows:
an aorta opening characteristic point identification method in pulse waves comprises the following steps: calculating a first-order difference signal of the original pulse wave signal, performing low-pass filtering on the first-order backward difference signal to obtain a signal serving as a first-order difference signal, and calculating a second-order difference signal; identifying a main wave peak point by using a first-order difference signal and a second-order difference signal, and removing an error point in the main wave peak point; carrying out backward difference on the dominant wave peak to obtain a periodic sequence, and taking a value between quantiles [ Q1 and Q3] in the sequence to calculate a periodic mean value T for reducing errors, wherein Q1 is a first quartile and Q3 is a third quartile in the periodic sequence; calculating an aorta opening characteristic point: and searching a right zero point corresponding to the maximum value of the first-order difference data in a period interval 10% before each main wave point, and searching positions of minimum values at two points on the left and right of the identified right zero point as main artery opening characteristic points for signals with lower sampling rates.
By adopting the technical scheme, the invention has the beneficial effects that: 1. the acquired original pulse wave signals are not filtered, so that the information of the acquired signals is retained to the maximum extent, and the possibility of follow-up disease tracking and other analysis is created; 2. the innovative identification step of the aorta opening point effectively reduces the complexity of calculation, thereby reducing the requirement on processing equipment and being more beneficial to the extraction of characteristic points; meanwhile, the pulse diagnosis device is convenient to wear, and the later-stage commercialization and popularization are facilitated.
Drawings
The invention is described in further detail below with reference to the accompanying drawings and the detailed description;
FIG. 1 is a schematic flow chart of a method for identifying aortic opening feature points in a pulse wave according to the present invention;
fig. 2 is a schematic diagram of an identification result of the method for identifying an aortic opening feature point in a pulse wave according to the present invention.
Detailed Description
To make the objects, technical solutions and advantages of the embodiments of the present application clearer, the technical solutions of the embodiments of the present application will be clearly and completely described below with reference to the drawings in the embodiments of the present application
In order to make the objects, technical solutions and advantages of the present invention more apparent, the present invention is further described in detail below with reference to the accompanying drawings and examples.
The system used by the invention is not limited, and can be applied to a traditional Chinese medicine pulse diagnosis instrument based on a pressure sensor, a traditional Chinese medicine pulse diagnosis instrument based on photoelectric volume pulse waves and the like, and the signal sampling rate is also not required. Meanwhile, the method has the advantages that the use platform is not limited, and the multiple languages are realized.
As shown in fig. 1 or fig. 2, the method for identifying the aortic opening feature point in a pulse wave according to the present invention is described in detail by using a traditional Chinese medicine pulse diagnosis apparatus using a pressure sensor to acquire pulse wave data with a sampling frequency of 100 Hz:
1) calculating a first order backward difference and a second order backward difference of the original pulse signal:
first order difference:
Diff[t]＝OriSignal[t+1]-OriSignal[t]
wherein, t is 1,2,3.. n, n is the number of points of the original pulse signal, Diff is a first-order backward difference signal, and OriSignal is the original pulse wave signal;
and marking a signal obtained by low-pass filtering the first-order backward difference signal Diff as firDiff, carrying out backward difference on the filtered signal, and calculating a second-order difference signal:
secDiff[t]＝firDiff[t+1]-firDiff[t]
2) utilizing the first order difference signal and the second order difference signal to identify the main peak point of the pulse wave:
and performing threshold processing on the firDiff, setting the part above the threshold value as 1, and setting the part below the threshold value as 0:
wherein the threshold value threshold is 0.5 (median) (firdiff)
Calculating the main wave peak point of the signal:
l:where(diff(SfirDiff·SsecDiff))＝＝-2
wherein point l is the main wave peak point of the signal, SsecDiffSymbol value for secDiff:
the higher derivative tide point is removed in point l:
extracting the amplitude OriSignal (l) of the original signal at the position of the point l, removing the point h with the median value outside three standard deviations from the OriSignal (l), and obtaining the main wave peak point i ═ l-h
3) Calculating a period from the identified main point:
and carrying out backward difference on the main wave peak to obtain a periodic sequence, and taking values between quantiles [ Q1 and Q3] in the sequence to calculate a periodic mean value T, wherein Q1 is a first quartile and Q3 is a third quartile in the periodic sequence.
4) Calculating the aorta opening point:
and searching a right zero point corresponding to the maximum value of the first-order difference data in a period interval 10% before each main wave point, and searching positions of minimum values at two points on the left and right of the identified right zero point in a signal with a low sampling rate to serve as aorta opening characteristic points.
Compared with the prior art, the invention has the beneficial effects that: 1. the acquired original pulse wave signals are not filtered, so that the information of the acquired signals is retained to the maximum extent, and the possibility of follow-up disease tracking and other analysis is created; 2. the innovative identification step of the aorta opening point effectively reduces the complexity of calculation, thereby reducing the requirement on processing equipment and being more beneficial to the extraction of characteristic points; meanwhile, the pulse diagnosis device is convenient to wear, and the later-stage commercialization and popularization are facilitated.
It is to be understood that the embodiments described are only a few embodiments of the present application and not all embodiments. The embodiments and features of the embodiments in the present application may be combined with each other without conflict. The components of the embodiments of the present application, generally described and illustrated in the figures herein, can be arranged and designed in a wide variety of different configurations. Thus, the detailed description of the embodiments of the present application is not intended to limit the scope of the claimed application, but is merely representative of selected embodiments of the application. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present application.
Claims (4)
1. A method for identifying aorta opening characteristic points in pulse waves is characterized by comprising the following steps: which comprises the following steps:
s1, difference processing: after a first-order posterior difference signal of the signal is obtained, filtering is carried out, and after the first-order difference after filtering is filtered, a second-order posterior difference is obtained;
s2, dominant wave crest and period identification: calculating pulse wave main wave peak point and period based on first order difference signal and second order difference signal
S3, aorta opening point identification: searching a right zero point corresponding to the maximum value of the first-order differential data in a period interval 10% before each main wave point; in the signal with lower sampling rate, the position of the minimum value is searched at the left and right positions of the identified right zero point to be used as the aorta opening characteristic point.
2. The method for identifying aortic opening feature points in pulse waves as claimed in claim 1, wherein: the calculation process of the first order difference signal and the second order difference signal in S1 is as follows:
calculate the first order difference of the original pulse signal:
Diff[t]＝OriSignal[t+1]-OriSignal[t]
wherein, t is 1,2,3.. n, n is the number of points of the original pulse signal, Diff is a first-order backward difference signal, and OriSignal is the original pulse wave signal;
and marking a signal obtained by low-pass filtering the first-order backward difference signal Diff as firDiff, carrying out backward difference on the filtered signal, and calculating a second-order difference signal:
secDiff[t]＝firDiff[t+1]-firDiff[t]。
3. the method for identifying aortic opening feature points in pulse waves as claimed in claim 2, wherein: the calculation steps of the period of the pulse wave and the main wave peak point in the S2 are as follows:
s2-1, identifying a main peak point:
s2-1-1, performing threshold processing on the firDiff, setting the part above the threshold value as 1, setting the part below the threshold value as 0:
wherein the threshold value threshold is 0.5 (median) (firdiff)
S2-1-2, calculating the main wave peak point of the signal,
l:where(diff(SfirDiff·SsecDiff))＝＝-2，
wherein point l is the main wave peak point of the signal, SsecDiffIs the symbol value of secDiff,
the tide point with larger derivative value is removed from the point l, and the steps are as follows: extracting the amplitude OriSignal (l) of the original signal at the position of the point l, removing the point h with the median value outside three standard deviations from the OriSignal (l), and obtaining the main wave peak point i ═ l-h
S2-2, calculating a period value:
and carrying out backward difference on the main wave peak to obtain a periodic sequence, and taking a value between a first quartile Q1 and a third quartile Q3 quantile in the periodic sequence to calculate a periodic mean value T so as to reduce errors.
4. The method for identifying the pulse wave aorta opening feature point according to claim 1, wherein: in S2, in the signal with the lower sampling rate, the position of the minimum value is found at the two left and right positions of the identified right zero point as the aorta opening feature point.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010463628.0A CN111513689A (en) | 2020-05-27 | 2020-05-27 | Method for identifying aorta opening characteristic points in pulse waves |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010463628.0A CN111513689A (en) | 2020-05-27 | 2020-05-27 | Method for identifying aorta opening characteristic points in pulse waves |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN111513689A true CN111513689A (en) | 2020-08-11 |
Family
ID=71908187
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010463628.0A Pending CN111513689A (en) | 2020-05-27 | 2020-05-27 | Method for identifying aorta opening characteristic points in pulse waves |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111513689A (en) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112244802A (en) * | 2020-10-16 | 2021-01-22 | 山东省人工智能研究院 | BCG signal heart rate extraction method based on Stockwell transformation |
| CN114469018A (en) * | 2022-02-22 | 2022-05-13 | 平安科技（深圳）有限公司 | Pulse-taking fingering recommendation method and device, electronic equipment and storage medium |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6331162B1 (en) * | 1999-02-01 | 2001-12-18 | Gary F. Mitchell | Pulse wave velocity measuring device |
| CN101732033A (en) * | 2008-11-07 | 2010-06-16 | 中国科学院计算技术研究所 | Method and device for extracting characteristic parameter in human body waveform |
| CN104161510A (en) * | 2014-08-27 | 2014-11-26 | 电子科技大学 | Multistage lead electrocardiograph signal QRS waveform identification method |
| CN107822608A (en) * | 2017-10-26 | 2018-03-23 | 中国民航大学 | Pulse wave feature extracting method based on gauss hybrid models |
| CN107898443A (en) * | 2017-11-21 | 2018-04-13 | 深圳先进技术研究院 | Dicrotic pulse wave detecting method, device and computer-readable storage medium |
-
2020
- 2020-05-27 CN CN202010463628.0A patent/CN111513689A/en active Pending
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6331162B1 (en) * | 1999-02-01 | 2001-12-18 | Gary F. Mitchell | Pulse wave velocity measuring device |
| CN101732033A (en) * | 2008-11-07 | 2010-06-16 | 中国科学院计算技术研究所 | Method and device for extracting characteristic parameter in human body waveform |
| CN104161510A (en) * | 2014-08-27 | 2014-11-26 | 电子科技大学 | Multistage lead electrocardiograph signal QRS waveform identification method |
| CN107822608A (en) * | 2017-10-26 | 2018-03-23 | 中国民航大学 | Pulse wave feature extracting method based on gauss hybrid models |
| CN107898443A (en) * | 2017-11-21 | 2018-04-13 | 深圳先进技术研究院 | Dicrotic pulse wave detecting method, device and computer-readable storage medium |
Non-Patent Citations (1)
| Title |
|---|
| RONNY: "图像分析：投影曲线的波峰查找", 《思维之际》 * |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112244802A (en) * | 2020-10-16 | 2021-01-22 | 山东省人工智能研究院 | BCG signal heart rate extraction method based on Stockwell transformation |
| CN114469018A (en) * | 2022-02-22 | 2022-05-13 | 平安科技（深圳）有限公司 | Pulse-taking fingering recommendation method and device, electronic equipment and storage medium |
| CN114469018B (en) * | 2022-02-22 | 2023-07-21 | 平安科技（深圳）有限公司 | Pulse diagnosis fingering recommendation method and device, electronic equipment and storage medium |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| WO2021208490A1 (en) | Blood pressure measuring method and device based on deep neural network | |
| CN106974631B (en) | Blood pressure measuring method and device based on pulse wave waveform and electrocardiosignal | |
| Kurylyak et al. | A Neural Network-based method for continuous blood pressure estimation from a PPG signal | |
| CN107432736B (en) | Method for identifying pulse waveform signal | |
| CN103479383B (en) | Device for analyzing heart sound signals, and intelligent heart stethoscope provided with device for analyzing heart sound signals | |
| EP1744666A1 (en) | Non-invasive measurement of second heart sound components | |
| CN101732033A (en) | Method and device for extracting characteristic parameter in human body waveform | |
| CN101919704B (en) | Heart sound signal positioning and segmenting method | |
| Argha et al. | A novel automated blood pressure estimation algorithm using sequences of Korotkoff sounds | |
| CN111528821A (en) | Method for identifying characteristic points of counterpulsation waves in pulse waves | |
| CN111513689A (en) | Method for identifying aorta opening characteristic points in pulse waves | |
| CN106073783A (en) | A kind of method extracting breathing rate from photoplethysmographic | |
| CN110840428B (en) | Noninvasive blood pressure estimation method based on one-dimensional U-Net network | |
| CN104068841B (en) | A kind of measuring method and device measuring Indices of Systolic Time parameter | |
| Foroozan et al. | Robust beat-to-beat detection algorithm for pulse rate variability analysis from wrist photoplethysmography signals | |
| Yang et al. | Removal of pulse waveform baseline drift using cubic spline interpolation | |
| CN113243902B (en) | Feature extraction method based on photoplethysmography | |
| RU2552685C2 (en) | Method, device and programme for automatic processing of blood pressure signals | |
| Pan et al. | Improved blood pressure estimation using photoplethysmography based on ensemble method | |
| CN108814642A (en) | A kind of heart sound location and rate calculation method of electronic auscultation device | |
| CN112587104A (en) | Method for filtering invalid pulse waveform | |
| Hikmah et al. | A signal processing framework for multimodal cardiac analysis | |
| Hoseinzadeh et al. | Design and Implementation of a blood pressure device with high sampling frequency to analyze cardiovascular diseases in LabVIEW | |
| CN110123369A (en) | A kind of blood pressure quick calculation method based on multichannel cardiechema signals | |
| Chen et al. | Automatic Onsets and Systolic Peaks Detection and Segmentation of Arterial Blood Pressure Waveforms using Fully Convolutional Neural Networks |
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 | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20200811 |
|
| RJ01 | Rejection of invention patent application after publication |