CN114601440A - Impedance rheogram processing method for blood pressure detection - Google Patents
Impedance rheogram processing method for blood pressure detectionDownload PDFInfo
- Publication number
- CN114601440A CN114601440ACN202210244270.1ACN202210244270ACN114601440ACN 114601440 ACN114601440 ACN 114601440ACN 202210244270 ACN202210244270 ACN 202210244270ACN 114601440 ACN114601440 ACN 114601440A
- Authority
- CN
- China
- Prior art keywords
- impedance
- point
- interval
- wave
- blood flow
- 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
- 230000036772blood pressureEffects0.000titleclaimsabstractdescription33
- 238000001514detection methodMethods0.000titleclaimsabstractdescription15
- 238000003672processing methodMethods0.000titleabstractdescription7
- 230000017531blood circulationEffects0.000claimsabstractdescription52
- 230000000737periodic effectEffects0.000claimsabstractdescription21
- 238000000034methodMethods0.000claimsabstractdescription20
- 230000011218segmentationEffects0.000claimsabstractdescription7
- 238000013519translationMethods0.000claimsabstract2
- 238000010586diagramMethods0.000claimsdescription34
- 230000010363phase shiftEffects0.000claimsdescription9
- 230000035487diastolic blood pressureEffects0.000claimsdescription7
- 230000035488systolic blood pressureEffects0.000claimsdescription7
- 238000004364calculation methodMethods0.000claimsdescription3
- 238000012417linear regressionMethods0.000claimsdescription3
- 238000012545processingMethods0.000claimsdescription3
- 230000001419dependent effectEffects0.000claimsdescription2
- 238000012549trainingMethods0.000claimsdescription2
- 238000002372labellingMethods0.000claims1
- 238000012544monitoring processMethods0.000claims1
- 239000000523sampleSubstances0.000description23
- 238000005259measurementMethods0.000description12
- 238000000605extractionMethods0.000description11
- 210000004369bloodAnatomy0.000description4
- 239000008280bloodSubstances0.000description4
- 238000009530blood pressure measurementMethods0.000description4
- 230000000747cardiac effectEffects0.000description3
- 230000000694effectsEffects0.000description3
- 239000000284extractSubstances0.000description3
- 230000000004hemodynamic effectEffects0.000description3
- 230000002861ventricularEffects0.000description3
- 238000005086pumpingMethods0.000description2
- 230000000630rising effectEffects0.000description2
- 230000002159abnormal effectEffects0.000description1
- 230000003321amplificationEffects0.000description1
- 238000004458analytical methodMethods0.000description1
- 230000001746atrial effectEffects0.000description1
- 230000009286beneficial effectEffects0.000description1
- 238000005516engineering processMethods0.000description1
- 238000001914filtrationMethods0.000description1
- 230000006870functionEffects0.000description1
- 230000004217heart functionEffects0.000description1
- 238000003199nucleic acid amplification methodMethods0.000description1
- 238000005192partitionMethods0.000description1
- 210000005259peripheral bloodAnatomy0.000description1
- 239000011886peripheral bloodSubstances0.000description1
- 238000013186photoplethysmographyMethods0.000description1
- 238000011160researchMethods0.000description1
- 230000000717retained effectEffects0.000description1
- 230000001629suppressionEffects0.000description1
- 230000001360synchronised effectEffects0.000description1
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 for evaluating the cardiovascular system, e.g. pulse, heart rate, blood pressure or blood flow
- A61B5/026—Measuring blood flow
- A61B5/029—Measuring blood output from the heart, e.g. minute volume
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/0033—Features or image-related aspects of imaging apparatus, e.g. for MRI, optical tomography or impedance tomography apparatus; Arrangements of imaging apparatus in a room
- A61B5/004—Features or image-related aspects of imaging apparatus, e.g. for MRI, optical tomography or impedance tomography apparatus; Arrangements of imaging apparatus in a room adapted for image acquisition of a particular organ or body part
- A61B5/0044—Features or image-related aspects of imaging apparatus, e.g. for MRI, optical tomography or impedance tomography apparatus; Arrangements of imaging apparatus in a room adapted for image acquisition of a particular organ or body part for the heart
- 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 for evaluating the cardiovascular system, e.g. pulse, heart rate, blood pressure or blood flow
- A61B5/021—Measuring pressure in heart or blood vessels
- 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 for evaluating the cardiovascular system, e.g. pulse, heart rate, blood pressure or blood flow
- A61B5/026—Measuring blood flow
- A61B5/0295—Measuring blood flow using plethysmography, i.e. measuring the variations in the volume of a body part as modified by the circulation of blood therethrough, e.g. impedance plethysmography
- 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
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Cardiology (AREA)
- Engineering & Computer Science (AREA)
- Molecular Biology (AREA)
- General Health & Medical Sciences (AREA)
- Biophysics (AREA)
- Pathology (AREA)
- Veterinary Medicine (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Medical Informatics (AREA)
- Public Health (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- Physics & Mathematics (AREA)
- Physiology (AREA)
- Hematology (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Radiology & Medical Imaging (AREA)
- Vascular Medicine (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Psychiatry (AREA)
- Signal Processing (AREA)
- Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)
Abstract
Description
Translated fromChinese技术领域technical field
本发明专利涉及医疗器械技术领域，尤其涉及一种用于血压检测的阻抗血流图处理方法。The patent of the present invention relates to the technical field of medical devices, in particular to an impedance flow diagram processing method for blood pressure detection.
背景技术Background technique
阻抗血流图是一种用电生物阻抗技术测定心输出量，评判心脏功能的无创性方法。其通过测量随血容量变化而变化的电阻抗来测定特定身体部分的血容量的变化。阻抗血流图作为一种容积描记法，随着每个心动周期引起的血容量的增减，利用有规律性阻抗变化的容积变化曲线，通过识别相关特征点获得心脏输出的相关血流动力学参数。Impedance flow diagram is a non-invasive method to measure cardiac output and evaluate cardiac function by electrical bioimpedance technique. It measures changes in blood volume in a particular body part by measuring the electrical impedance as a function of blood volume. As a plethysmography method, impedance flow diagram uses the volume change curve of regular impedance change with the increase or decrease of blood volume caused by each cardiac cycle, and obtains the relevant hemodynamics of cardiac output by identifying relevant characteristic points. parameter.
血压是血流动力学参数中一项重要的指标，血压能体现心脏的泵血能力、周围血管的阻力和弹性，日常生活中备受人们关注。现有的血压测量技术主要分为侵入式测血压和非侵入式测血压两种方式，侵入式测量血压存在测量不便，易感染和不能进行长时间连续监测的缺点，而针对无创测量血压方式，国内外大都采用光电容积描记法检测手臂上的脉搏波，获取脉搏波传导时间，通过脉搏波传导时间提取血压参数，血压参数与实测存在一定误差，且检测准确度仍有待提高。Blood pressure is an important indicator of hemodynamic parameters. Blood pressure can reflect the pumping ability of the heart, the resistance and elasticity of peripheral blood vessels, and it has attracted much attention in daily life. Existing blood pressure measurement technologies are mainly divided into invasive blood pressure measurement and non-invasive blood pressure measurement. At home and abroad, photoplethysmography is mostly used to detect the pulse wave on the arm, obtain the pulse wave transit time, and extract the blood pressure parameters through the pulse wave transit time. There is a certain error between the blood pressure parameters and the actual measurement, and the detection accuracy still needs to be improved.
阻抗血流图反映了动脉血管的容积变化，体现心脏的泵血趋势。阻抗血流图存在快速射血期和缓慢射血期，其血流特征点与血压参数有着紧密联系。现阶段，国内大部分研究者针对阻抗血流图的血压检测研究并不多，人们在进行检测时不仅存在个体性的差异，还会不可避免的因为测量抖动产生噪声，伪影，进而影响了特征点提取的准确性，造成血流动力学参数的误差。另外国内外研究者采用同步采集阻抗血流信号与心电信号的形式，利用两者在时域上的对应关系，通过识别心电信号特征点来辅助血流信号特征点的识别，且得到了较高的识别率。然而同时采集心电图和阻抗血流图，两者间必然会存在一定的相位差，容易造成测量误差，同时这样也增加了阻抗血流图设备的复杂性，减少了测量的便捷性。The impedance flow map reflects the volume change of the arterial vessels and reflects the pumping trend of the heart. There are fast ejection and slow ejection phases in impedance flow diagram, and the characteristic points of blood flow are closely related to blood pressure parameters. At this stage, most domestic researchers do not have much research on blood pressure detection with impedance blood flow diagrams. People not only have individual differences in the detection, but also inevitably generate noise and artifacts due to measurement jitter, which in turn affects the The accuracy of feature point extraction results in errors in hemodynamic parameters. In addition, domestic and foreign researchers adopted the form of synchronous acquisition of impedance blood flow signal and ECG signal, and used the corresponding relationship between the two in the time domain to assist the identification of blood flow signal feature points by identifying ECG signal feature points, and obtained higher recognition rate. However, when the electrocardiogram and the impedance flow chart are collected at the same time, there will inevitably be a certain phase difference between the two, which is likely to cause measurement errors. At the same time, this also increases the complexity of the impedance flow chart equipment and reduces the convenience of measurement.
一个阻抗血流信号的周期是从A波(房缩波)开始，到O波(室舒波)之后结束。国内研究者在进行阻抗血流图的样本划分时，划分完不能直观体现一个血流信号周期的始末这样不利于对特征点的识别提取。同时忽略了阻抗血流图中A波和O波对样本选择的影响以及特征点B点位置偏离和如何准确识别的问题。虽然说A波和O波不是作为特征点提取的参考波段，但是周期样本是一个整体，倘若这两个波段不正常容易对特征点的提取产生影响。当两个波峰值比较高的时候也会有可能是因为测量时抖动所引起的噪声变化。大部分人在对特征点B点定位的时候，采用临床上C波高度的15％来定位该点的方式，但是在实际测量过程中不可避免产生噪声干扰，检测出来的波形不是稳定的曲线，会出现多个符合这个条件的点，因此使用该方法容易造成识别B点错误，导致特征点B点位置的偏移。B点作为心室开始射血的起始点，代表一种血液喷射的趋势，主要表现为阻抗导数(dZ/dt)快速上升的反转、拐点或者快速斜率变化。但是在实际测量中由于受到噪声和伪影的影响，容易导致dZ/dt出现多个反转、拐点和快速斜率变化，进而影响到B点的识别。当然在实际测量中，也有很多样本不存在这样的趋势，它可能是一条平滑的上升曲线，因此需要分情况判断。针对以上问题，本发明对阻抗血流图进行相位平移，对样本数据分割和分区搜索，约束A波和O波幅值，限制B点位置区间以及考虑BC段出现拐点的情况，实现对阻抗血流图中噪声的抑制，提高血流特征点提取的准确性。The cycle of an impedance blood flow signal starts from the A wave (atrial systolic wave) and ends after the O wave (ventricular systolic wave). When domestic researchers divide the samples of the impedance blood flow map, the division cannot directly reflect the beginning and end of a blood flow signal cycle, which is not conducive to the identification and extraction of feature points. At the same time, the influence of the A wave and O wave on the sample selection in the impedance blood flow diagram and the position deviation of the feature point B and how to accurately identify the problem are ignored. Although A wave and O wave are not used as reference bands for feature point extraction, the periodic sample is a whole, and if these two bands are abnormal, it will easily affect the extraction of feature points. When the two peaks are relatively high, it may also be due to noise variation caused by jitter during measurement. When most people locate the feature point B, they use 15% of the clinical C wave height to locate the point, but in the actual measurement process, noise interference is inevitable, and the detected waveform is not a stable curve. There will be multiple points that meet this condition, so using this method is easy to cause errors in identifying point B, resulting in the offset of the position of point B of the feature point. Point B, as the starting point of ventricular ejection, represents a trend of blood ejection, which is mainly manifested as the inversion of the rapid rise of the impedance derivative (dZ/dt), the inflection point or the rapid slope change. However, due to the influence of noise and artifacts in actual measurement, it is easy to cause multiple inversions, inflection points and rapid slope changes in dZ/dt, which in turn affects the identification of point B. Of course, in actual measurement, there are also many samples that do not have such a trend. It may be a smooth upward curve, so it needs to be judged according to the situation. In view of the above problems, the present invention performs phase shift on the impedance blood flow map, divides and partitions the sample data, constrains the amplitudes of A and O waves, limits the position interval of point B, and considers the inflection point of the BC segment, so as to realize the analysis of the impedance blood flow. The suppression of noise in the flow graph improves the accuracy of the extraction of blood flow feature points.
发明内容SUMMARY OF THE INVENTION
本发明的目的在于提供一种用于血压检测的阻抗血流图处理方法，旨在解决针对阻抗血流图中重要特征点B点(心室射血起始点)需要准确的定位与标记，需要对在使用阻抗血流图设备过程中因为测量抖动引起的噪声进行有效剔除。以及在保证较高识别率的同时，需要减少阻抗血流图设备的复杂性，降低成本，提高测量的便捷性以及通过提取阻抗血流图中的血流特征参数获取血压数值的问题，包括：The purpose of the present invention is to provide an impedance blood flow diagram processing method for blood pressure detection, which aims to solve the need for accurate positioning and marking of the important feature point B (the starting point of ventricular ejection) in the impedance flow diagram. The noise caused by measurement jitter is effectively eliminated in the process of using the impedance flow diagram equipment. And while ensuring a high recognition rate, it is necessary to reduce the complexity of impedance blood flow map equipment, reduce costs, improve the convenience of measurement, and obtain blood pressure values by extracting blood flow characteristic parameters in impedance blood flow maps, including:
步骤1：确定所述阻抗血流图中阻抗微分图的起始点和终止点；Step 1: determine the starting point and the ending point of the impedance differential graph in the impedance blood flow graph;
步骤2：对所述阻抗微分图的周期样本区间进行相位平移；Step 2: performing phase shift on the periodic sample interval of the impedance differential map;
步骤3：对所述阻抗微分图的周期样本进行数据分割和区间搜索；Step 3: perform data segmentation and interval search on the periodic samples of the impedance differential graph;
步骤4：对所述阻抗微分图中的A波和O波进行幅值约束，剔除无效样本，并且重新标记起始点和终止点；Step 4: Constrain the amplitude of the A wave and the O wave in the impedance differential diagram, remove invalid samples, and re-mark the starting point and the ending point;
步骤5：对所述阻抗微分图的特征点B点进行约束，限制B点位置在QC间期，剔除无效样本；Step 5: Constrain the feature point B of the impedance differential map, limit the position of point B in the QC interval, and eliminate invalid samples;
步骤6：判断所述阻抗微分图中QC间期出现拐点的情况，对特征点B，C，X点进行标记。Step 6: Determine the inflection point in the QC interval in the impedance differential diagram, and mark the feature points B, C, and X.
其中，确定所述阻抗血流图中阻抗微分图的起始点和终止点，包括：Wherein, determining the starting point and the ending point of the impedance differential graph in the impedance blood flow graph, including:
利用阻抗变化量波形来辅助阻抗微分量波形划分样本。遍历阻抗血流图中阻抗变化量每个周期的极小值点，标记其在时域上对应于阻抗微分图各个周期的点即为起始点和终止点。The impedance variation waveform is used to assist in dividing the samples by the impedance derivative waveform. Traverse the minimum value points of each cycle of the impedance variation in the impedance blood flow diagram, and mark the points corresponding to each cycle of the impedance differential diagram in the time domain as the start point and the end point.
其中，对所述阻抗微分图的周期样本区间进行相位平移，包括：Wherein, performing phase shift on the periodic sample interval of the impedance differential map includes:
对所述阻抗微分图的周期样本区间在时域上进行相位平移，划分一个完整的阻抗微分图周期在一个样本区间内。Phase shift is performed on the periodic sample interval of the impedance differential map in the time domain, and a complete impedance differential map cycle is divided into one sample interval.
其中，对所述阻抗微分图的周期样本进行数据分割和区间搜索，包括：Wherein, performing data segmentation and interval search on the periodic samples of the impedance differential graph, including:
对所述阻抗微分图进行数据分割，以Q，C，X点为间隔线，将阻抗微分图的QX间期划分为QC间期和CX间期，实现特征点在固定区间进行搜索。The impedance differential map is divided into data, and the Q, C, and X points are used as interval lines, and the QX interval of the impedance differential map is divided into the QC interval and the CX interval, so that the feature points can be searched in a fixed interval.
其中，对所述阻抗微分图中的A波和O波进行幅值约束，剔除无效样本，并且重新标记起始点和终止点，包括：Wherein, amplitude constraints are performed on the A-wave and O-wave in the impedance differential diagram, invalid samples are eliminated, and the start and end points are re-marked, including:
所述阻抗微分图在进行区间划分后，遍历QC间期之前和CX间期之后的极大值为所述阻抗微分图的A波和O波，重新标记A波为所述阻抗微分图的周期起始点，O波为终止点。约束A波和O波的幅值<C波幅值的1/3，使两波的峰值在限制在一个正常的范围内，剔除不符合该条件的无效样本周期，尽可能在一定程度上还原真实准确的周期样本。After the impedance differential map is divided into intervals, the maximum values before traversing the QC interval and after the CX interval are the A waves and O waves of the impedance differential map, and the A wave is re-marked as the period of the impedance differential map. The starting point, the O wave is the ending point. Constrain the amplitude of wave A and wave O to be less than 1/3 of the amplitude of wave C, so that the peaks of the two waves are limited to a normal range, and invalid sample periods that do not meet this condition are eliminated, and restore to a certain extent as much as possible True and accurate cycle samples.
其中，对所述阻抗微分图的特征点B点进行约束，限制B点位置在QC间期，剔除无效样本周期，包括：Among them, the characteristic point B of the impedance differential map is constrained to limit the position of point B to the QC interval, and the invalid sample period is eliminated, including:
遍历C波幅值的15％的点，将其标记为所述阻抗微分图的特征点B点，限制B点位置在QC间期，剔除B点偏离QC间期的无效样本周期。Traverse the 15% point of the C wave amplitude, mark it as the characteristic point B point of the impedance differential map, limit the position of point B to the QC interval, and eliminate the invalid sample period when the point B deviates from the QC interval.
其中，判断所述阻抗微分图中QC间期出现拐点的情况，对特征点B，C，X点进行标记，包括：Wherein, it is judged that the inflection point occurs in the QC interval in the impedance differential diagram, and the characteristic points B, C, and X are marked, including:
判断QC间期是否出现拐点，如果没有则考虑将B点定位为C波振幅的15％，如果出现了拐点，判断出现的拐点的数量，若出现两个或者以上的拐点则判断该区间为无效样本周期，若只出现了一个拐点则判断该点为特征点B点。遍历所述阻抗微分图周期样本的最大值和最小值，分别标记为阻抗微分图的C点和X点。Determine whether there is an inflection point in the QC interval. If not, consider positioning the B point as 15% of the C wave amplitude. If there is an inflection point, determine the number of inflection points that appear. If there are two or more inflection points, the interval is judged to be invalid. Sample period, if there is only one inflection point, it is judged that this point is the feature point B point. Traverse the maximum and minimum values of the periodic samples of the differential impedance graph, marked as points C and X of the differential impedance graph, respectively.
本发明提供一种基于阻抗血流图的血压检测方法，与现阶段常规的阻抗血流图处理方法和无创血压检测方法相比，本发明的有效益果体现在：The present invention provides a blood pressure detection method based on impedance blood flow diagram. Compared with the current conventional impedance flow diagram processing method and non-invasive blood pressure detection method, the present invention has the following beneficial effects:
采用数据分割和区间划分的方式对阻抗微分图的周期样本进行处理，考虑了A波和O波幅值对周期样本的影响，有效剔除因噪声引起的无效样本，提高抗干扰性。The periodic samples of the impedance differential graph are processed by means of data division and interval division, and the influence of A-wave and O-wave amplitudes on periodic samples is considered, which effectively eliminates invalid samples caused by noise and improves anti-interference.
对重要特征点B点进行区间约束和QC间期上拐点数量的判断，排除了在实际测量过程因为噪声和伪影而造成的B点识别错误的情况，提高B点识别的准确率。The interval constraint of important feature point B and the judgment of the number of inflection points on the QC interval eliminates the error of B point recognition caused by noise and artifacts in the actual measurement process, and improves the accuracy of B point recognition.
在保证较高识别率的同时，仅使用阻抗血流图的相关波形进行特征点提取，减少阻抗血流图设备的复杂性，降低成本，提高测量的便捷性。While ensuring a high recognition rate, only the relevant waveforms of the impedance blood flow diagram are used for feature point extraction, which reduces the complexity of the impedance flow diagram equipment, reduces the cost, and improves the convenience of measurement.
利用阻抗血流图的相关血流特征参数建立血压线性回归模型，提高测量血压数值的准确性。The blood pressure linear regression model was established by using the relevant blood flow characteristic parameters of the impedance flow chart to improve the accuracy of blood pressure measurement.
附图说明Description of drawings
图1是本发明的一种阻抗血流图的血压检测方法的流程图；Fig. 1 is a flow chart of a blood pressure detection method of an impedance blood flow diagram of the present invention;
图2是本发明的硬件平台所采集的阻抗血流信号；Fig. 2 is the impedance blood flow signal collected by the hardware platform of the present invention;
图3是本发明的初始样本划分的效果图；Fig. 3 is the effect diagram of the initial sample division of the present invention;
图4是本发明的经过相位平移后的样本划分效果图；Fig. 4 is the sample division effect diagram after phase shift of the present invention;
图5是本发明的数据分割与区间划分效果图；Fig. 5 is the data division and interval division effect diagram of the present invention;
图6是本发明的特征点提取流程图；Fig. 6 is the characteristic point extraction flow chart of the present invention;
具体实施方式Detailed ways
为了使本申请的目的、技术方案及优点更加清楚明白，下面结合附图以及实施例，对本发明进一步详细说明。In order to make the objectives, technical solutions and advantages of the present application clearer, the present invention will be further described in detail below with reference to the accompanying drawings and embodiments.
请参阅图1，本发明提供一种基于阻抗血流图的特征点提取处理方法，包括：Referring to FIG. 1, the present invention provides a feature point extraction processing method based on an impedance blood flow diagram, including:
S101、确定阻抗血流图中阻抗微分图的起始点和终止点。S101. Determine the starting point and the ending point of the impedance differential graph in the impedance blood flow graph.
具体的，参阅图2，如图2所示为本发明的硬件平台所采集的阻抗血流信号，其经过差分放大，检波，滤波以及去基线处理，得到的实测波形与理想波形的相似度较高，符合客观规律，具有参考价值。阻抗血流信号主要包括阻抗变化量和阻抗微分量，其中阻抗变化量周期起始点相对于周期内其他信号点而言，具有更低的幅值，易于区分信号的周期，因此利用阻抗变化量波形来辅助阻抗微分量波形划分样本。实现步骤：通过逐点遍历比较周期内最小值的方式提取阻抗血流图中阻抗变化量每个周期的极小值点，标记其在时域上对应于阻抗微分图各个周期的点即为起始点和终止点。Specifically, referring to FIG. 2 , as shown in FIG. 2 is the impedance blood flow signal collected by the hardware platform of the present invention. After differential amplification, wave detection, filtering and de-baseline processing, the obtained measured waveform is more similar to the ideal waveform. High, in line with objective laws, with reference value. The impedance blood flow signal mainly includes impedance variation and impedance differential. The starting point of the impedance variation cycle has a lower amplitude compared to other signal points in the cycle, which is easy to distinguish the cycle of the signal. Therefore, the impedance variation waveform is used. to assist in dividing the impedance derivative waveform into samples. Implementation steps: Extract the minimum value point of each cycle of the impedance change in the impedance blood flow diagram by traversing and comparing the minimum value in the cycle point by point, and mark the point corresponding to each cycle of the impedance differential diagram in the time domain. start and end points.
S102、对阻抗微分图的周期样本区间进行相位平移。S102. Perform phase shift on the periodic sample interval of the impedance differential graph.
具体的，参阅图3和图4，对样本区间在时域上进行相位平移，设定一个样本移动的阈值，阈值<平均周期样本时间，使得周期样本在阈值前后移动，直到整个完整的阻抗微分图周期划分到样本中。Specifically, referring to Fig. 3 and Fig. 4, the phase shift is performed on the sample interval in the time domain, and a threshold value of sample movement is set, the threshold value < average period sample time, so that the period sample moves before and after the threshold value until the entire impedance differential is complete The graph period is divided into samples.
S103、对所述阻抗微分图的周期样本进行数据分割和区间搜索。S103. Perform data division and interval search on the periodic samples of the impedance differential map.
具体的，参阅图5，为后面更好实现特征点在固定区间进行搜索，提高特征点提取准确性。对阻抗微分图进行数据分割，以Q，C，X点为间隔线，将阻抗微分图的QX间期划分为QC间期和CX间期，通过逐点遍历比较周期内最小值和最大值的方式提取每个周期样本的极大值点和极小值点，即为所需要提取的特征点C点和X点。然后将B至X区间数据以B，X特征点为界进行区间数据分割。在B点前段遍历极小值点，标记该点为Q点，令QB时间区间TQB<1/3*TBC。使用分割线对Q点，B点，X点为时间轴进行区间划分。Specifically, referring to FIG. 5 , in order to better implement the feature point search in a fixed interval later, the feature point extraction accuracy is improved. Divide the data of the impedance differential graph, take Q, C, and X points as the interval lines, divide the QX interval of the impedance differential graph into the QC interval and the CX interval, and compare the minimum and maximum values in the period by traversing point by point. The maximum value point and the minimum value point of each cycle sample are extracted in this way, that is, the feature points C and X that need to be extracted. Then, the interval data from B to X is divided into interval data with B and X feature points as the boundary. Traverse the minimum point before point B, mark this point as point Q, and set the QB time interval TQB <1/3*TBC . Use the dividing line to divide the Q point, B point, and X point as the time axis.
S104、对阻抗微分图中的A波和O波进行幅值约束，剔除无效样本，并且重新标记起始点和终止点。S104: Constrain the amplitude of the A wave and the O wave in the impedance differential diagram, remove invalid samples, and remark the start point and the end point.
具体的，排除因为噪声引起A波和O波幅值偏大的样本，保留与标准样本的共性，提高特征点提取的准确性。通过逐点遍历比较最大值找到以Q点为界限之前的极大值点，即为A波，令A波幅值≤C波幅值的1/3。同理遍历X点界限之后的区间找到极大值点，即为O波，令O波幅值≤C波幅值的1/3。最后重新标记A点和O点为周期样本的起始点和终止点。Specifically, the samples with large amplitudes of A-wave and O-wave due to noise are excluded, the commonality with the standard samples is retained, and the accuracy of feature point extraction is improved. By traversing and comparing the maximum point by point, the maximum value point before the Q point is found, which is the A wave, and the amplitude of the A wave is ≤ 1/3 of the amplitude of the C wave. Similarly, traverse the interval after the boundary of point X to find the maximum value point, which is the O wave, so that the amplitude of the O wave is ≤ 1/3 of the amplitude of the C wave. Finally, point A and O are re-marked as the start and end points of the periodic sample.
S105、对阻抗微分图的特征点B点进行约束，限制B点位置在QC间期，剔除无效样本。S105 , constrain point B of the characteristic point of the impedance differential graph, limit the position of point B in the QC interval, and eliminate invalid samples.
具体的，在以C点为界限的前段区间，通过逐点遍历找到C波幅值的15％的点，将其标记为阻抗微分图的特征点B点，判断B点位置是否位于经过数据分割后的QC间期，剔除B点偏离QC间期的无效样本周期。Specifically, in the front section bounded by point C, find the point with 15% of the amplitude of C wave through point-by-point traversal, mark it as point B of the characteristic point of the impedance differential map, and determine whether the position of point B is located after the data segmentation. After the QC interval, the invalid sample cycle that deviates from the QC interval at point B is eliminated.
S106、判断阻抗微分图中QC间期出现拐点的情况，对特征点B，C，X点进行标记。S106 , judging that an inflection point occurs in the QC interval in the impedance differential diagram, and marking the feature points B, C, and X.
具体的，参阅图6，遍历阻抗微分图周期样本的最大值和最小值，分别为特征点C点和X点。判断QC间期中是否出现拐点(d2Z/dt2＝0的点)，即diff[i]<0,diff[i+1]>0，若不存在拐点，则在保证B点位于QC间期的条件下，标记C波幅值的15％的位置为特征点B点；若存在拐点，则判断QC间期中出现拐点的数量，若大于1个则对这样的无效样本进行剔除，若仅出现一个拐点，则该点即为特征点B点。Specifically, referring to FIG. 6 , the maximum and minimum values of the periodic samples of the traversal impedance differential graph are the feature points C and X, respectively. Determine whether there is an inflection point in the QC interval (d2Z/dt2=0 point), that is, diff[i]<0, diff[i+1]>0, if there is no inflection point, then ensure that point B is located in the QC interval. If there is an inflection point, the number of inflection points in the QC interval is judged. If more than 1, such invalid samples will be eliminated. If there is only one inflection point , then this point is the feature point B.
S107、利用特征点形成阻抗血流特征参数，并建立血压模型。S107 , using the feature points to form impedance blood flow characteristic parameters, and establish a blood pressure model.
具体的，利用特征点B，C，X点构建与血压相关的阻抗血流特征参数包括：C波幅值HC，B点幅值HB，X点幅值HX，射血时间TBX，BC上升期时间TBC，CX下降期时间TCX，射血时间占周期比RBX，BC上升期时间占周期比RBC，CX下降期时间占周期比RCX，BC段斜率KBC，CX段斜率KCX。Specifically, using the characteristic points B, C, and X to construct the blood pressure-related impedance blood flow characteristic parameters, including: C wave amplitude HC , B point amplitude HB , X point amplitude HX , and ejection time TBX , BC rising phase time TBC , CX falling phase time TCX , ejection time to cycle ratio RBX , BC rising phase time to cycle ratio RBC , CX fall time to cycle ratio RCX , BC segment slope KBC , CX segment slope KCX .
然后通过多元回归的方法，根据特征参数建立血压多元线性模型，以测量的血压数值为因变量，以阻抗血流参数为自变量构建多元线性回归方程，最后得到舒张压DBP的计算表达式为：Then, through the method of multiple regression, a multiple linear model of blood pressure is established according to the characteristic parameters, the measured blood pressure value is used as the dependent variable, and the multiple linear regression equation is constructed with the impedance blood flow parameter as the independent variable. Finally, the calculation expression of the diastolic blood pressure DBP is obtained as:
DBP＝a1+b1*HB+c1*HX+d1*TBX+e1*TCX+f1*RCX+g1*RBX+h1*KCXDBP=a1 +b1 *HB +c1 *HX +d1 *TBX +e1 *TCX +f1 *RCX +g1 *RBX +h1 *KCX
收缩压SBP的计算表达式为：The calculation expression of systolic blood pressure SBP is:
其通过提取阻抗血流图的阻抗血流参数和使用血压计同步测量收缩压和舒张压，将测量者的阻抗血流参数数据组和对应的收缩压、舒张压代入血压模型训练，便可得到与之匹配的模型参数值。It extracts the impedance blood flow parameters of the impedance blood flow map and uses the sphygmomanometer to measure systolic blood pressure and diastolic blood pressure simultaneously, and substitutes the measurer's impedance blood flow parameter data set and the corresponding systolic blood pressure and diastolic blood pressure into the blood pressure model training. The model parameter value to match with.
本发明的一种用于血压检测的阻抗血流图处理方法，通过确定阻抗微分图的起始点和终止点，对阻抗微分图的周期样本区间进行相位平移，数据分割和区间搜索来实现A波和O波进行幅值约束，剔除无效样本，并且重新标记起始点和终止点，特征点B点进行约束，限制B点位置在QC间期，判断阻抗微分图中QC间期出现拐点的情况，对特征点B，C，X点进行标记，利用特征点形成阻抗血流特征参数，并建立血压模型。最终提高了特征点提取的准确性，通过血流参数建立的血压模型与用电子血压计标定的相比，收缩压的绝对误差为3.85mmHg，相对误差为3.81％；舒张压的绝对误差为2.85mmHg，相对误差为3.85％。The present invention provides an impedance blood flow diagram processing method for blood pressure detection, which realizes the A wave by determining the starting point and the ending point of the impedance differential diagram, performing phase shift, data division and interval search on the periodic sample interval of the impedance differential diagram. Constrain the amplitude with the O wave, remove the invalid samples, and re-mark the start and end points. The feature point B is constrained to limit the position of point B to the QC interval, and determine the inflection point in the QC interval in the impedance differential graph. Mark the characteristic points B, C and X, and use the characteristic points to form the characteristic parameters of impedance blood flow, and establish the blood pressure model. Finally, the accuracy of feature point extraction is improved. Compared with the blood pressure model established by blood flow parameters, the absolute error of systolic blood pressure is 3.85mmHg, and the relative error is 3.81%; the absolute error of diastolic blood pressure is 2.85. mmHg, the relative error is 3.85%.
Claims (8)
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210244270.1ACN114601440A (en) | 2022-03-14 | 2022-03-14 | Impedance rheogram processing method for blood pressure detection |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210244270.1ACN114601440A (en) | 2022-03-14 | 2022-03-14 | Impedance rheogram processing method for blood pressure detection |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN114601440Atrue CN114601440A (en) | 2022-06-10 |
Family
ID=81863357
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202210244270.1APendingCN114601440A (en) | 2022-03-14 | 2022-03-14 | Impedance rheogram processing method for blood pressure detection |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114601440A (en) |
Citations (15)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101332084A (en)* | 2007-06-28 | 2008-12-31 | 罗晓民 | Area barycenter displacement analytical method of pulse wave |
| CN101791216A (en)* | 2008-12-26 | 2010-08-04 | 罗晓民 | System for detecting blood flow parameters based on pulse wave measurement and analysis |
| CN101925331A (en)* | 2007-11-28 | 2010-12-22 | 爱德华兹生命科学Iprm公司 | Method and apparatus to determine the end of systolic part of pressure curve |
| CN103027667A (en)* | 2011-09-30 | 2013-04-10 | Ge医疗系统环球技术有限公司 | Characteristic parameter extraction of pulse wave |
| CN103169456A (en)* | 2013-03-29 | 2013-06-26 | 深圳职业技术学院 | Processing method and processing system for pulse wave signals |
| CN104173030A (en)* | 2014-09-09 | 2014-12-03 | 北京航空航天大学 | Pulse wave starting point real-time detection method resisting waveform change interference and application thereof |
| CN105433916A (en)* | 2014-09-19 | 2016-03-30 | 株式会社日立制作所 | Pulse feature extraction device and method |
| CN106646029A (en)* | 2016-11-18 | 2017-05-10 | 武汉电联电力电气技术有限公司 | Capacitor measurement device and reactor characteristic testing method |
| CN109171677A (en)* | 2018-09-20 | 2019-01-11 | 芯海科技（深圳）股份有限公司 | A kind of Portable type measurement unit and the hemodynamic parameter measurement method based on the device |
| CN109691993A (en)* | 2018-12-07 | 2019-04-30 | 芯海科技（深圳）股份有限公司 | A kind of method and human body balance measuring heart rate variability |
| CN110495863A (en)* | 2018-10-08 | 2019-11-26 | 合云力健康科技(北京)有限公司 | The method and apparatus for identifying radial artery wave shape dicrotic notch characteristic point |
| CN111202510A (en)* | 2020-01-21 | 2020-05-29 | 桂林电子科技大学 | Data processing method of impedance hemogram |
| CN111248893A (en)* | 2020-01-21 | 2020-06-09 | 桂林电子科技大学 | A method for establishing blood flow map based on bioelectrical impedance technology |
| CN111938640A (en)* | 2020-08-05 | 2020-11-17 | 深圳大学 | Method, device and storage medium for C-point positioning of cardiac impedance differential signal |
| CN113128350A (en)* | 2021-03-25 | 2021-07-16 | 同济大学 | Time domain identification and positioning method, equipment and medium for PPG signal feature points |
- 2022
- 2022-03-14CNCN202210244270.1Apatent/CN114601440A/enactivePending
Patent Citations (15)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101332084A (en)* | 2007-06-28 | 2008-12-31 | 罗晓民 | Area barycenter displacement analytical method of pulse wave |
| CN101925331A (en)* | 2007-11-28 | 2010-12-22 | 爱德华兹生命科学Iprm公司 | Method and apparatus to determine the end of systolic part of pressure curve |
| CN101791216A (en)* | 2008-12-26 | 2010-08-04 | 罗晓民 | System for detecting blood flow parameters based on pulse wave measurement and analysis |
| CN103027667A (en)* | 2011-09-30 | 2013-04-10 | Ge医疗系统环球技术有限公司 | Characteristic parameter extraction of pulse wave |
| CN103169456A (en)* | 2013-03-29 | 2013-06-26 | 深圳职业技术学院 | Processing method and processing system for pulse wave signals |
| CN104173030A (en)* | 2014-09-09 | 2014-12-03 | 北京航空航天大学 | Pulse wave starting point real-time detection method resisting waveform change interference and application thereof |
| CN105433916A (en)* | 2014-09-19 | 2016-03-30 | 株式会社日立制作所 | Pulse feature extraction device and method |
| CN106646029A (en)* | 2016-11-18 | 2017-05-10 | 武汉电联电力电气技术有限公司 | Capacitor measurement device and reactor characteristic testing method |
| CN109171677A (en)* | 2018-09-20 | 2019-01-11 | 芯海科技（深圳）股份有限公司 | A kind of Portable type measurement unit and the hemodynamic parameter measurement method based on the device |
| CN110495863A (en)* | 2018-10-08 | 2019-11-26 | 合云力健康科技(北京)有限公司 | The method and apparatus for identifying radial artery wave shape dicrotic notch characteristic point |
| CN109691993A (en)* | 2018-12-07 | 2019-04-30 | 芯海科技（深圳）股份有限公司 | A kind of method and human body balance measuring heart rate variability |
| CN111202510A (en)* | 2020-01-21 | 2020-05-29 | 桂林电子科技大学 | Data processing method of impedance hemogram |
| CN111248893A (en)* | 2020-01-21 | 2020-06-09 | 桂林电子科技大学 | A method for establishing blood flow map based on bioelectrical impedance technology |
| CN111938640A (en)* | 2020-08-05 | 2020-11-17 | 深圳大学 | Method, device and storage medium for C-point positioning of cardiac impedance differential signal |
| CN113128350A (en)* | 2021-03-25 | 2021-07-16 | 同济大学 | Time domain identification and positioning method, equipment and medium for PPG signal feature points |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106974631B (en) | A kind of blood pressure measurement method and device based on pulse wave waveform and ECG signal | |
| CN111839488B (en) | Device and method for non-invasive continuous blood pressure measurement based on pulse wave | |
| CN104116503B (en) | A kind of measurement apparatus of noinvasive continuous blood pressure | |
| CN102270264B (en) | Physiological signal quality assessment system and method | |
| CN114145724A (en) | Method for dynamically monitoring blood pressure based on ECG (electrocardiogram) and PPG (photoplethysmography) multiple physiological characteristic parameters | |
| CN112089405B (en) | Pulse wave characteristic parameter measuring and displaying device | |
| CN109124610B (en) | Anti-interference method and device for non-invasive blood pressure measurement | |
| KR20120108575A (en) | Apparatus and method improving accuracy of wrist blood pressure by using multiple bio-signal | |
| CN106413534A (en) | Continuous blood pressure measurement device, measurement model establishment method and system | |
| CN112545472B (en) | PPG signal quality evaluation method, device, equipment and storage medium | |
| CN110840428B (en) | Noninvasive blood pressure estimation method based on one-dimensional U-Net network | |
| CN104382571A (en) | Method and device for measuring blood pressure upon radial artery pulse wave conduction time | |
| CN111000544B (en) | Construction method and system of hybrid continuous blood pressure measurement model based on PPG waveform | |
| CN113598737A (en) | Blood pressure prediction method based on feature fusion | |
| CN102885616A (en) | Method for removing baseline drift from pulse wave signal | |
| CN108354597A (en) | A kind of rapid blood pressure computational methods based on the extraction of optimal wave | |
| CN113576438A (en) | Non-invasive blood pressure extraction method and system | |
| WO2023226293A1 (en) | Doppler ultrasound-based continuous blood pressure measuring apparatus and electronic device | |
| Lin et al. | Investigation on pulse wave forward peak detection and its applications in cardiovascular health | |
| CN102755154B (en) | Calculation method for extracting conduction time from pulse wave | |
| CN113243902B (en) | Feature extraction method based on photoplethysmography | |
| CN208677373U (en) | A kind of pulse analysis system | |
| CN117694863B (en) | Hemodynamic parameter evaluation device and medium | |
| CN114601440A (en) | Impedance rheogram processing method for blood pressure detection | |
| CN116570259A (en) | Non-invasive blood pressure measurement method based on pulse wave characteristics |
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 |