[go: up one dir, main page]

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 PDF

Info

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
Application number
CN2013102395557A
Other languages
Chinese (zh)
Inventor
胡晓娟
董军
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Suzhou Institute of Nano Tech and Nano Bionics of CAS
Original Assignee
Suzhou Institute of Nano Tech and Nano Bionics of CAS
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Suzhou Institute of Nano Tech and Nano Bionics of CAS filed Critical Suzhou Institute of Nano Tech and Nano Bionics of CAS
Priority to CN2013102395557A priority Critical patent/CN103284702A/en
Publication of CN103284702A publication Critical patent/CN103284702A/en
Pending legal-status Critical Current

Links

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

The invention relates to the technical field of measurement for diagnosis, and provides an electrocardiogram and pulse wave relation analysis method and a method and device of fusion analysis based on the electrocardiogram and pulse wave relation analysis method. The electrocardiogram and pulse wave relation analysis method includes the steps that step1, obtained electrocardiogram signals and pulse wave signals of a known patient are preprocessed respectively, and the preprocess comprises base line horizontal-moving and denoising; step2, the preprocessed electrocardiogram and the pulse wave are visually showed in a two-dimensional and/or three-dimensional image; step3, characteristics of the electrocardiogram and the pulse wave are extracted; step4, the characteristics of the electrocardiogram and the pulse wave are compared and analyzed to obtain the mutual relation of the characteristics of the electrocardiogram and the pulse wave, and fusion analysis inspection is carried out according to the mutual relation of the characteristics. Compared with analysis of singular signals, the accuracy rate of classification of results is improved.

Description

心电图和脉搏波关系分析及融合分析的方法和装置Method and device for relationship analysis and fusion analysis of electrocardiogram and pulse wave

【技术领域】【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 Embodiment 3 of the present invention;

图4是本发明实施例4一种心电图和脉搏波融合分析装置结构框图。Fig. 4 is a structural block diagram of an electrocardiogram and pulse wave fusion analysis device according to Embodiment 4 of the present invention.

【具体实施方式】【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 feature extraction module 3 , a pulse wave feature extraction module 4 and a comparative analysis module 6 . Among them, the preprocessing module 1 performs preprocessing on the acquired ECG and pulse wave signals respectively, including baseline translation and denoising; the visualization presentation module 2 performs two-dimensional and/or three-dimensional visual presentation of the preprocessed ECG and pulse wave signals The electrocardiogram feature extraction module 3 extracts the electrocardiogram feature; the pulse wave feature extraction module 4 extracts the pulse wave feature; the comparative analysis module 6 compares and analyzes the electrocardiogram feature and the pulse wave feature, and draws the relationship between each feature of the electrocardiogram and the pulse wave signal.

心电图特征主要包括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 feature extraction module 3 can further include ECG R wave detection unit 3a and other time domain features of ECG Feature extraction unit 3b, wherein, electrocardiogram R wave detection unit 3a detects electrocardiogram R wave, extracts R wave position in electrocardiogram; other time domain feature extraction unit 3b of electrocardiogram performs other time domain feature detection except R wave in electrocardiogram, extracts electrocardiogram other temporal features in addition to the R wave.

脉搏波特征提取模块4还可以进一步包括脉搏波主波检测单元4a和脉搏波时域特征提取单元4b,其中,脉搏波主波检测单元4a进行脉搏波主波检测,提取脉搏波中主波位置;脉搏波时域特征提取单元4b进行脉搏波时域特征检测,包括:以主波位置为依据,进行脉搏波周期起始点检测,确定脉搏波周期,检测脉搏波周期是否为噪声周期,对非噪声周期检测重博前波、重博波位置幅值和主波及周期形态。The pulse wave feature extraction module 4 can further include a pulse wave main wave detection unit 4a and a pulse wave time-domain feature extraction unit 4b, wherein the pulse wave main wave detection unit 4a detects the pulse wave main wave and extracts the main wave position in the pulse wave The pulse wave time-domain feature extraction unit 4b carries out pulse wave time-domain feature detection, including: taking the position of the main wave as the basis, carrying out pulse wave cycle starting point detection, determining the pulse wave cycle, and detecting whether the pulse wave cycle is a noise cycle, for non The noise cycle detects the front wave of the dichroic wave, the position and amplitude of the dichotomic wave, the main wave and the periodic shape.

