CN103284702A - Electrocardiogram and pulse wave relation analysis method and method and device of fusion analysis - Google Patents
Electrocardiogram and pulse wave relation analysis method and method and device of fusion analysis Download PDFInfo
- Publication number
- CN103284702A CN103284702A CN2013102395557A CN201310239555A CN103284702A CN 103284702 A CN103284702 A CN 103284702A CN 2013102395557 A CN2013102395557 A CN 2013102395557A CN 201310239555 A CN201310239555 A CN 201310239555A CN 103284702 A CN103284702 A CN 103284702A
- Authority
- CN
- China
- Prior art keywords
- pulse wave
- electrocardiogram
- wave
- ripple
- cycle
- 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
- 238000004458 analytical method Methods 0.000 title claims abstract description 53
- 230000004927 fusion Effects 0.000 title claims abstract description 43
- 238000000034 method Methods 0.000 title claims abstract description 41
- 238000001514 detection method Methods 0.000 claims description 41
- 238000000605 extraction Methods 0.000 claims description 39
- 239000000284 extract Substances 0.000 claims description 27
- 230000002123 temporal effect Effects 0.000 claims description 20
- 201000010099 disease Diseases 0.000 claims description 17
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 claims description 17
- 238000007781 pre-processing Methods 0.000 claims description 11
- 238000013135 deep learning Methods 0.000 claims description 9
- 230000000007 visual effect Effects 0.000 claims description 8
- 238000013519 translation Methods 0.000 claims description 6
- 238000003786 synthesis reaction Methods 0.000 claims 3
- 230000007935 neutral effect Effects 0.000 claims 2
- 238000005259 measurement Methods 0.000 abstract description 2
- 238000003745 diagnosis Methods 0.000 abstract 1
- 238000007689 inspection Methods 0.000 abstract 1
- 206010003658 Atrial Fibrillation Diseases 0.000 description 10
- 238000010835 comparative analysis Methods 0.000 description 6
- 230000000877 morphologic effect Effects 0.000 description 6
- 238000004364 calculation method Methods 0.000 description 5
- 238000001914 filtration Methods 0.000 description 5
- 238000013528 artificial neural network Methods 0.000 description 4
- 208000024172 Cardiovascular disease Diseases 0.000 description 3
- 230000006870 function Effects 0.000 description 3
- 230000000737 periodic effect Effects 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 206010003662 Atrial flutter Diseases 0.000 description 2
- 238000004422 calculation algorithm Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000010606 normalization Methods 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- 238000012800 visualization Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000017531 blood circulation Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000007635 classification algorithm Methods 0.000 description 1
- 238000003759 clinical diagnosis Methods 0.000 description 1
- 230000000295 complement effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000007717 exclusion Effects 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 230000001902 propagating effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
Images
Landscapes
- Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
- Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)
Abstract
Description
【技术领域】【Technical field】
本发明涉及用于诊断目的的测量技术领域,特别是涉及一种心电图和脉搏波关系分析及在此基础上融合分析的方法和装置。The invention relates to the technical field of measurement for diagnostic purposes, in particular to a method and device for analyzing the relationship between an electrocardiogram and a pulse wave and based on the fusion analysis.
【背景技术】【Background technique】
心电图是指把心脏周期活动产生的电信号记录下来并按规定的格式绘制出来的图形。脉搏波是心脏的搏动沿动脉血管和血流向外周传播而形成的。二者产生的形式虽然不同,但二者的来源都是心脏,所以研究二者之间互斥、一致和互补的关系有着重要的意义。An electrocardiogram is a graph that records the electrical signals generated by the heart cycle and draws them in a prescribed format. The pulse wave is formed by the pulse of the heart propagating along the arterial vessels and blood flow to the periphery. Although the forms produced by the two are different, the source of both is the heart, so it is of great significance to study the mutually exclusive, consistent and complementary relationship between the two.
心电图已经被广泛应用于心血管疾病的临床诊断中。近些年,为了减少远程监护、体检中心等医生的工作量,心电图自动分析的需求越来越迫切。然而,现有的心电图分类方法在实际应用中的分类准确率还是不尽如人意。除了心电图分类算法本身的问题之外,单一的心电图信号所包含的信息量也有限。所以,融合多种信号的诊断分析越来越被重视,但融合的前提是要了解信号之间的关系,才能给出更准确地分析结论。ECG has been widely used in the clinical diagnosis of cardiovascular diseases. In recent years, in order to reduce the workload of doctors in remote monitoring and medical examination centers, the demand for automatic ECG analysis has become more and more urgent. However, the classification accuracy of existing ECG classification methods in practical applications is still not satisfactory. In addition to the problems of the ECG classification algorithm itself, the amount of information contained in a single ECG signal is also limited. Therefore, more and more attention is paid to the diagnostic analysis of fusion of multiple signals, but the premise of fusion is to understand the relationship between signals in order to give more accurate analysis conclusions.
现有技术中,尚未有对心电图和脉搏波关系进行分析,以及在此分析基础上对采样数据进行融合分析的方法和装置。In the prior art, there is no method and device for analyzing the relationship between the electrocardiogram and the pulse wave, and performing fusion analysis on the sampled data based on the analysis.
鉴于此,克服该现有技术所存在的缺陷是本技术领域亟待解决的问题。In view of this, it is an urgent problem to be solved in this technical field to overcome the defects in the prior art.
【发明内容】【Content of invention】
本发明要解决的技术问题是提供一种融合多种信号的心电图和脉搏波关系分析方法和装置。The technical problem to be solved by the present invention is to provide a method and device for analyzing the relationship between an electrocardiogram and a pulse wave that fuses multiple signals.
本发明进一步要解决的技术问题是提供一种在前述关系分析基础上的心电图和脉搏波融合分析方法和装置。The further technical problem to be solved by the present invention is to provide a method and device for fusion analysis of electrocardiogram and pulse wave based on the aforementioned relationship analysis.
本发明采用如下技术方案:The present invention adopts following technical scheme:
一种心电图和脉搏波关系分析方法,包括:A method for analyzing the relationship between electrocardiogram and pulse wave, comprising:
步骤S1:对获取的一已知疾病对象的心电图和脉搏波信号分别进行预处理,包括基线平移和去噪;Step S1: Perform preprocessing on the acquired ECG and pulse wave signals of a subject with a known disease, including baseline translation and denoising;
步骤S2:对预处理后的心电图和脉搏波进行二维和/或三维图像的可视化呈现;Step S2: performing visual presentation of two-dimensional and/or three-dimensional images on the preprocessed electrocardiogram and pulse wave;
步骤S3:提取心电图特征和脉搏波特征;Step S3: extracting ECG features and pulse wave features;
步骤S4:对心电图特征和脉搏波特征进行对比分析,得出心电图和脉搏波信号各特征之间的关系。Step S4: Perform comparative analysis on the features of the electrocardiogram and pulse wave signals to obtain the relationship between the features of the electrocardiogram and pulse wave signals.
本发明还提供了一种心电图和脉搏波关系分析装置,包括:The present invention also provides a device for analyzing the relationship between electrocardiogram and pulse wave, comprising:
预处理模块,用于对获取的心电图和脉搏波信号分别进行预处理,包括基线平移和去噪;The preprocessing module is used to preprocess the acquired ECG and pulse wave signals, including baseline translation and denoising;
可视化呈现模块,用于对预处理后的心电图和脉搏波进行二维和/或三维图像的可视化呈现;A visualization presentation module, used for visual presentation of two-dimensional and/or three-dimensional images of the preprocessed electrocardiogram and pulse wave;
心电图特征提取模块,用于提取心电图特征;The electrocardiogram feature extraction module is used to extract the electrocardiogram feature;
脉搏波特征提取模块,用于提取脉搏波特征;The pulse wave feature extraction module is used to extract the pulse wave feature;
对比分析模块,用于对心电图特征和脉搏波特征进行对比分析,得出心电图和脉搏波信号各特征之间的关系。The comparative analysis module is used for comparatively analyzing the characteristics of the electrocardiogram and the pulse wave, and obtaining the relationship between the characteristics of the electrocardiogram and the pulse wave signal.
本发明还提供了一种心电图和脉搏波融合分析方法,包括:The present invention also provides a kind of electrocardiogram and pulse wave fusion analysis method, comprising:
步骤Q1:对获取的一未知疾病对象的心电图和脉搏波信号进行决策层融合;Step Q1: Decision-making fusion of the acquired ECG and pulse wave signals of an unknown disease subject;
步骤Q2:采用深度学习分类器,利用深层神经网络对所述未知疾病对象的心电图和脉搏波信号进行数据层和特征层融合;Step Q2: Using a deep learning classifier, using a deep neural network to perform data layer and feature layer fusion on the ECG and pulse wave signals of the unknown disease subject;
步骤Q3:对所述决策层融合、数据层和特征层融合后的心电图和脉搏波信号各特征之间的关系进行综合分析,并根据上述关系分析方法得到的心电图和脉搏波信号各特征之间的关系进行识别分类。Step Q3: Comprehensively analyze the relationship between the features of the electrocardiogram and the pulse wave signal after the fusion of the decision-making layer, the data layer and the feature layer, and obtain the relationship between the features of the electrocardiogram and the pulse wave signal according to the above relationship analysis method relationship to identify and classify.
本发明还提供了一种心电图和脉搏波融合分析装置,包括:The present invention also provides an electrocardiogram and pulse wave fusion analysis device, comprising:
决策层融合模块,用于对获取的一未知疾病对象的心电图和脉搏波信号进行决策层融合;The decision-making layer fusion module is used to carry out decision-making layer fusion to the electrocardiogram and pulse wave signal of an unknown disease object acquired;
数据层与特征层融合模块,包括一深度学习分类器,所述深度学习分类器利用深层神经网络对所述未知疾病对象的心电图和脉搏波信号进行数据层和特征层融合;The data layer and feature layer fusion module includes a deep learning classifier, and the deep learning classifier utilizes a deep neural network to perform data layer and feature layer fusion on the electrocardiogram and pulse wave signal of the unknown disease object;
综合分析模块,用于对所述决策层融合、数据层和特征层融合后的心电图和脉搏波信号各特征之间的关系进行综合分析,并根据上述关系分析装置得到的心电图和脉搏波信号各特征之间的关系进行识别分类。The comprehensive analysis module is used to comprehensively analyze the relationship between the characteristics of the electrocardiogram and the pulse wave signal after the fusion of the decision-making layer, the data layer and the feature layer, and analyze the relationship between the electrocardiogram and the pulse wave signal obtained by the above-mentioned relationship analysis device. The relationship between features is identified and classified.
与现有技术相比,本发明的有益效果在于:本发明对心电图和脉搏波的信号特征进行关系分析,得出其各特征之间的相互关系,并依据相互的关系进行融合分析检测,与单一信号的分析相比,提高了结果的分类准确率。Compared with the prior art, the beneficial effect of the present invention is that: the present invention carries out relational analysis to the signal characteristic of electrocardiogram and pulse wave, obtains the interrelationship between its each characteristic, and carries out fusion analysis detection according to mutual relation, and Compared with the analysis of a single signal, the classification accuracy of the results is improved.
【附图说明】【Description of drawings】
图1是本发明实施例1一种心电图和脉搏波关系分析方法流程图;Fig. 1 is a flow chart of a method for analyzing the relationship between an electrocardiogram and a pulse wave in Embodiment 1 of the present invention;
图2是本发明实施例2一种心电图和脉搏波关系分析装置结构框图;Fig. 2 is a structural block diagram of a device for analyzing the relationship between electrocardiogram and pulse wave in Embodiment 2 of the present invention;
图3是本发明实施例3一种心电图和脉搏波融合分析方法流程图;Fig. 3 is a flow chart of an electrocardiogram and pulse wave fusion analysis method according to
图4是本发明实施例4一种心电图和脉搏波融合分析装置结构框图。Fig. 4 is a structural block diagram of an electrocardiogram and pulse wave fusion analysis device according to
【具体实施方式】【Detailed ways】
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。In order to make the object, technical solution and advantages of the present invention clearer, the present invention will be further described in detail below in conjunction with the accompanying drawings and embodiments. It should be understood that the specific embodiments described here are only used to explain the present invention, not to limit the present invention.
此外,下面所描述的本发明各个实施方式中所涉及到的技术特征只要彼此之间未构成冲突就可以相互组合。In addition, the technical features involved in the various embodiments of the present invention described below can be combined with each other as long as they do not constitute a conflict with each other.
本发明提出的心电图和脉搏波关系分析方法,得到的心电图和脉搏波信号各特征之间的相互关系包括互斥性、一致性和互补性,According to the electrocardiogram and pulse wave relationship analysis method proposed by the present invention, the relationship between the obtained electrocardiogram and pulse wave signal features includes mutual exclusion, consistency and complementarity,
本发明提出的心电图和脉搏波融合分析的方法,是基于本发明的分析心电图和脉搏波信号之间关系的方法,在确定已知疾病对象的心电图和脉搏波确切关系的基础上进行心电图和脉搏波融合分析,与单一信号的分析相比,提高了结果的分类准确率。The method for electrocardiogram and pulse wave fusion analysis proposed by the present invention is based on the method for analyzing the relationship between electrocardiogram and pulse wave signals of the present invention, and performs electrocardiogram and pulse wave analysis on the basis of determining the exact relationship between electrocardiogram and pulse wave signals of known disease objects. Wave fusion analysis, compared with the analysis of a single signal, improves the classification accuracy of the results.
实施例1Example 1
本发明实施例1提供了一种心电图和脉搏波关系分析方法。本实施例中的心电图和脉搏波数据为针对某一具体已知的心血管疾病获取的心电图和脉搏波数据,例如但不限于房颤疾病,此处不一一列举。以房颤疾病为例,首先需获得房颤疾病的脉搏波和心电图数据。Embodiment 1 of the present invention provides a method for analyzing the relationship between an electrocardiogram and a pulse wave. The electrocardiogram and pulse wave data in this embodiment are the electrocardiogram and pulse wave data obtained for a specific known cardiovascular disease, such as but not limited to atrial fibrillation, which are not listed here. Taking atrial fibrillation as an example, it is first necessary to obtain the pulse wave and electrocardiogram data of atrial fibrillation.
如图1所示,该方法包括如下步骤:As shown in Figure 1, the method includes the following steps:
步骤S1:对获取的一已知疾病对象的心电图和脉搏波信号分别进行预处理,包括基线平移和去噪;Step S1: Perform preprocessing on the acquired ECG and pulse wave signals of a subject with a known disease, including baseline translation and denoising;
步骤S2:对预处理后的心电图和脉搏波进行二维和/或三维图像的可视化呈现;Step S2: performing visual presentation of two-dimensional and/or three-dimensional images on the preprocessed electrocardiogram and pulse wave;
步骤S3:提取心电图特征和脉搏波特征;Step S3: extracting ECG features and pulse wave features;
步骤S4:对心电图特征和脉搏波特征进行对比分析,得出心电图和脉搏波信号各特征之间的关系。例如:在房颤疾病数据下,脉率与心率不一致等。Step S4: Perform comparative analysis on the features of the electrocardiogram and pulse wave signals to obtain the relationship between the features of the electrocardiogram and pulse wave signals. For example: in the data of atrial fibrillation, the pulse rate is inconsistent with the heart rate, etc.
下面对上述步骤进行详细说明:The above steps are described in detail below:
步骤S1预处理时,以房颤疾病为例,可选择1-30Hz的带通滤波对原始脉搏波x[n]进行滤波。带通滤波器采用Matlab设计,采用6阶、两个通带截止点分别为1Hz和30Hz的切比雪夫滤波器(chebyshev)I型无限响应滤波器。In step S1 preprocessing, taking atrial fibrillation as an example, a band-pass filter of 1-30 Hz can be selected to filter the original pulse wave x[n]. The band-pass filter is designed by Matlab, using a 6th-order Chebyshev filter (chebyshev) type I infinite response filter with two pass-band cut-off points of 1Hz and 30Hz respectively.
经过步骤S1预处理后的心电图和脉搏波,在步骤S2中进行可视化呈现,通过二维或三维图像显示,以此确定采集位置的准确程度,即是否获得有效波形信号。此外,还可以找出直观特征之间心电图与脉搏波二者之间的关系。以房颤疾病为例,获得12点阵脉搏波波形,以确定采集位置的准确性,根据脉搏波在脉管中央波动较大的特点,由振幅波动强弱判断检测位置,从而保证信号质量。The electrocardiogram and pulse wave preprocessed in step S1 are visualized in step S2 and displayed in two-dimensional or three-dimensional images to determine the accuracy of the acquisition position, that is, whether to obtain valid waveform signals. In addition, the relationship between the electrocardiogram and the pulse wave among the intuitive features can also be found out. Taking atrial fibrillation as an example, the 12-dot matrix pulse wave waveform is obtained to determine the accuracy of the acquisition position. According to the characteristics of the pulse wave fluctuating greatly in the center of the vessel, the detection position is judged by the strength of the amplitude fluctuation, so as to ensure the signal quality.
步骤S3可通过如下方法提取心电图特征和脉搏波特征:Step S3 can extract electrocardiogram feature and pulse wave feature by following method:
步骤S3a:进行心电图R波检测,提取心电图中R波位置;Step S3a: performing electrocardiogram R-wave detection, and extracting the R-wave position in the electrocardiogram;
步骤S3b:进行心电图中除R波以外的其他时域特征检测,提取心电图中除R波以外的其他时域特征,包括但不限于QRS波等的提取;Step S3b: Perform other time-domain feature detection in the electrocardiogram except the R wave, and extract other time-domain features in the electrocardiogram except the R wave, including but not limited to the extraction of QRS waves;
步骤S3c:进行脉搏波主波检测,提取脉搏波中主波位置;Step S3c: detecting the main wave of the pulse wave, and extracting the position of the main wave in the pulse wave;
步骤S3d:进行脉搏波时域特征检测,包括:以主波位置为依据,进行脉搏波周期起始点检测,确定脉搏波周期,检测脉搏波周期是否为噪声周期,对非噪声周期检测重博前波、重博波位置幅值和主波及周期形态。Step S3d: Perform pulse wave time-domain feature detection, including: based on the position of the main wave, detect the starting point of the pulse wave cycle, determine the pulse wave cycle, detect whether the pulse wave cycle is a noise cycle, and detect whether the pulse wave cycle is a non-noise cycle. wave, double wave position amplitude and main wave and cycle form.
检测脉搏波周期是否为噪声周期可通过如下方法:获取平均脉搏波周期;以平均脉搏波周期为基础,计算每一个经过插值后的脉搏波周期的余弦形相似度,与平均脉搏波周期的相似度值小于一设定值S的脉搏波周期为噪声周期,并直接排除在外。该设定值S可根据实际情况设定,例如取S=0.8。The following method can be used to detect whether the pulse wave cycle is a noise cycle: obtain the average pulse wave cycle; based on the average pulse wave cycle, calculate the cosine similarity of each interpolated pulse wave cycle, which is similar to the average pulse wave cycle The pulse wave period whose degree value is less than a set value S is a noise period and is directly excluded. The setting value S can be set according to actual conditions, for example, S=0.8.
其中,对非噪声周期检测重博前波、重博波位置幅值和主波及周期形态可根据医生经验规则进行,能提高时域特征提取的准确率,尤其是重博前波。根据中医生经验和中医脉象研究的结论,目前脉搏波结束点前的大约0.05秒意义尚不明确,另外,如果从主波位置开始检测,可能会将主波误检为其他特征点,所以此处从主波后t秒(例如取t=0.03)开始检测其他时域特征点。最终确定重博波和重博前波检测区间t1,例如可设t1=[p[m]+0.03:pend[m]-0.05]为时域特征检测范围,其中p[m]表示第m个主波位置,pend[m]表示对应的脉搏波周期结束点位置。根据统计发现,重搏前波若存在则出现范围大都位于t1范围的前1/3内,所以取重博前波检测范围t2,具体识别方法则是根据一次、二次差分信号找到波峰及波谷。该方法不涉及复杂的步骤,在相同的数据集上仅需提取时域特征点,将提取的结果与标注结果进行对比,准确率高。主波、重搏前波、重搏波三个明显时域特征准确率分别为99.71%、88.68%、89.03%,高于其他方法。Among them, the non-noise periodic detection of the dichroic front wave, the position amplitude of the dichroic wave, the main wave and the periodic shape can be carried out according to the doctor's empirical rules, which can improve the accuracy of time-domain feature extraction, especially the dichroic front wave. According to the experience of TCM doctors and the conclusions of TCM pulse research, the meaning of about 0.05 seconds before the end of the pulse wave is not clear. In addition, if the main wave is detected from the position of the main wave, the main wave may be misdetected as other feature points, so this Detect other time-domain feature points from t seconds after the main wave (for example, t=0.03). Finally determine the detection interval t1 of the double wave and the heavy wave front wave. For example, t1=[p[m]+0.03:pend[m]-0.05] can be set as the time domain feature detection range, where p[m] represents the mth Main wave position, pend[m] indicates the position of the end point of the corresponding pulse wave cycle. According to the statistics, it is found that if the dicrotic front wave exists, the range of occurrence is mostly within the first 1/3 of the range of t1, so the detection range of the dicrotic front wave is t2, and the specific identification method is to find the peak and trough based on the primary and secondary differential signals . This method does not involve complicated steps, and only needs to extract time-domain feature points on the same data set, and compare the extracted results with the labeled results, with high accuracy. The accuracy rates of three obvious time-domain features of main wave, dicrotic wave and dicrotic wave are 99.71%, 88.68%, and 89.03%, respectively, which are higher than other methods.
具体地,可采用如下方法进行心电图R波检测:Specifically, the following methods can be used to detect the R wave of the electrocardiogram:
1a)带通滤波:使用5-18Hz的前向滤波器对对输入的心电图信号进行滤波并做相位延迟补偿;1a) Band-pass filtering: use a 5-18Hz forward filter to filter the input ECG signal and perform phase delay compensation;
1b)差分:对前向滤波器输出的信号进行差分处理,形成差分信号;1b) Difference: Perform differential processing on the signal output by the forward filter to form a differential signal;
1c)数据整理:对差分信号进行变换:输出=输入的绝对值/G1,G1为假定的人类R波差分值最大值,可取值0.1~0.5,若输出等于1则设置成1,若输出小于0.01也置1;该数据整理过程不使用归一化方式,而是用一个绝对值,这样在后面检测中可以计算信号幅度,从而检测出停搏;1c) Data sorting: transform the differential signal: output = absolute value of input/G1, G1 is the assumed maximum value of human R-wave differential value, and can take a value of 0.1 to 0.5. If the output is equal to 1, set it to 1. If the output It is also set to 1 if it is less than 0.01; the data sorting process does not use the normalization method, but uses an absolute value, so that the signal amplitude can be calculated in the subsequent detection, so as to detect the arrest;
2a)对经数据整理后的信号使用公式d(n)*d(n)*log(d(n)*d(n))进行香农能量转换;2a) Use the formula d(n)*d(n)*log(d(n)*d(n)) to perform Shannon energy conversion on the signal after data arrangement;
2b)平均滤波:使用M=55~75点(153~208ms)的前向滤波器进行滤波并做相位延迟补偿;2b) Average filtering: Use a forward filter with M=55~75 points (153~208ms) for filtering and phase delay compensation;
3a)检测极大/小点:极大点是指值大于左边又大于右边的点,极小点是指值小于左边又小于右边的点;3a) Detection of maximum/small points: a maximum point refers to a point whose value is greater than the left and greater than the right, and a minimum point refers to a point whose value is smaller than the left and smaller than the right;
3b)排除假R点:假设n点是极大点,它只需满足以下任一个条件就将被排除:3b) Exclude false R points: Assuming that point n is a maximum point, it will be excluded as long as it meets any of the following conditions:
(a)如果在n点停搏时限范围内差分信号小于G2a时,则此处应为停搏,n点对应的极大点不是R波,G2a为R波差分后最小允许有效值,可取值为0.01~0.06,其取值跟信号噪声有关;(a) If the differential signal is less than G2a within the pause time limit of point n, it should be pause here, the maximum point corresponding to point n is not R wave, and G2a is the minimum allowable effective value after R wave difference, which is desirable The value is 0.01~0.06, and its value is related to the signal noise;
(b)如果n点在平均滤波后的数据值小于G2b,则排除该点,G2b为允许R波最小峰值,可取0.005~0.01;(b) If the data value of point n after average filtering is less than G2b, then exclude this point. G2b is the minimum peak value of the allowed R wave, which can be 0.005-0.01;
(c)如果n点与旁边极小点的差值比n点停搏时限范围内的最大值的G3倍小则被排除,G3为允许R波幅度突然变小的最大比例,可以取0.06~0.2;(c) If the difference between point n and the minimum point next to it is smaller than G3 times the maximum value within the arrest time limit of point n, it will be excluded. G3 is the maximum ratio that allows the amplitude of R wave to suddenly decrease, and it can be 0.06~ 0.2;
3c)纠正误排除点:该操作可以纠正突然R波变小的情况,进一步提高准确率,其过程为:假设n点是被排除的极大点,如果它全部满足以下条件则认为是误排除的R波:3c) Correct mis-excluded points: This operation can correct the situation where the R wave becomes smaller suddenly, and further improve the accuracy rate. The process is as follows: Assume that n points are the maximum points to be excluded, and if all of them meet the following conditions, it is considered to be mis-excluded R wave:
(a)n点处于两个R波之间,并且n点到前述两个R波时间间隔都大于前述两个R波的时间间隔的2/3;(a) Point n is between two R waves, and the time interval from point n to the two preceding R waves is greater than 2/3 of the time interval between the preceding two R waves;
(b)n点与前一个R波之间不存在多个无效的极大点;(b) There are no multiple invalid maximum points between point n and the previous R wave;
(c)n点到前一个R点和到后一个R点间隔都大于1/3秒;(c) The interval from point n to the previous point R and to the next point R is greater than 1/3 second;
(d)在差分信号上,n点幅度与旁边极小点幅度差大于n点前一个R波与旁边极小点幅度差0.1倍;(d) On the differential signal, the difference between the amplitude of point n and the minimum point next to it is greater than 0.1 times the difference between the amplitude of the R wave before point n and the minimum point next to it;
(e)在差分信号上,n点幅度与旁边极小点幅度差大于n点后一个R波与旁边极小点幅度差0.1倍;(e) On the differential signal, the amplitude difference between the n-point amplitude and the adjacent minimum point is greater than 0.1 times the amplitude difference between the R wave after n point and the adjacent minimum point;
4a)在近似R波位置周围±25点范围内寻找到真正R位置。4a) Find the true R position within ±25 points around the approximate R wave position.
具体地,可采用如下方法进行心电图其他时域特征提取:Specifically, the following methods can be used to extract other time-domain features of the ECG:
时域特征计算,分别计算QT间期、QRS波的斜率、ST段斜率以及相邻两个QRS波的间隔;Calculation of time-domain features, respectively calculating the QT interval, the slope of the QRS wave, the slope of the ST segment, and the interval between two adjacent QRS waves;
形态提取预处理,分别对P波、QRS波和T波的[start,end]范围内采样数据计算所有拐点,进一步计算出所有趋势转折点,并定位趋势谷与趋势顶;趋势谷是指那些趋势由下向上转变的转折点,趋势顶是指那些趋势由上向下转变的转折点;Form extraction preprocessing, respectively calculate all inflection points for the sampling data in the [start, end] range of P wave, QRS wave and T wave, and further calculate all trend turning points, and locate trend valleys and trend tops; trend valleys refer to those trends The turning point that changes from bottom to top, and the top of the trend refers to the turning point of those trends that change from top to bottom;
QRS波形态特征点计算,根据预处理步骤得到的趋势转折点,配合基线定位Q波、R波、S波、R’波和S’波,同时计算上述的Q波、R波、S波、R’波和S’波的振幅;Calculation of QRS wave morphological feature points, according to the trend turning point obtained in the preprocessing step, coordinate with the baseline to locate Q wave, R wave, S wave, R' wave and S' wave, and calculate the above-mentioned Q wave, R wave, S wave, R wave at the same time 'wave and S' wave amplitudes;
形态识别,根据形态提取预处理步骤的结果得到P波和T波的形态特征;根据QRS波形态特征点计算步骤的计算结果,识别QRS波形态模式,得到QRS波的形态特征。Morphological recognition, according to the results of the shape extraction preprocessing step to obtain the morphological features of the P wave and T wave; according to the calculation results of the QRS wave morphological feature point calculation step, identify the QRS wave morphological pattern, and obtain the morphological features of the QRS wave.
具体地,可通过如下方法进行脉搏波主波检测:Specifically, pulse wave main wave detection can be performed by the following methods:
1)、脉搏波先经1Hz~30Hz的带通滤波,去除低频和高频噪音;1) The pulse wave is first filtered by a 1Hz-30Hz band-pass filter to remove low-frequency and high-frequency noise;
2)、然后对滤波后信号进行幅值归一化;2), and then perform amplitude normalization on the filtered signal;
3)、计算信号幅值归一化后的香农能量,并进行低通滤波求取包络线;3) Calculate the Shannon energy after normalizing the signal amplitude, and perform low-pass filtering to obtain the envelope;
4)、对香农能量包络线信号进行希尔伯特变换;4) Hilbert transform the Shannon energy envelope signal;
5)、利用希尔伯特变换后的信号提取香农能量包络线的峰值点t,即希尔伯特变换后的信号从负到正过零点的值与香农能量包络线峰值点对应;5) Use the Hilbert transformed signal to extract the peak point t of the Shannon energy envelope, that is, the value of the Hilbert transformed signal from negative to positive zero crossing corresponds to the peak point of the Shannon energy envelope;
6)确定原始脉搏波信号中的真正主波位置,以香农能量峰值点的时间位置t为中心确定时间范围T,脉搏波原始信号在T范围找到真正的脉搏波主波位置;0.2<T<0.3。6) Determine the real main wave position in the original pulse wave signal, and determine the time range T centered on the time position t of the Shannon energy peak point, and find the real main pulse wave position in the T range of the original pulse wave signal; 0.2<T< 0.3.
当然,还可以采用其他已知的对心电图R波检测、其他时域特征提取,以及对脉搏波主波检测的方法,上述描述仅为举例说明,本发明对此不做限定。Of course, other known methods for detecting the R wave of the electrocardiogram, extracting other time-domain features, and detecting the main wave of the pulse wave can also be used. The above description is only for illustration, and the present invention is not limited thereto.
优选的,在提取心电图特征和脉搏波特征时,还可以通过采用频域特征分析方法,提取心电图和脉搏波的频域特征。频域特征分析方法包括但不限于傅里叶变换、小波变换、小波包和高阶统计量等频域特征分析方法,其方法选择根据具体情况确定。Preferably, when extracting the features of the electrocardiogram and the pulse wave, the frequency domain features of the electrocardiogram and the pulse wave can also be extracted by adopting a frequency domain feature analysis method. Frequency domain feature analysis methods include but are not limited to frequency domain feature analysis methods such as Fourier transform, wavelet transform, wavelet packet, and high-order statistics, and the method selection is determined according to the specific situation.
步骤S4对比分析时,以房颤疾病为例,在某一具体实施例中,心电图符合房扑和房颤特征,计算心率(RR间期)和脉率(脉搏波主波间期),发现房颤疾病脉率与心率不一致,而房扑则一致。In the comparative analysis of step S4, taking atrial fibrillation as an example, in a specific embodiment, the electrocardiogram conforms to the characteristics of atrial flutter and atrial fibrillation, and the heart rate (RR interval) and pulse rate (pulse wave main wave interval) are calculated, and it is found that In atrial fibrillation the pulse rate does not coincide with the heart rate, whereas in atrial flutter it does.
实施例2Example 2
本发明实施例2提供了一种心电图和脉搏波关系分析装置,该装置采用实施例1提供的方法进行心电图和脉搏波关系分析。本实施例中的心电图和脉搏波数据为针对某一具体已知的心血管疾病获取的心电图和脉搏波数据,例如但不限于房颤疾病,此处不一一列举。Embodiment 2 of the present invention provides a device for analyzing the relationship between electrocardiogram and pulse wave. The device adopts the method provided in embodiment 1 to analyze the relationship between electrocardiogram and pulse wave. The electrocardiogram and pulse wave data in this embodiment are the electrocardiogram and pulse wave data obtained for a specific known cardiovascular disease, such as but not limited to atrial fibrillation, which are not listed here.
如图2所示,该装置包括预处理模块1、可视化呈现模块2、心电图特征提取模块3、脉搏波特征提取模块4和对比分析模块6。其中,预处理模块1对获取的心电图和脉搏波信号分别进行预处理,包括基线平移和去噪;可视化呈现模块2对预处理后的心电图和脉搏波进行二维和/或三维图像的可视化呈现;心电图特征提取模块3提取心电图特征;脉搏波特征提取模块4提取脉搏波特征;对比分析模块6对心电图特征和脉搏波特征进行对比分析,得出心电图和脉搏波信号各特征之间的关系。As shown in FIG. 2 , the device includes a preprocessing module 1 , a visual presentation module 2 , an electrocardiogram
心电图特征主要包括R波和心电图其他时域特征,脉搏波特征主要包括主波和脉搏波其他时域特征,因此,心电图特征提取模块3还可以进一步包括心电图R波检测单元3a和心电图其他时域特征提取单元3b,其中,心电图R波检测单元3a进行心电图R波检测,提取心电图中R波位置;心电图其他时域特征提取单元3b进行心电图中除R波以外的其他时域特征检测,提取心电图中除R波以外的其他时域特征。ECG features mainly include R wave and other time domain features of ECG, pulse wave features mainly include main wave and other time domain features of pulse wave, therefore, ECG
脉搏波特征提取模块4还可以进一步包括脉搏波主波检测单元4a和脉搏波时域特征提取单元4b,其中,脉搏波主波检测单元4a进行脉搏波主波检测,提取脉搏波中主波位置;脉搏波时域特征提取单元4b进行脉搏波时域特征检测,包括:以主波位置为依据,进行脉搏波周期起始点检测,确定脉搏波周期,检测脉搏波周期是否为噪声周期,对非噪声周期检测重博前波、重博波位置幅值和主波及周期形态。The pulse wave
优选的,该装置还可以包括频域特征提取模块5,频域特征提取模块5通过频域特征分析方法,提取心电图和脉搏波的频域特征。Preferably, the device may further include a frequency-domain feature extraction module 5, which extracts the frequency-domain features of the electrocardiogram and pulse wave through a frequency-domain feature analysis method.
实施例3Example 3
本发明实施例3提供了一种心电图和脉搏波关系分析方法。本实施例中需要处理的心电图和脉搏波数据为针对某一未知疾病对象获取的心电图和脉搏波数据,该未知疾病对象可能患有心血管疾病,例如但不限于房颤疾病,此处不一一列举。该方法的实施依赖于实施例1,也即其需要在实施例1的基础上进行。
如图3所示,该方法包括如下步骤:As shown in Figure 3, the method includes the following steps:
步骤Q1:对获取的一未知疾病对象的心电图和脉搏波信号进行决策层融合;Step Q1: Decision-making fusion of the acquired ECG and pulse wave signals of an unknown disease subject;
步骤Q2:采用深度学习分类器,利用深层神经网络对未知疾病对象的心电图和脉搏波信号进行数据层和特征层融合;Step Q2: Using a deep learning classifier, using a deep neural network to perform data layer and feature layer fusion on the ECG and pulse wave signals of unknown disease subjects;
步骤Q3:对决策层融合、数据层和特征层融合后的心电图和脉搏波信号各特征之间的关系进行综合分析,并根据实施例1得到的心电图和脉搏波信号各特征之间的关系进行识别分类。Step Q3: Comprehensively analyze the relationship between the features of the electrocardiogram and the pulse wave signal after the fusion of the decision-making layer, the data layer and the feature layer, and perform a comprehensive analysis according to the relationship between the features of the electrocardiogram and the pulse wave signal obtained in Example 1 Identify categories.
下面对上述步骤进行详细说明:The above steps are described in detail below:
步骤Q1具体包括:Step Q1 specifically includes:
步骤Q1a:进行心电图R波检测,提取心电图中R波位置;Step Q1a: performing electrocardiogram R-wave detection, and extracting the R-wave position in the electrocardiogram;
步骤Q1b:进行心电图中除R波以外的其他时域特征检测,提取心电图中除R波以外的其他时域特征;Step Q1b: Perform other time-domain feature detection in the electrocardiogram except the R wave, and extract other time-domain features in the electrocardiogram except the R wave;
步骤Q1c:进行脉搏波主波检测,提取脉搏波中主波位置;Step Q1c: Perform pulse wave main wave detection, and extract the main wave position in the pulse wave;
步骤Q1d:进行脉搏波时域特征检测,包括:以主波位置为依据,进行脉搏波周期起始点检测,确定脉搏波周期,检测脉搏波周期是否为噪声周期,对非噪声周期检测重博前波、重博波位置幅值和主波及周期形态;Step Q1d: Perform pulse wave time-domain feature detection, including: based on the position of the main wave, detect the starting point of the pulse wave cycle, determine the pulse wave cycle, detect whether the pulse wave cycle is a noise cycle, and detect whether the pulse wave cycle is a non-noise cycle. wave, double wave position amplitude and main wave and cycle form;
步骤Q1e:根据实施例1得到的心电图和脉搏波信号各特征之间的关系,对步骤Q1a-Q1d中提取的心电图特征和脉搏波特征之间的关系进行规则推理。Step Q1e: According to the relationship between the features of the electrocardiogram and the pulse wave signal obtained in embodiment 1, rule-based reasoning is performed on the relationship between the features of the electrocardiogram and the pulse wave signal extracted in steps Q1a-Q1d.
其中,步骤Q1a-Q1d所提供的方法与实施例1中步骤S3a-S3d的方法类似。此处不再对其中的细节作详细描述,请参考实施例1。Wherein, the methods provided in steps Q1a-Q1d are similar to the methods in steps S3a-S3d in Embodiment 1. The details will not be described in detail here, please refer to Embodiment 1.
步骤Q2数据层和特征层融合可采用如下方法:在确定脉搏波的主波的位置后,计算脉搏波特征的直接特征值,并对分段的数据段进行卷积和取样,得到内部特征值,再结合直接特征值和内部特征值根据预定算法进行计算,得到分类结果。不提取各种准确率不太高、容易受噪声干扰的特征值(如:重搏前波),而是直接提取准确率很高的直接特征值(如:主波间期)纳入算法进行计算,可以提高最后分类准确率,输出更准确的脉搏波分类结果。可以理解的是,还可以采用其他已知的方法进行数据层和特征层融合,上述描述仅为举例说明,本发明对此不做限定。In step Q2, the fusion of the data layer and the feature layer can adopt the following method: after determining the position of the main wave of the pulse wave, calculate the direct eigenvalue of the pulse wave feature, and perform convolution and sampling on the segmented data segments to obtain the internal eigenvalue , combined with direct eigenvalues and internal eigenvalues to calculate according to a predetermined algorithm to obtain classification results. Instead of extracting various eigenvalues with low accuracy and being easily disturbed by noise (such as dicrotic front wave), direct eigenvalues with high accuracy (such as main wave interval) are directly extracted and included in the algorithm for calculation , can improve the final classification accuracy and output more accurate pulse wave classification results. It can be understood that other known methods can also be used to fuse the data layer and the feature layer, and the above description is only for illustration, and the present invention is not limited thereto.
实施例4Example 4
本发明实施例4提供了一种心电图和脉搏波关系分析装置。本实施例中需要处理的心电图和脉搏波数据为针对某一未知疾病对象获取的心电图和脉搏波数据,该未知疾病对象可能患有心血管疾病,例如但不限于房颤疾病,此处不一一列举。该装置的实现依赖于实施例2,也即其需要在实施例2的基础上进行。
如图4所示,该装置包括决策层融合模块7、数据层与特征层融合模块8和综合分析模块9,其中,决策层融合模块7对获取的一未知疾病对象的心电图和脉搏波信号进行决策层融合;数据层与特征层融合模块8包括一深度学习分类器,深度学习分类器利用深层神经网络对未知疾病对象的心电图和脉搏波信号进行数据层特征层融合;综合分析模块9对决策层融合、数据层和特征层融合后的心电图和脉搏波信号各特征之间的关系进行综合分析,并根据实施例2得到的心电图和脉搏波信号各特征之间的关系进行识别分类。As shown in Figure 4, the device includes a decision-making layer fusion module 7, a data layer and feature layer fusion module 8 and a comprehensive analysis module 9, wherein the decision-making layer fusion module 7 performs an electrocardiogram and a pulse wave signal on an acquired unknown disease object. Decision-making layer fusion; data layer and feature layer fusion module 8 includes a deep learning classifier, and the deep learning classifier utilizes deep neural network to carry out data layer feature layer fusion to the electrocardiogram and pulse wave signal of unknown disease object; comprehensive analysis module 9 pairs of decision-making Layer fusion, the relationship between the electrocardiogram and the pulse wave signal features after the fusion of the data layer and the feature layer is analyzed comprehensively, and the relationship between the electrocardiogram and the pulse wave signal features is identified and classified according to embodiment 2.
优选地,决策层融合模块7具体包括心电图R波检测单元7a、心电图其他时域特征提取单元7b、脉搏波主波检测单元7c、脉搏波时域特征提取单元7d和规则推理单元7e,其中,心电图R波检测单元7a进行心电图R波检测,提取心电图中R波位置;心电图其他时域特征提取单元7b进行心电图中除R波以外的其他时域特征检测,提取心电图中除R波以外的其他时域特征;脉搏波主波检测单元7c进行脉搏波主波检测,提取脉搏波中主波位置;脉搏波时域特征提取单元7d进行脉搏波时域特征检测,包括:以主波位置为依据,进行脉搏波周期起始点检测,确定脉搏波周期,检测脉搏波周期是否为噪声周期,对非噪声周期检测重博前波、重博波位置幅值和主波及周期形态;规则推理单元7e根据实施例2得到的心电图和脉搏波信号各特征之间的关系,对心电图R波检测单元7a、心电图其他时域特征提取单元7b、脉搏波主波检测单元7c和脉搏波时域特征提取单元7d提取的心电图特征和脉搏波特征之间的关系进行规则推理。Preferably, the decision-making layer fusion module 7 specifically includes an electrocardiogram R
其中,心电图R波检测单元7a、心电图其他时域特征提取单元7b、脉搏波主波检测单元7c和脉搏波时域特征提取单元7d的功能分别与实施例2中的心电图R波检测单元3a、心电图其他时域特征提取单元3b、脉搏波主波检测单元3c和脉搏波时域特征提取单元3d类似。具体实施时,可通过一个单元实现两种功能,也可以分别通过两个单元来实现两种功能,例如在本实施例的装置中只设置一个心电图R波检测单元,既检测已知疾病对象的心电图R波,也检测未知疾病对象的R波;也可以在本实施例的装置中分别设置两个心电图R波检测单元:心电图R波检测单元7a和心电图R波检测单元3a,分别检测未知疾病对象和已知疾病对象的心电图R波。Wherein, the functions of the electrocardiogram R
值得说明的是,上述装置和系统内的模块、单元之间的信息交互、执行过程等内容,由于与本发明的处理方法实施例基于同一构思,具体内容可参见本发明方法实施例中的叙述,此处不再赘述。It is worth noting that the information interaction and execution process between the above-mentioned devices and modules and units in the system are based on the same idea as the embodiment of the processing method of the present invention, and the specific content can refer to the description in the embodiment of the method of the present invention , which will not be repeated here.
本领域普通技术人员可以理解实施例的各种方法中的全部或部分步骤是可以通过程序来指令相关的硬件来完成,该程序可以存储于一计算机可读存储介质中,存储介质可以包括:只读存储器(ROM,Read Only Memory)、随机存取存储器(RAM,Random Access Memory)、磁盘或光盘等。Those of ordinary skill in the art can understand that all or part of the steps in the various methods of the embodiments can be completed by instructing related hardware through a program. The program can be stored in a computer-readable storage medium, and the storage medium can include: only Read memory (ROM, Read Only Memory), random access memory (RAM, Random Access Memory), disk or CD, etc.
以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。The above descriptions are only preferred embodiments of the present invention, and are not intended to limit the present invention. Any modifications, equivalent replacements and improvements made within the spirit and principles of the present invention should be included in the protection of the present invention. within range.
Claims (12)
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN2013102395557A CN103284702A (en) | 2013-06-17 | 2013-06-17 | Electrocardiogram and pulse wave relation analysis method and method and device of fusion analysis |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN2013102395557A CN103284702A (en) | 2013-06-17 | 2013-06-17 | Electrocardiogram and pulse wave relation analysis method and method and device of fusion analysis |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN103284702A true CN103284702A (en) | 2013-09-11 |
Family
ID=49086644
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN2013102395557A Pending CN103284702A (en) | 2013-06-17 | 2013-06-17 | Electrocardiogram and pulse wave relation analysis method and method and device of fusion analysis |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN103284702A (en) |
Cited By (34)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104739384A (en) * | 2013-12-30 | 2015-07-01 | 财团法人工业技术研究院 | Arterial wave analysis method and system |
| CN104939821A (en) * | 2014-03-27 | 2015-09-30 | 中国科学院苏州纳米技术与纳米仿生研究所 | Electrocardiogram waveform peak value detecting method and system |
| CN105708441A (en) * | 2014-12-03 | 2016-06-29 | 中国科学院苏州纳米技术与纳米仿生研究所 | Wearable finger stall and electrocardiogram and pulse condition acquisition terminal |
| CN105832533A (en) * | 2016-06-06 | 2016-08-10 | 厦门精图信息技术有限公司 | Intelligent wearable heart protector |
| CN106333648A (en) * | 2016-09-18 | 2017-01-18 | 京东方科技集团股份有限公司 | Sleep asphyxia monitoring method based on wearable device and wearable device |
| CN106388788A (en) * | 2016-09-27 | 2017-02-15 | 中国科学院微电子研究所 | Pulse wave diagnosis system with cardiac function information acquisition function |
| CN107106840A (en) * | 2014-10-07 | 2017-08-29 | 纽罗路普有限公司 | components that can be implanted |
| US9824287B2 (en) | 2015-09-29 | 2017-11-21 | Huami Inc. | Method, apparatus and system for biometric identification |
| CN107874753A (en) * | 2016-09-29 | 2018-04-06 | 中国科学院微电子研究所 | Pulse condition recognition method and device |
| US9948642B2 (en) | 2015-09-29 | 2018-04-17 | Anhui Huami Information Technology Co., Ltd. | Multi-modal biometric identification |
| CN108403094A (en) * | 2018-03-23 | 2018-08-17 | 中国科学院微电子研究所 | A method for identifying pulse wave peaks |
| CN108511055A (en) * | 2017-02-27 | 2018-09-07 | 中国科学院苏州纳米技术与纳米仿生研究所 | Ventricular premature beat identifying system and method based on Multiple Classifier Fusion and diagnostic rule |
| CN108670209A (en) * | 2018-03-29 | 2018-10-19 | 中国科学院微电子研究所 | A method and system for automatically identifying pulse conditions in traditional Chinese medicine |
| CN109044299A (en) * | 2018-09-25 | 2018-12-21 | 北方民族大学 | A kind of physiological signal collection device and method of cardiovascular health detection |
| CN109157202A (en) * | 2018-09-18 | 2019-01-08 | 东北大学 | A kind of cardiovascular disease early warning system based on more physiological signal depth integrations |
| CN109871808A (en) * | 2019-02-21 | 2019-06-11 | 天津惊帆科技有限公司 | Atrial fibrillation model training and detecting method and device |
| US10426364B2 (en) | 2015-10-27 | 2019-10-01 | Cardiologs Technologies Sas | Automatic method to delineate or categorize an electrocardiogram |
| US10467548B2 (en) | 2015-09-29 | 2019-11-05 | Huami Inc. | Method, apparatus and system for biometric identification |
| CN111091092A (en) * | 2019-12-16 | 2020-05-01 | 桂林电子科技大学 | A Data Augmentation Method Based on Optimizing Convolutional Neural Network Classification |
| US10779744B2 (en) | 2015-10-27 | 2020-09-22 | Cardiologs Technologies Sas | Automatic method to delineate or categorize an electrocardiogram |
| CN111839494A (en) * | 2020-09-04 | 2020-10-30 | 广东电网有限责任公司电力科学研究院 | Heart rate monitoring method and system |
| US10827938B2 (en) | 2018-03-30 | 2020-11-10 | Cardiologs Technologies Sas | Systems and methods for digitizing electrocardiograms |
| CN111956211A (en) * | 2020-07-29 | 2020-11-20 | 鲁东大学 | Automatic detection method for atrial fibrillation of single lead electrocardiosignal |
| CN112004470A (en) * | 2018-04-12 | 2020-11-27 | 日本电信电话株式会社 | Maximum heart rate estimation method and device |
| CN113499059A (en) * | 2021-06-01 | 2021-10-15 | 武汉理工大学 | BCG signal processing system and method based on optical fiber sensing non-contact |
| CN114081502A (en) * | 2021-11-11 | 2022-02-25 | 中国科学院深圳先进技术研究院 | A non-invasive heart disease diagnosis method and device based on machine learning |
| US11331034B2 (en) | 2015-10-27 | 2022-05-17 | Cardiologs Technologies Sas | Automatic method to delineate or categorize an electrocardiogram |
| CN114916923A (en) * | 2022-05-19 | 2022-08-19 | 南京中医药大学 | Electro-mechanical interconnection electrocardio pulse signal analysis method and system |
| CN115590469A (en) * | 2021-07-09 | 2023-01-13 | 特拉量子股份公司(Ch) | Method and system for converting physiological signals |
| US11672464B2 (en) | 2015-10-27 | 2023-06-13 | Cardiologs Technologies Sas | Electrocardiogram processing system for delineation and classification |
| US11678831B2 (en) | 2020-08-10 | 2023-06-20 | Cardiologs Technologies Sas | Electrocardiogram processing system for detecting and/or predicting cardiac events |
| CN116763311A (en) * | 2023-06-26 | 2023-09-19 | 山东大学 | Depression state identifying system based on multi-level characteristic fusion |
| US11826150B2 (en) | 2017-08-25 | 2023-11-28 | Koninklijke Philips N.V. | User interface for analysis of electrocardiograms |
| US12016694B2 (en) | 2019-02-04 | 2024-06-25 | Cardiologs Technologies Sas | Electrocardiogram processing system for delineation and classification |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101897578A (en) * | 2010-06-30 | 2010-12-01 | 重庆大学 | A beat-by-beat segmentation method of arterial pressure signal |
| CN102008296A (en) * | 2010-12-24 | 2011-04-13 | 吉林大学 | Device and method for measuring arterial blood pressures based on pulse wave signals and electrocardiosignals |
| CN102270264A (en) * | 2010-06-04 | 2011-12-07 | 中国科学院深圳先进技术研究院 | Physiological signal quality evaluation system and method |
| CN102824166A (en) * | 2012-07-05 | 2012-12-19 | 华东师范大学 | Device for sorting treatment of pulse wave |
-
2013
- 2013-06-17 CN CN2013102395557A patent/CN103284702A/en active Pending
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102270264A (en) * | 2010-06-04 | 2011-12-07 | 中国科学院深圳先进技术研究院 | Physiological signal quality evaluation system and method |
| CN101897578A (en) * | 2010-06-30 | 2010-12-01 | 重庆大学 | A beat-by-beat segmentation method of arterial pressure signal |
| CN102008296A (en) * | 2010-12-24 | 2011-04-13 | 吉林大学 | Device and method for measuring arterial blood pressures based on pulse wave signals and electrocardiosignals |
| CN102824166A (en) * | 2012-07-05 | 2012-12-19 | 华东师范大学 | Device for sorting treatment of pulse wave |
Non-Patent Citations (1)
| Title |
|---|
| 孙松松: "心电脉搏相关性及其信息融合方法研究", 《中国优秀硕士学位论文全文数据库信息科技辑》, no. 10, 15 October 2012 (2012-10-15), pages 136 - 464 * |
Cited By (46)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104739384A (en) * | 2013-12-30 | 2015-07-01 | 财团法人工业技术研究院 | Arterial wave analysis method and system |
| CN104939821A (en) * | 2014-03-27 | 2015-09-30 | 中国科学院苏州纳米技术与纳米仿生研究所 | Electrocardiogram waveform peak value detecting method and system |
| CN107106840A (en) * | 2014-10-07 | 2017-08-29 | 纽罗路普有限公司 | components that can be implanted |
| CN105708441A (en) * | 2014-12-03 | 2016-06-29 | 中国科学院苏州纳米技术与纳米仿生研究所 | Wearable finger stall and electrocardiogram and pulse condition acquisition terminal |
| US10467548B2 (en) | 2015-09-29 | 2019-11-05 | Huami Inc. | Method, apparatus and system for biometric identification |
| US9946942B2 (en) | 2015-09-29 | 2018-04-17 | Huami Inc. | Method, apparatus and system for biometric identification |
| US9948642B2 (en) | 2015-09-29 | 2018-04-17 | Anhui Huami Information Technology Co., Ltd. | Multi-modal biometric identification |
| US9824287B2 (en) | 2015-09-29 | 2017-11-21 | Huami Inc. | Method, apparatus and system for biometric identification |
| US11134880B2 (en) | 2015-10-27 | 2021-10-05 | Cardiologs Technologies Sas | Automatic method to delineate or categorize an electrocardiogram |
| US10758139B2 (en) | 2015-10-27 | 2020-09-01 | Cardiologs Technologies Sas | Automatic method to delineate or categorize an electrocardiogram |
| US11331034B2 (en) | 2015-10-27 | 2022-05-17 | Cardiologs Technologies Sas | Automatic method to delineate or categorize an electrocardiogram |
| US11147500B2 (en) | 2015-10-27 | 2021-10-19 | Cardiologs Technologies Sas | Electrocardiogram processing system for delineation and classification |
| US10959660B2 (en) | 2015-10-27 | 2021-03-30 | Cardiologs Technologies Sas | Electrocardiogram processing system for delineation and classification |
| US10426364B2 (en) | 2015-10-27 | 2019-10-01 | Cardiologs Technologies Sas | Automatic method to delineate or categorize an electrocardiogram |
| US11672464B2 (en) | 2015-10-27 | 2023-06-13 | Cardiologs Technologies Sas | Electrocardiogram processing system for delineation and classification |
| US10779744B2 (en) | 2015-10-27 | 2020-09-22 | Cardiologs Technologies Sas | Automatic method to delineate or categorize an electrocardiogram |
| CN105832533A (en) * | 2016-06-06 | 2016-08-10 | 厦门精图信息技术有限公司 | Intelligent wearable heart protector |
| CN106333648A (en) * | 2016-09-18 | 2017-01-18 | 京东方科技集团股份有限公司 | Sleep asphyxia monitoring method based on wearable device and wearable device |
| CN106388788A (en) * | 2016-09-27 | 2017-02-15 | 中国科学院微电子研究所 | Pulse wave diagnosis system with cardiac function information acquisition function |
| CN107874753A (en) * | 2016-09-29 | 2018-04-06 | 中国科学院微电子研究所 | Pulse condition recognition method and device |
| CN108511055B (en) * | 2017-02-27 | 2021-10-12 | 中国科学院苏州纳米技术与纳米仿生研究所 | Ventricular premature beat recognition system and method based on classifier fusion and diagnosis rules |
| CN108511055A (en) * | 2017-02-27 | 2018-09-07 | 中国科学院苏州纳米技术与纳米仿生研究所 | Ventricular premature beat identifying system and method based on Multiple Classifier Fusion and diagnostic rule |
| US11826150B2 (en) | 2017-08-25 | 2023-11-28 | Koninklijke Philips N.V. | User interface for analysis of electrocardiograms |
| CN108403094A (en) * | 2018-03-23 | 2018-08-17 | 中国科学院微电子研究所 | A method for identifying pulse wave peaks |
| CN108670209A (en) * | 2018-03-29 | 2018-10-19 | 中国科学院微电子研究所 | A method and system for automatically identifying pulse conditions in traditional Chinese medicine |
| US10827938B2 (en) | 2018-03-30 | 2020-11-10 | Cardiologs Technologies Sas | Systems and methods for digitizing electrocardiograms |
| CN112004470A (en) * | 2018-04-12 | 2020-11-27 | 日本电信电话株式会社 | Maximum heart rate estimation method and device |
| CN112004470B (en) * | 2018-04-12 | 2024-06-04 | 日本电信电话株式会社 | Maximum heart rate estimation method and device |
| CN109157202A (en) * | 2018-09-18 | 2019-01-08 | 东北大学 | A kind of cardiovascular disease early warning system based on more physiological signal depth integrations |
| CN109044299B (en) * | 2018-09-25 | 2021-03-30 | 北方民族大学 | Physiological signal acquisition device and method for cardiovascular health detection |
| CN109044299A (en) * | 2018-09-25 | 2018-12-21 | 北方民族大学 | A kind of physiological signal collection device and method of cardiovascular health detection |
| US12016694B2 (en) | 2019-02-04 | 2024-06-25 | Cardiologs Technologies Sas | Electrocardiogram processing system for delineation and classification |
| CN109871808B (en) * | 2019-02-21 | 2021-03-02 | 天津惊帆科技有限公司 | Atrial fibrillation model training and detecting method and device |
| CN109871808A (en) * | 2019-02-21 | 2019-06-11 | 天津惊帆科技有限公司 | Atrial fibrillation model training and detecting method and device |
| CN111091092A (en) * | 2019-12-16 | 2020-05-01 | 桂林电子科技大学 | A Data Augmentation Method Based on Optimizing Convolutional Neural Network Classification |
| CN111956211A (en) * | 2020-07-29 | 2020-11-20 | 鲁东大学 | Automatic detection method for atrial fibrillation of single lead electrocardiosignal |
| US11678831B2 (en) | 2020-08-10 | 2023-06-20 | Cardiologs Technologies Sas | Electrocardiogram processing system for detecting and/or predicting cardiac events |
| CN111839494A (en) * | 2020-09-04 | 2020-10-30 | 广东电网有限责任公司电力科学研究院 | Heart rate monitoring method and system |
| CN113499059A (en) * | 2021-06-01 | 2021-10-15 | 武汉理工大学 | BCG signal processing system and method based on optical fiber sensing non-contact |
| CN115590469A (en) * | 2021-07-09 | 2023-01-13 | 特拉量子股份公司(Ch) | Method and system for converting physiological signals |
| WO2023082402A1 (en) * | 2021-11-11 | 2023-05-19 | 中国科学院深圳先进技术研究院 | Non-invasive heart disease diagnosis method and device based on machine learning |
| CN114081502A (en) * | 2021-11-11 | 2022-02-25 | 中国科学院深圳先进技术研究院 | A non-invasive heart disease diagnosis method and device based on machine learning |
| CN114916923B (en) * | 2022-05-19 | 2023-04-11 | 南京中医药大学 | Electro-mechanical interconnection electrocardio pulse signal analysis method and system |
| CN114916923A (en) * | 2022-05-19 | 2022-08-19 | 南京中医药大学 | Electro-mechanical interconnection electrocardio pulse signal analysis method and system |
| CN116763311A (en) * | 2023-06-26 | 2023-09-19 | 山东大学 | Depression state identifying system based on multi-level characteristic fusion |
| CN116763311B (en) * | 2023-06-26 | 2026-01-30 | 山东大学 | Depression State Recognition System Based on Multi-Level Feature Fusion |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN103284702A (en) | Electrocardiogram and pulse wave relation analysis method and method and device of fusion analysis | |
| Chen et al. | Atrial fibrillation detection using a feedforward neural network | |
| Sahoo et al. | De-noising of ECG signal and QRS detection using Hilbert transform and adaptive thresholding | |
| Gupta | Wavelet transform and vector machines as emerging tools for computational medicine | |
| CN111481192B (en) | A method for detecting R wave of ECG signal based on improved U-Net | |
| CN104173043B (en) | Electrocardio data analysis method suitable for mobile platform | |
| CN105286853B (en) | Disease detection method and device based on wearable device, wearable device | |
| CN105411565B (en) | Heart rate variability tagsort method based on broad sense multi-scale wavelet entropy | |
| Xu et al. | Rule-based method for morphological classification of ST segment in ECG signals | |
| CN106473750A (en) | Personal identification method based on photoplethysmographic optimal period waveform | |
| Orhan | Real-time CHF detection from ECG signals using a novel discretization method | |
| CN107688553B (en) | Method for detecting electrocardiographic waveform characteristics based on wavelet transform and logistic regression algorithm | |
| CN110826631B (en) | Cardiac arrhythmia disease classification method, device, storage medium and processor | |
| Tseng et al. | Sliding large kernel of deep learning algorithm for mobile electrocardiogram diagnosis | |
| Beyramienanlou et al. | An Efficient Teager Energy Operator‐Based Automated QRS Complex Detection | |
| Zhong et al. | Maternal ECG removal using short time Fourier transform and convolutional auto-encoder | |
| CN116392136A (en) | A disease early warning method, device, equipment and medium | |
| Dash et al. | A novel approach for detecting fetal QRS and estimating fetal heart rate from abdominal ECG using EMD and wavelet decomposition | |
| Hussein et al. | Detection of electrocardiogram QRS complex based on modified adaptive threshold | |
| JP7199673B2 (en) | Radar heartbeat detection method and system | |
| Zhang et al. | Adaptive pulse wave extraction method for radial artery video | |
| CN118749990A (en) | Cardiopulmonary coupling feature extraction method and system based on EMD and dMI | |
| CN106963368A (en) | QRS complex detection method and device | |
| CN102715897A (en) | Method for determining average heart rates of fetuses on basis of fetal electrocardio spectrograms | |
| Gao et al. | Research on Heart Rate Signal Recognition Based on Deep Learning |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| C02 | Deemed withdrawal of patent application after publication (patent law 2001) | ||
| WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20130911 |