优选的,该装置还可以包括频域特征提取模块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的基础上进行。Embodiment 3 of the present invention provides a method for analyzing the relationship between electrocardiogram and pulse wave. The electrocardiogram and pulse wave data that need to be processed in this embodiment are the electrocardiogram and pulse wave data obtained for a subject with an unknown disease. enumerate. The implementation of this method relies on Example 1, that is, it needs to be carried out on the basis of Example 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的基础上进行。Embodiment 4 of the present invention provides a device for analyzing the relationship between electrocardiogram and pulse wave. The electrocardiogram and pulse wave data that need to be processed in this embodiment are the electrocardiogram and pulse wave data obtained for a subject with an unknown disease. enumerate. The implementation of the device depends on Embodiment 2, that is, it needs to be performed on the basis of Embodiment 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 wave detection unit 7a, other time-domain feature extraction units 7b of the electrocardiogram, a main pulse wave detection unit 7c, a pulse wave time-domain feature extraction unit 7d and a rule reasoning unit 7e, wherein, The electrocardiogram R wave detection unit 7a detects the electrocardiogram R wave, and extracts the R wave position in the electrocardiogram; the other time domain feature extraction unit 7b of the electrocardiogram performs other time domain feature detection except the R wave in the electrocardiogram, and extracts other time domain features except the R wave in the electrocardiogram. Time-domain features; pulse wave main wave detection unit 7c detects pulse wave main wave and extracts main wave position in pulse wave; pulse wave time-domain feature extraction unit 7d performs pulse wave time-domain feature detection, including: based on main wave position , to detect the starting point of the pulse wave cycle, to determine the pulse wave cycle, to detect whether the pulse wave cycle is a noise cycle, to detect the non-noise cycle of the front wave of the pulse wave, the position amplitude of the pulse wave, the main wave and the cycle shape; the rule reasoning unit 7e according to The relationship between the electrocardiogram obtained in embodiment 2 and each feature of the pulse wave signal, for the electrocardiogram R wave detection unit 7a, other time domain feature extraction units 7b of the electrocardiogram, pulse wave main wave detection unit 7c and pulse wave time domain feature extraction unit 7d The relationship between the extracted ECG features and pulse wave features is used for rule reasoning.

其中,心电图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 wave detection unit 7a, other time domain feature extraction units 7b of the electrocardiogram, the pulse wave main wave detection unit 7c and the pulse wave time domain feature extraction unit 7d are respectively the same as those of the electrocardiogram R wave detection unit 3a, the pulse wave time domain feature extraction unit 7d in embodiment 2. The other time-domain feature extraction unit 3b of the electrocardiogram, the main pulse wave detection unit 3c and the pulse wave time-domain feature extraction unit 3d are similar. During specific implementation, two functions can be realized by one unit, and two functions can also be realized by two units respectively. The electrocardiogram R wave also detects the R wave of the unknown disease object; two electrocardiogram R wave detection units can also be respectively set in the device of this embodiment: the electrocardiogram R wave detection unit 7a and the electrocardiogram R wave detection unit 3a detect unknown diseases respectively ECG R-waves of subjects and subjects with known disease.

值得说明的是,上述装置和系统内的模块、单元之间的信息交互、执行过程等内容,由于与本发明的处理方法实施例基于同一构思,具体内容可参见本发明方法实施例中的叙述,此处不再赘述。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)

1. an electrocardiogram and pulse wave relationship analysis method is characterized in that, comprising:
Step S1: electrocardiogram and pulse wave signal to the known disease object that obtains carry out pretreatment respectively, comprise baseline translation and denoising;
Step S2: pretreated electrocardiogram and pulse wave are carried out the visual of two dimension and/or 3-D view present;
Step S3: extract Characteristics of electrocardiogram and pulse wave feature;
Step S4: Characteristics of electrocardiogram and pulse wave feature are analyzed, draw the relation between each feature of electrocardiogram and pulse wave signal.
2. electrocardiogram according to claim 1 and pulse wave relationship analysis method is characterized in that, among the described step S3, the method for extracting Characteristics of electrocardiogram specifically comprises:
Step S3a: carry out ECG R wave and detect, extract R ripple position in the electrocardiogram;
Step S3b: carry out other temporal signatures detections except the R ripple in the electrocardiogram, extract other temporal signatures except the R ripple in the electrocardiogram;
The method of extracting the pulse wave feature specifically comprises:
Step S3c: carry out pulse wave master ripple and detect, extract main ripple position in the pulse wave;
Step S3d: carry out the pulse wave temporal signatures and detect, comprise: be foundation with described main ripple position, carrying out the pulse wave cycle starting point detects, determine pulse wave cycle, whether be noise periods, non-noise periods is detected heavily rich prewave, heavily rich ripple position amplitude and the main cycle form that involves if detecting described pulse wave cycle.
3. electrocardiogram according to claim 2 and pulse wave relationship analysis method is characterized in that, in described step S3d, whether the described pulse wave cycle of described detection is that the step of noise periods is specially:
Obtain average pulse period of wave;
Based on described average pulse period of wave, calculate each through the cosine-shaped similarity of the pulse wave cycle after interpolation, with the similarity value of average pulse period of wave be noise periods less than the pulse wave cycle of a setting value S.
4. electrocardiogram according to claim 2 and pulse wave relationship analysis method is characterized in that, among the described step S3, the method for extracting Characteristics of electrocardiogram and pulse wave feature also comprises:
Adopt the frequency domain character analytical method, extract the frequency domain character of electrocardiogram and pulse wave.
5. an electrocardiogram and pulse wave relationship analysis apparatus is characterized in that, comprising:
Pre-processing module is used for electrocardiogram and the pulse wave signal that obtains carried out pretreatment respectively, comprises baseline translation and denoising;
The visual module that presents is used for that pretreated electrocardiogram and pulse wave are carried out the visual of two dimension and/or 3-D view and presents;
The Characteristics of electrocardiogram extraction module is used for extracting Characteristics of electrocardiogram;
The pulse wave characteristic extracting module is used for extracting the pulse wave feature;
The relative analysis module is used for Characteristics of electrocardiogram and pulse wave feature are analyzed, and draws the relation between each feature of electrocardiogram and pulse wave signal.
6. electrocardiogram according to claim 5 and pulse wave relationship analysis apparatus is characterized in that, described Characteristics of electrocardiogram extraction module specifically comprises:
The ECG R wave detecting unit is used for carrying out ECG R wave and detects, and extracts R ripple position in the electrocardiogram;
Other temporal signatures extraction units of electrocardiogram are used for carrying out electrocardiogram other temporal signatures except the R ripple and detect, and extract other temporal signatures except the R ripple in the electrocardiogram;
Described pulse wave characteristic extracting module specifically comprises:
Pulse wave master ripple detecting unit is used for carrying out pulse wave master ripple and detects, and extracts main ripple position in the pulse wave;
Pulse wave temporal signatures extraction unit, being used for carrying out the pulse wave temporal signatures detects, comprise: be foundation with described main ripple position, carrying out the pulse wave cycle starting point detects, determine pulse wave cycle, whether be noise periods, non-noise periods is detected heavily rich prewave, heavily rich ripple position amplitude and the main cycle form that involves if detecting described pulse wave cycle.
7. electrocardiogram according to claim 6 and pulse wave relationship analysis apparatus is characterized in that, whether described pulse wave temporal signatures extraction unit detects described pulse wave cycle is that the method for noise periods is specially:
Obtain average pulse period of wave;
Based on described average pulse period of wave, calculate each through the cosine-shaped similarity of the pulse wave cycle after interpolation, with the similarity value of average pulse period of wave be noise periods less than the pulse wave cycle of a setting value S.
8. electrocardiogram according to claim 5 and pulse wave relationship analysis apparatus is characterized in that, described device also comprises:
The frequency domain character extraction module is used for extracting the frequency domain character of electrocardiogram and pulse wave by the frequency domain character analytical method.
9. an electrocardiogram and pulse wave convergence analysis method is characterized in that, comprising:
Step Q1: the electrocardiogram of the unidentified illness object that obtains and pulse wave signal are carried out decision-making level merge;
Step Q2: adopt degree of deep learning classification device, utilize the deep layer neutral net to the electrocardiogram of described unidentified illness object and pulse wave signal carries out data Layer and characteristic layer merges;
Step Q3: the electrocardiogram after the fusion of described decision-making level, data Layer and the characteristic layer fusion and the relation between each feature of pulse wave signal are carried out analysis-by-synthesis, and carry out discriminator according to the relation between each electrocardiogram that obtains of claim 1-4 and each feature of pulse wave signal.
10. electrocardiogram according to claim 9 and pulse wave convergence analysis method is characterized in that described step Q1 specifically comprises:
Step Q1a: carry out ECG R wave and detect, extract R ripple position in the electrocardiogram;
Step Q1b: carry out other temporal signatures detections except the R ripple in the electrocardiogram, extract other temporal signatures except the R ripple in the electrocardiogram;
Step Q1c: carry out pulse wave master ripple and detect, extract main ripple position in the pulse wave;
Step Q1d: carry out the pulse wave temporal signatures and detect, comprise: be foundation with described main ripple position, carrying out the pulse wave cycle starting point detects, determine pulse wave cycle, whether be noise periods, non-noise periods is detected heavily rich prewave, heavily rich ripple position amplitude and the main cycle form that involves if detecting described pulse wave cycle;
Step Q1e: according to the relation between each electrocardiogram that obtains of claim 1-4 and each feature of pulse wave signal, the Characteristics of electrocardiogram extracted among the step Q1a-Q1d and the relation between the pulse wave feature are carried out rule-based reasoning.
11. an electrocardiogram and pulse wave convergence analysis device is characterized in that, comprising:
Decision-making level's Fusion Module is used for that the electrocardiogram of the unidentified illness object that obtains and pulse wave signal are carried out decision-making level and merges;
Data Layer and characteristic layer Fusion Module comprise a degree of deep learning classification device, and described degree of deep learning classification device utilizes the deep layer neutral net to the electrocardiogram of described unidentified illness object and pulse wave signal carries out data Layer and characteristic layer merges;
The analysis-by-synthesis module, be used for the electrocardiogram after the fusion of described decision-making level, data Layer and the characteristic layer fusion and the relation between each feature of pulse wave signal are carried out analysis-by-synthesis, and carry out discriminator according to the relation between each electrocardiogram that obtains of claim 5-8 and each feature of pulse wave signal.
12. electrocardiogram according to claim 11 and pulse wave convergence analysis device is characterized in that, described decision-making level Fusion Module specifically comprises:
The ECG R wave detecting unit is used for carrying out ECG R wave and detects, and extracts R ripple position in the electrocardiogram;
Other temporal signatures extraction units of electrocardiogram are used for carrying out electrocardiogram other temporal signatures except the R ripple and detect, and extract other temporal signatures except the R ripple in the electrocardiogram;
Pulse wave master ripple detecting unit is used for carrying out pulse wave master ripple and detects, and extracts main ripple position in the pulse wave;
Pulse wave temporal signatures extraction unit, being used for carrying out the pulse wave temporal signatures detects, comprise: be foundation with described main ripple position, carrying out the pulse wave cycle starting point detects, determine pulse wave cycle, whether be noise periods, non-noise periods is detected heavily rich prewave, heavily rich ripple position amplitude and the main cycle form that involves if detecting described pulse wave cycle;
The rule-based reasoning unit, be used for according to the relation between each electrocardiogram that obtains of claim 5-8 and each feature of pulse wave signal, ECG R wave detecting unit, other temporal signatures extraction units of electrocardiogram, pulse wave master ripple detecting unit and the Characteristics of electrocardiogram of pulse wave temporal signatures extraction unit extraction and the relation between the pulse wave feature are carried out rule-based reasoning.
CN2013102395557A 2013-06-17 2013-06-17 Electrocardiogram and pulse wave relation analysis method and method and device of fusion analysis Pending CN103284702A (en)

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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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

Patent Citations (4)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
Title
孙松松: "心电脉搏相关性及其信息融合方法研究", 《中国优秀硕士学位论文全文数据库信息科技辑》, no. 10, 15 October 2012 (2012-10-15), pages 136 - 464 *

Cited By (46)

* Cited by examiner, † Cited by third party
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