CN115270848A - 一种ppg与ecg自动转换智能算法、存储介质和计算机系统 - Google Patents
一种ppg与ecg自动转换智能算法、存储介质和计算机系统 Download PDFInfo
- Publication number
- CN115270848A CN115270848A CN202210688131.8A CN202210688131A CN115270848A CN 115270848 A CN115270848 A CN 115270848A CN 202210688131 A CN202210688131 A CN 202210688131A CN 115270848 A CN115270848 A CN 115270848A
- Authority
- CN
- China
- Prior art keywords
- signal
- ecg
- ppg
- signals
- neural network
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000006243 chemical reaction Methods 0.000 title claims abstract description 66
- 238000003860 storage Methods 0.000 title claims abstract description 33
- 238000012549 training Methods 0.000 claims abstract description 37
- 238000013473 artificial intelligence Methods 0.000 claims abstract description 23
- 238000007781 pre-processing Methods 0.000 claims abstract description 22
- 238000005520 cutting process Methods 0.000 claims abstract description 15
- 230000002452 interceptive effect Effects 0.000 claims abstract description 7
- 230000010354 integration Effects 0.000 claims abstract description 4
- 230000001360 synchronised effect Effects 0.000 claims abstract description 4
- 238000013186 photoplethysmography Methods 0.000 claims description 119
- 238000013528 artificial neural network Methods 0.000 claims description 60
- 230000006870 function Effects 0.000 claims description 40
- 238000000034 method Methods 0.000 claims description 25
- 238000005457 optimization Methods 0.000 claims description 24
- 238000005070 sampling Methods 0.000 claims description 22
- 238000012795 verification Methods 0.000 claims description 18
- 238000010606 normalization Methods 0.000 claims description 15
- 238000004590 computer program Methods 0.000 claims description 13
- ORILYTVJVMAKLC-UHFFFAOYSA-N Adamantane Natural products C1C(C2)CC3CC1CC2C3 ORILYTVJVMAKLC-UHFFFAOYSA-N 0.000 claims description 12
- 238000012952 Resampling Methods 0.000 claims description 12
- 238000001914 filtration Methods 0.000 claims description 12
- 238000009826 distribution Methods 0.000 claims description 6
- 238000012360 testing method Methods 0.000 claims description 6
- 238000012935 Averaging Methods 0.000 claims description 3
- 238000010276 construction Methods 0.000 claims description 3
- 238000009795 derivation Methods 0.000 claims description 3
- 238000002372 labelling Methods 0.000 claims description 3
- 230000008569 process Effects 0.000 claims description 3
- 230000000694 effects Effects 0.000 description 5
- 210000005036 nerve Anatomy 0.000 description 4
- 230000017531 blood circulation Effects 0.000 description 3
- 230000008859 change Effects 0.000 description 3
- 238000013507 mapping Methods 0.000 description 3
- 230000003287 optical effect Effects 0.000 description 3
- 230000033764 rhythmic process Effects 0.000 description 3
- 238000013461 design Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 239000000835 fiber Substances 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 208000024172 Cardiovascular disease Diseases 0.000 description 1
- 230000004913 activation Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000008280 blood Substances 0.000 description 1
- 210000004369 blood Anatomy 0.000 description 1
- 210000004204 blood vessel Anatomy 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 201000010099 disease Diseases 0.000 description 1
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 1
- 210000000056 organ Anatomy 0.000 description 1
- 238000005086 pumping Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/02—Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
- A61B5/024—Detecting, measuring or recording pulse rate or heart rate
- A61B5/02416—Detecting, measuring or recording pulse rate or heart rate using photoplethysmograph signals, e.g. generated by infrared radiation
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/24—Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
- A61B5/316—Modalities, i.e. specific diagnostic methods
- A61B5/318—Heart-related electrical modalities, e.g. electrocardiography [ECG]
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7203—Signal processing specially adapted for physiological signals or for diagnostic purposes for noise prevention, reduction or removal
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7235—Details of waveform analysis
- A61B5/725—Details of waveform analysis using specific filters therefor, e.g. Kalman or adaptive filters
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7235—Details of waveform analysis
- A61B5/7264—Classification of physiological signals or data, e.g. using neural networks, statistical classifiers, expert systems or fuzzy systems
- A61B5/7267—Classification of physiological signals or data, e.g. using neural networks, statistical classifiers, expert systems or fuzzy systems involving training the classification device
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/08—Learning methods
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Biophysics (AREA)
- Artificial Intelligence (AREA)
- General Health & Medical Sciences (AREA)
- Molecular Biology (AREA)
- Biomedical Technology (AREA)
- Animal Behavior & Ethology (AREA)
- Veterinary Medicine (AREA)
- Pathology (AREA)
- Heart & Thoracic Surgery (AREA)
- Medical Informatics (AREA)
- Public Health (AREA)
- Surgery (AREA)
- Physiology (AREA)
- Signal Processing (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Psychiatry (AREA)
- Cardiology (AREA)
- Mathematical Physics (AREA)
- Evolutionary Computation (AREA)
- Theoretical Computer Science (AREA)
- Computational Linguistics (AREA)
- General Physics & Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Software Systems (AREA)
- Computing Systems (AREA)
- Data Mining & Analysis (AREA)
- Fuzzy Systems (AREA)
- Character Discrimination (AREA)
Abstract
本发明公开了一种PPG与ECG自动转换智能算法、存储介质和计算机系统,其中算法为:S1,从开源的生理信号库中获取同步的PPG和ECG及其他生理信号,并将不同信号进行重命名后单独放置于文件夹中;S2,将步骤1获取的数据进行人工标注和划分,并进行裁剪;S3,对步骤2裁剪后的数据进行预处理;S4,分别构建并训练信号类型识别模型以及信号转换模型;S5,将步骤4训练完成的模型进行级联整合得到PPG与ECG自动转换的人工智能算法;S6,通过交互界面获得用户输入的信号数据;S7,将步骤6获取的信号数据作为步骤5中得到的PPG与ECG自动转换的人工智能算法的输入,从而输出用户所需PPG与ECG信号转换后的结果。本发明实现了PPG信号与ECG信号之间的互相转换。
Description
技术领域
本发明涉及智能医疗领域,特别涉及一种PPG与ECG自动转换智能算法、存储介质和计算机系统。
背景技术
心脏是人体最重要的器官之一,它为全身血液循环提供源源不断的动力。心脏跳动的节律是由神经电信号控制的,神经电信号的激动会进一步引发心脏的机械活动。通过神经电信号及其控制的机械活动就能够实现“泵血”的功能。心电图(electrocardiogram,ECG)就是一种记录神经电信号发生时电位变化情况的技术。当心血管疾病发生时,控制心脏节律的神经电信号会有较明显的变化。因此,通过观察ECG中心脏的节律就能大致判断疾病的类型。随着心脏机械活动的执行,人体血管内的血液流动情况会发生变化。这使得通过光信号和传感器就能记录到血液流速变化。光电容积脉搏波描记图(photoplethysmopraphy,PPG)就是这种技术。ECG和PPG分别是基于心脏电生理特性和机械活动所获得的,是心脏状态的不同表现形式。因此两者之间有着紧密的联系。目前,已有研究者对两者间的转换进行了研究,但是难以忽略的是现有的转换方法存在有一些缺陷:
(1)这些方法都不能实现PPG与ECG的互相转换。这使用户在采集到PPG时可以进行转换,而当用户采集到ECG时需要使用其他方法来转换。方法之间的差异可能会对分析结果造成不可忽略的影响。
(2)现有的部分方法是通过提取信号参数来实现信号的转换,这种方式极大的限制了转换方法的使用范围。
(3)从信号语义的角度来看,PPG信号与ECG信号存在于不同的语义空间中。现有的方法难以从空间域的角度实现两者间关系的等效映射。
发明内容
为了解决现有问题,本发明提供了一种PPG与ECG自动转换智能算法、存储介质和计算机系统,具体方案如下:
一种PPG与ECG自动转换的人工智能算法,包括以下步骤:
S1,从开源的生理信号库中获取同步的PPG信号和ECG信号及其他生理信号,并根据需求将不同信号数据进行重命名后单独放置于文件夹中;
S2,将步骤1获取的数据进行人工标注和划分,并对划分好的数据进行裁剪;
S3,对步骤2裁剪后的数据进行预处理;
S4,分别构建并训练信号类型识别模型以及信号转换模型;
S5,将步骤4训练完成的模型进行级联整合得到PPG与ECG自动转换的人工智能算法;
S6,通过交互界面获得用户输入的信号数据;
S7,将步骤6获取的信号数据作为步骤5中得到的PPG与ECG自动转换的人工智能算法的输入,从而输出用户所需PPG与ECG信号转换后的结果。
优选的,步骤2具体包括以下步骤:
S2.1,将步骤1所述文件夹中的信号数据进行人工标注,共分为3类,分别是PPG、ECG以及其他生理信号;
S2.2,遍历步骤2.1标注好的信号,并按比例随机划分为信号类型识别数据集和信号转换数据集,所述信号类型识别数据集又进一步按比例划分为训练集一、验证集一和测试集一,所述信号转换数据集进一步按比例划分为训练集二、验证集二和测试集二;
S2.3,对所有信号进行长度相同的分段,并按照预设步长进行裁剪;其中信号总长度为个体采样时间*采样频率,分段长度为预设时间*采样频率,当分段到最后一段,且长度小于所述分段长度时,将最后一段舍弃。
优选的,步骤3中的所述预处理的步骤包括对PPG进行预处理的步骤、对ECG和其他信号进行预处理的步骤;
其中,对PPG进行预处理包括以下步骤:
SA3.1,利用两个低通滤波器将PPG信号的基线拉平和均值归零;
SA3.2,结合步骤1获取信号样本的采样频率对滤波后的信号进行重采样至固定频率;
SA3.3,使用归一化函数对PPG信号进行归一化;
对ECG和其他信号进行预处理包括以下步骤:
SB3.1,利用带通滤波和陷波滤波去除ECG和其他信号中的噪音干扰;
SB3.2,结合步骤1获取信号样本的采样频率对滤波后的信号进行重采样至固定频率;
SB3.3,使用归一化函数对ECG信号和其他信号进行归一化。
优选的,步骤4构建并训练信号类型识别模型的步骤包括:
S4.1,构建用于信号类型识别的深度神经网络,输入为PPG、ECG以及其他信号中的一个信号,输出为输入信号分别为PPG信号、ECG信号以及其他信号的概率;
S4.2,初始化用于信号类型识别的深度神经网络参数并定义损失函数和优化算法;
S4.3,训练并验证用于信号类型识别的深度神经网络。
优选的,步骤4.2具体包括以下步骤:
S4.2.1,采用均值为0、方差为1的正态分布随机初始化用于信号类型识别的深度神经网络;
S4.2.2,定义损失函数来衡量真实类别和预测结果之间的差异;
S4.2.3,定义Adam优化算法为用于信号类型识别的深度神经网络的优化算法;
步骤4.3具体包括以下步骤:
S4.3.1,在训练集上通过前向传播算法来实现对信号特征的提取以预测信号的类别;
S4.3.2,通过损失函数计算预测结果和真实结果之间的损失值;
S4.3.3,结合基于链式求导法则的反向传播算法和Adam优化算法来对用于信号类型识别的深度神经网络中的权重参数进行梯度更新;
S4.3.4,在验证集上通过前向传播算法来实现对信号特征的提取以预测信号的类别;
S4.3.5,通过损失函数计算预测结果和真实结果之间的损失值;
S4.3.6,根据验证集上的损失值保存用于信号类型识别的深度神经网络;
S4.3.7,重复步骤4.3.1-步骤4.3.6,直至训练集上用于信号类型识别的深度神经网络的损失值不再降低,最终保存的深度神经网络就是用于信号类型识别的深度神经网络。
优选的,步骤4中所述信号转换模型包括PPG转ECG模型,以及ECG转PPG模型;
其中,构建并训练所述PPG转ECG模型包括以下步骤:
S4.1’,构建用于PPG转ECG的深度神经网络,输入为PPG信号,输出为由PPG生成的ECG信号;
S4.2’,初始化用于PPG转ECG的深度神经网络参数并定义损失函数和优化算法;具体包括以下步骤:
S4.2.1’,采用均值为0、方差为1的正态分布随机初始化深度神经网络;
S4.2.2’,定义损失函数来衡量原始ECG和生成ECG之间的差异;
S4.2.3’,定义Adam优化算法为用于PPG转ECG的深度神经网络的优化算法;
S4.3’,训练并验证用于PPG转ECG的深度神经网络;具体包括以下步骤:
S4.3.1’,在训练集上通过前向传播算法来实现对信号特征的提取以生成目标信号;
S4.3.2’,通过损失函数计算生成结果和真实结果之间的损失值;
S4.3.3’,结合基于链式求导法则的反向传播算法和Adam优化算法来对用于PPG转ECG的深度神经网络中的权重参数进行梯度更新;
S4.3.4’,在验证集上通过前向传播算法来实现对信号特征的提取以生成目标信号;
S4.3.5’,通过损失函数计算生成结果和真实结果之间的损失值;
S4.3.6’,根据验证集上的损失值保存用于PPG转ECG的深度神经网络;
S4.3.7’,重复步骤4.3.1’-步骤4.3.6’,直至训练集上用于PPG转ECG的深度神经网络的损失值不再降低,最终保存的深度神经网络就是用于PPG转ECG的深度神经网络;
其中,构建并训练所述ECG转PPG模型,只需将输入与输出对调位置后,再重复步骤4.1’-步骤4.3’的所有过程即可得到所述ECG转PPG模型。
优选的,步骤7中的具体步骤包括:
S7.1,将步骤6获取的信号数据作为输入数据,并通过交互界面获得用户输入的信号、频率和转换结果的保存位置;
S7.2,将所述输入数据进行裁剪;即对所述输入数据进行长度相同的分段,并按照预设步长进行裁剪;其中信号总长度为个体采样时间*采样频率,分段长度为预设时间*采样频率,当分段到最后一段,且长度小于所述分段长度时,将最后一段舍弃;
S7.3,将裁剪后的数据进行预处理;
其中,对PPG进行预处理包括以下步骤:
SA7.3.1,利用两个低通滤波器将PPG信号的基线拉平和均值归零;
SA7.3.2,结合用户输入的采样频率对滤波后的信号进行重采样至固定频率;
SA7.3.3,使用归一化函数对PPG信号进行归一化;
对ECG和其他信号进行预处理包括以下步骤:
SB7.3.1,利用带通滤波和陷波滤波去除ECG和其他信号中的噪音干扰;
SB7.3.2,结合用户输入的采样频率对滤波后的信号进行重采样至固定频率;
SB7.3.3,使用归一化函数对ECG信号和其他信号进行归一化;
S7.4,对预处理后的信号数据输入至训练好后的所述信号类型识别模型中进行信号类型识别,若信号为PPG或ECG,则进行下一步操作,若信号不为PPG或ECG,则返回给用户一个提示;
S7.5,利用已训练好的信号转换模型对PPG、ECG信号进行转换;
S7.6,将步骤7.5转换后的结果保存至用户输入的保存位置中,并随机选取一段转换结果呈现给用户。
本发明还揭示了一种计算机可读存储介质,介质上存有计算机程序,计算机程序运行后,执行上述的PPG与ECG自动转换的人工智能算法。
本发明还揭示了一种计算机系统,包括处理器、存储介质,存储介质上存有计算机程序,处理器从存储介质上读取并运行计算机程序以执行上述的PPG与ECG自动转换的人工智能算法。
本发明的有益效果在于:
本发明从空间域的角度实现了PPG信号与ECG信号两者间关系的等效映射,实现了PPG信号与ECG信号之间的互相转换。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作一简单地介绍,显而易见地,下面描述中的附图是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其它的附图。
图1为本发明的算法流程图;
图2为本发明的具体实施例的流程图;
图3为本发明的ECG转PPG的应用图例;
图4为本发明的PPG转ECG的应用图例。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地说明,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其它实施例,都属于本发明保护的范围。
如图1和图2,一种PPG与ECG自动转换的人工智能算法,包括以下步骤:
S1,从开源的生理信号库中获取同步的PPG信号和ECG信号及其他生理信号,并根据需求将不同信号数据进行重命名后单独放置于文件夹中;
S2,将步骤1获取的数据进行人工标注和划分,并对划分好的数据进行裁剪;
具体包括以下步骤:
S2.1,将步骤1所述文件夹中的信号数据进行人工标注,共分为3类,分别是PPG、ECG以及其他生理信号;所有的标注类别被添加文件夹的第一个字符,假设PPG信号、ECG信号以及其他信号所在文件夹第一个字符分别为0、1、2,则所有PPG信号标注0,所有ECG信号标注1,所有其他信号标注2。
S2.2,遍历步骤2.1标注好的信号,并按比例随机划分为信号类型识别数据集和信号转换数据集,所述信号类型识别数据集又进一步按比例,具体优选6:2:2,划分为训练集一、验证集一和测试集一,所述信号转换数据集进一步按比例,具体优选6:2:2,划分为训练集二、验证集二和测试集二;
S2.3,对所有信号进行长度相同的分段,并按照预设步长进行裁剪;其中信号总长度为个体采样时间*采样频率,分段长度为预设时间*采样频率,当分段到最后一段,且长度小于所述分段长度时,将最后一段舍弃。
假设预设时间为10s,预设步长为15s,采样频率为500HZ,个体采样时间为107s,则此时信号总长度为53500,分段长度为5000,也就是按照预设步长15s裁剪,从信号总长的起点取0-10s的一段,10s-15s的一段舍弃,再取15s-25s的一段,25s-30s的一段舍弃。依次类推,也就是每15s裁剪一段,取该15s内前10s长度的信号,舍弃后5s长度的信号,一共可以裁剪出7段长度为5000的信号,最后剩余2s的信号长度不足一个分段长度,舍弃。
S3,对步骤2裁剪后的数据进行预处理;
所述预处理的步骤包括对PPG进行预处理的步骤、对ECG和其他信号进行预处理的步骤;
其中,对PPG进行预处理包括以下步骤:
SA3.1,利用两个低通滤波器将PPG信号的基线拉平和均值归零;
SA3.2,结合步骤1获取信号样本的采样频率对滤波后的信号进行重采样至固定频率;也就是,不论采样频率是多少,该步骤都将信号进行重采样至固定的一个频率;
SA3.3,使用归一化函数对PPG信号进行归一化;
对ECG和其他信号进行预处理包括以下步骤:
SB3.1,利用带通滤波和陷波滤波去除ECG和其他信号中的噪音干扰;
SB3.2,结合步骤1获取信号样本的采样频率对滤波后的信号进行重采样至固定频率;
SB3.3,使用归一化函数对ECG信号和其他信号进行归一化。
S4,分别构建并训练信号类型识别模型以及信号转换模型;
其中,构建并训练信号类型识别模型的步骤包括:
S4.1,构建用于信号类型识别的深度神经网络,输入为PPG、ECG以及其他信号中的一个信号,输出为输入信号分别为PPG信号、ECG信号以及其他信号的概率;
S4.2,初始化用于信号类型识别的深度神经网络参数并定义损失函数和优化算法;具体包括以下步骤:
S4.2.1,采用均值为0、方差为1的正态分布随机初始化用于信号类型识别的深度神经网络;
S4.2.2,定义损失函数来衡量真实类别和预测结果之间的差异;
S4.2.3,定义Adam优化算法为用于信号类型识别的深度神经网络的优化算法;
S4.3,训练并验证用于信号类型识别的深度神经网络。具体包括以下步骤:
S4.3.1,在训练集上通过前向传播算法来实现对信号特征的提取以预测信号的类别;
S4.3.2,通过损失函数计算预测结果和真实结果之间的损失值;
S4.3.3,结合基于链式求导法则的反向传播算法和Adam优化算法来对用于信号类型识别的深度神经网络中的权重参数进行梯度更新;
S4.3.4,在验证集上通过前向传播算法来实现对信号特征的提取以预测信号的类别;
S4.3.5,通过损失函数计算预测结果和真实结果之间的损失值;
S4.3.6,根据验证集上的损失值保存用于信号类型识别的深度神经网络;
S4.3.7,重复步骤4.3.1-步骤4.3.6,直至训练集上用于信号类型识别的深度神经网络的损失值不再降低,最终保存的深度神经网络就是用于信号类型识别的深度神经网络。
其中,构建并训练所述PPG转ECG模型包括以下步骤:
S4.1’,构建用于PPG转ECG的深度神经网络,输入为PPG信号,输出为由PPG生成的ECG信号;
S4.2’,初始化用于PPG转ECG的深度神经网络参数并定义损失函数和优化算法;具体包括以下步骤:
S4.2.1’,采用均值为0、方差为1的正态分布随机初始化深度神经网络;
S4.2.2’,定义损失函数来衡量原始ECG和生成ECG之间的差异;
S4.2.3’,定义Adam优化算法为用于PPG转ECG的深度神经网络的优化算法;
S4.3’,训练并验证用于PPG转ECG的深度神经网络;具体包括以下步骤:
S4.3.1’,在训练集上通过前向传播算法来实现对信号特征的提取以生成目标信号;
S4.3.2’,通过损失函数计算生成结果和真实结果之间的损失值;
S4.3.3’,结合基于链式求导法则的反向传播算法和Adam优化算法来对用于PPG转ECG的深度神经网络中的权重参数进行梯度更新;
S4.3.4’,在验证集上通过前向传播算法来实现对信号特征的提取以生成目标信号;
S4.3.5’,通过损失函数计算生成结果和真实结果之间的损失值;
S4.3.6’,根据验证集上的损失值保存用于PPG转ECG的深度神经网络;
S4.3.7’,重复步骤4.3.1’-步骤4.3.6’,直至训练集上用于PPG转ECG的深度神经网络的损失值不再降低,最终保存的深度神经网络就是用于PPG转ECG的深度神经网络;
其中,构建并训练所述ECG转PPG模型,只需将输入与输出对调位置后,再重复步骤4.1’-步骤4.3’的所有过程即可得到所述ECG转PPG模型。
S5,将步骤4训练完成的模型进行级联整合得到PPG与ECG自动转换的人工智能算法;
S6,通过交互界面获得用户输入的信号数据;
S7,将步骤6获取的信号数据作为步骤5中得到的PPG与ECG自动转换的人工智能算法的输入,从而输出用户所需PPG与ECG信号转换后的结果。具体步骤包括:
S7.1,将步骤6获取的信号数据作为输入数据,并通过交互界面获得用户输入的信号、频率和转换结果的保存位置;
S7.2,将所述输入数据进行裁剪;即对所述输入数据进行长度相同的分段,并按照预设步长进行裁剪;其中信号总长度为个体采样时间*采样频率,分段长度为预设时间*采样频率,当分段到最后一段,且长度小于所述分段长度时,将最后一段舍弃;
S7.3,将裁剪后的数据进行预处理;
其中,对PPG进行预处理包括以下步骤:
SA7.3.1,利用两个低通滤波器将PPG信号的基线拉平和均值归零;
SA7.3.2,结合用户输入的采样频率对滤波后的信号进行重采样至固定频率;
SA7.3.3,使用归一化函数对PPG信号进行归一化;
对ECG和其他信号进行预处理包括以下步骤:
SB7.3.1,利用带通滤波和陷波滤波去除ECG和其他信号中的噪音干扰;
SB7.3.2,结合用户输入的采样频率对滤波后的信号进行重采样至固定频率;
SB7.3.3,使用归一化函数对ECG信号和其他信号进行归一化;
S7.4,对预处理后的信号数据输入至训练好后的所述信号类型识别模型中进行信号类型识别,若信号为PPG或ECG,则进行下一步操作,若信号不为PPG或ECG,则返回给用户一个提示:该信号不是PPG信号或ECG信号,请重新输入;
S7.5,利用已训练好的信号转换模型对PPG、ECG信号进行转换;
S7.6,将步骤7.5转换后的结果保存至用户输入的保存位置中,并随机选取一段转换结果呈现给用户,如图3和图4。
本发明从空间域的角度实现了PPG信号与ECG信号两者间关系的等效映射,实现了PPG信号与ECG信号之间的互相转换。
本发明还揭示了一种计算机可读存储介质,介质上存有计算机程序,计算机程序运行后,执行上述的PPG与ECG自动转换的人工智能算法。
本发明还揭示了一种计算机系统,包括处理器、存储介质,存储介质上存有计算机程序,处理器从存储介质上读取并运行计算机程序以执行上述的PPG与ECG自动转换的人工智能算法。
本领域技术人员将进一步领会,结合本文中所公开的实施例来描述的各种解说性逻辑板块、模块、电路、和算法步骤可实现为电子硬件、计算机软件、或这两者的组合。为清楚地解说硬件与软件的这一可互换性,各种解说性组件、框、模块、电路、和步骤在上面是以其功能性的形式作一般化描述的。此类功能性是被实现为硬件还是软件取决于具体应用和施加于整体系统的设计约束。技术人员对于每种特定应用可用不同的方式来实现所描述的功能性,但这样的实现决策不应被解读成导致脱离了本发明的范围。
结合本文所公开的实施例描述的各种解说性逻辑板块、模块、和电路可用通用处理器、数字信号处理器(DSP)、专用集成电路(ASIC)、现场可编程门阵列(FPGA)或其它可编程逻辑器件、分立的门或晶体管逻辑、分立的硬件组件、或其设计成执行本文所描述功能的任何组合来实现或执行。通用处理器可以是微处理器,但在替换方案中,该处理器可以是任何常规的处理器、控制器、微控制器、或状态机。处理器还可以被实现为计算设备的组合,例如DSP与微处理器的组合、多个微处理器、与DSP核心协作的一个或多个微处理器、或任何其他此类配置。
结合本文中公开的实施例描述的方法或算法的步骤可直接在硬件中、在由处理器执行的软件模块中或在这两者的组合中体现。软件模块可驻留在RAM存储器、闪存、ROM存储器、EPROM存储器、EEPROM存储器、寄存器、硬盘、可移动盘、CD-ROM、或本领域中所知的任何其他形式的存储介质中。示例性存储介质耦合到处理器以使得该处理器能从/向该存储介质读取和写入信息。在替换方案中,存储介质可以被整合到处理器。处理器和存储介质可驻留在ASIC中。ASIC可驻留在用户终端中。在替换方案中,处理器和存储介质可作为分立组件驻留在用户终端中。
在一个或多个示例性实施例中,所描述的功能可在硬件、软件、固件或其任何组合中实现。如果在软件中实现为计算机程序产品,则各功能可以作为一条或更多条指令或代码存储在计算机可读介质上或藉其进行传送。计算机可读介质包括计算机存储介质和通信介质两者,其包括促成计算机程序从一地向另一地转移的任何介质。存储介质可以是能被计算机访问的任何可用介质。作为示例而非限定,这样的计算机可读介质可包括RAM、ROM、EEPROM、CD-ROM或其它光盘存储、磁盘存储或其它磁存储设备、或能被用来携带或存储指令或数据结构形式的合意程序代码且能被计算机访问的任何其它介质。任何连接也被正当地称为计算机可读介质。例如,如果软件是使用同轴电缆、光纤电缆、双绞线、数字订户线(DSL)、或诸如红外、无线电、以及微波之类的无线技术从web网站、服务器、或其它远程源传送而来,则该同轴电缆、光纤电缆、双绞线、DSL、或诸如红外、无线电、以及微波之类的无线技术就被包括在介质的定义之中。如本文中所使用的盘(disk)和碟(disc)包括压缩碟(CD)、激光碟、光碟、数字多用碟(DVD)、软盘和蓝光碟,其中盘(disk)往往以磁的方式再现数据,而碟(disc)用激光以光学方式再现数据。上述的组合也应被包括在计算机可读介质的范围内。
提供对本公开的先前描述是为使得本领域任何技术人员皆能够制作或使用本公开。对本公开的各种修改对本领域技术人员来说都将是显而易见的,且本文中所定义的普适原理可被应用到其他变体而不会脱离本公开的精神或范围。由此,本公开并非旨在被限定于本文中所描述的示例和设计,而是应被授予与本文中所公开的原理和新颖性特征相一致的最广范围。
尽管参照前述实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。
Claims (9)
1.一种PPG与ECG自动转换的人工智能算法,其特征在于,包括以下步骤:
S1,从开源的生理信号库中获取同步的PPG信号和ECG信号及其他生理信号,并根据需求将不同信号数据进行重命名后单独放置于文件夹中;
S2,将步骤1获取的数据进行人工标注和划分,并对划分好的数据进行裁剪;
S3,对步骤2裁剪后的数据进行预处理;
S4,分别构建并训练信号类型识别模型以及信号转换模型;
S5,将步骤4训练完成的模型进行级联整合得到PPG与ECG自动转换的人工智能算法;
S6,通过交互界面获得用户输入的信号数据;
S7,将步骤6获取的信号数据作为步骤5中得到的PPG与ECG自动转换的人工智能算法的输入,从而输出用户所需PPG与ECG信号转换后的结果。
2.根据权利要求1所述的PPG与ECG自动转换的人工智能算法,其特征在于,步骤2具体包括以下步骤:
S2.1,将步骤1所述文件夹中的信号数据进行人工标注,共分为3类,分别是PPG、ECG以及其他生理信号;
S2.2,遍历步骤2.1标注好的信号,并按比例随机划分为信号类型识别数据集和信号转换数据集,所述信号类型识别数据集又进一步按比例划分为训练集一、验证集一和测试集一,所述信号转换数据集进一步按比例划分为训练集二、验证集二和测试集二;
S2.3,对所有信号进行长度相同的分段,并按照预设步长进行裁剪;其中信号总长度为个体采样时间*采样频率,分段长度为预设时间*采样频率,当分段到最后一段,且长度小于所述分段长度时,将最后一段舍弃。
3.根据权利要求1所述的PPG与ECG自动转换的人工智能算法,其特征在于,步骤3中的所述预处理的步骤包括对PPG进行预处理的步骤、对ECG和其他信号进行预处理的步骤;
其中,对PPG进行预处理包括以下步骤:
SA3.1,利用两个低通滤波器将PPG信号的基线拉平和均值归零;
SA3.2,结合步骤1获取信号样本的采样频率对滤波后的信号进行重采样至固定频率;
SA3.3,使用归一化函数对PPG信号进行归一化;
对ECG和其他信号进行预处理包括以下步骤:
SB3.1,利用带通滤波和陷波滤波去除ECG和其他信号中的噪音干扰;
SB3.2,结合步骤1获取信号样本的采样频率对滤波后的信号进行重采样至固定频率;
SB3.3,使用归一化函数对ECG信号和其他信号进行归一化。
4.根据权利要求1所述的PPG与ECG自动转换的人工智能算法,其特征在于,步骤4构建并训练信号类型识别模型的步骤包括:
S4.1,构建用于信号类型识别的深度神经网络,输入为PPG、ECG以及其他信号中的一个信号,输出为输入信号分别为PPG信号、ECG信号以及其他信号的概率;
S4.2,初始化用于信号类型识别的深度神经网络参数并定义损失函数和优化算法;
S4.3,训练并验证用于信号类型识别的深度神经网络。
5.根据权利要求4所述的PPG与ECG自动转换的人工智能算法,其特征在于,步骤4.2具体包括以下步骤:
S4.2.1,采用均值为0、方差为1的正态分布随机初始化用于信号类型识别的深度神经网络;
S4.2.2,定义损失函数来衡量真实类别和预测结果之间的差异;
S4.2.3,定义Adam优化算法为用于信号类型识别的深度神经网络的优化算法;
步骤4.3具体包括以下步骤:
S4.3.1,在训练集上通过前向传播算法来实现对信号特征的提取以预测信号的类别;
S4.3.2,通过损失函数计算预测结果和真实结果之间的损失值;
S4.3.3,结合基于链式求导法则的反向传播算法和Adam优化算法来对用于信号类型识别的深度神经网络中的权重参数进行梯度更新;
S4.3.4,在验证集上通过前向传播算法来实现对信号特征的提取以预测信号的类别;
S4.3.5,通过损失函数计算预测结果和真实结果之间的损失值;
S4.3.6,根据验证集上的损失值保存用于信号类型识别的深度神经网络;
S4.3.7,重复步骤4.3.1-步骤4.3.6,直至训练集上用于信号类型识别的深度神经网络的损失值不再降低,最终保存的深度神经网络就是用于信号类型识别的深度神经网络。
6.根据权利要求1所述的PPG与ECG自动转换的人工智能算法,其特征在于,步骤4中所述信号转换模型包括PPG转ECG模型,以及ECG转PPG模型;
其中,构建并训练所述PPG转ECG模型包括以下步骤:
S4.1’,构建用于PPG转ECG的深度神经网络,输入为PPG信号,输出为由PPG生成的ECG信号;
S4.2’,初始化用于PPG转ECG的深度神经网络参数并定义损失函数和优化算法;具体包括以下步骤:
S4.2.1’,采用均值为0、方差为1的正态分布随机初始化深度神经网络;
S4.2.2’,定义损失函数来衡量原始ECG和生成ECG之间的差异;
S4.2.3’,定义Adam优化算法为用于PPG转ECG的深度神经网络的优化算法;
S4.3’,训练并验证用于PPG转ECG的深度神经网络;具体包括以下步骤:
S4.3.1’,在训练集上通过前向传播算法来实现对信号特征的提取以生成目标信号;
S4.3.2’,通过损失函数计算生成结果和真实结果之间的损失值;
S4.3.3’,结合基于链式求导法则的反向传播算法和Adam优化算法来对用于PPG转ECG的深度神经网络中的权重参数进行梯度更新;
S4.3.4’,在验证集上通过前向传播算法来实现对信号特征的提取以生成目标信号;
S4.3.5’,通过损失函数计算生成结果和真实结果之间的损失值;
S4.3.6’,根据验证集上的损失值保存用于PPG转ECG的深度神经网络;
S4.3.7’,重复步骤4.3.1’-步骤4.3.6’,直至训练集上用于PPG转ECG的深度神经网络的损失值不再降低,最终保存的深度神经网络就是用于PPG转ECG的深度神经网络;
其中,构建并训练所述ECG转PPG模型,只需将输入与输出对调位置后,再重复步骤4.1’-步骤4.3’的所有过程即可得到所述ECG转PPG模型。
7.根据权利要求1所述的PPG与ECG自动转换的人工智能算法,其特征在于,步骤7中的具体步骤包括:
S7.1,将步骤6获取的信号数据作为输入数据,并通过交互界面获得用户输入的信号、频率和转换结果的保存位置;
S7.2,将所述输入数据进行裁剪;即对所述输入数据进行长度相同的分段,并按照预设步长进行裁剪;其中信号总长度为个体采样时间*采样频率,分段长度为预设时间*采样频率,当分段到最后一段,且长度小于所述分段长度时,将最后一段舍弃;
S7.3,将裁剪后的数据进行预处理;
其中,对PPG进行预处理包括以下步骤:
SA7.3.1,利用两个低通滤波器将PPG信号的基线拉平和均值归零;
SA7.3.2,结合用户输入的采样频率对滤波后的信号进行重采样至固定频率;
SA7.3.3,使用归一化函数对PPG信号进行归一化;
对ECG和其他信号进行预处理包括以下步骤:
SB7.3.1,利用带通滤波和陷波滤波去除ECG和其他信号中的噪音干扰;
SB7.3.2,结合用户输入的采样频率对滤波后的信号进行重采样至固定频率;
SB7.3.3,使用归一化函数对ECG信号和其他信号进行归一化;
S7.4,对预处理后的信号数据输入至训练好后的所述信号类型识别模型中进行信号类型识别,若信号为PPG或ECG,则进行下一步操作,若信号不为PPG或ECG,则返回给用户一个提示;
S7.5,利用已训练好的信号转换模型对PPG、ECG信号进行转换;
S7.6,将步骤7.5转换后的结果保存至用户输入的保存位置中,并随机选取一段转换结果呈现给用户。
8.一种计算机可读存储介质,其特征在于:介质上存有计算机程序,计算机程序运行后,执行如权利要求1至7中任一项所述的PPG与ECG自动转换的人工智能算法。
9.一种计算机系统,其特征在于:包括处理器、存储介质,存储介质上存有计算机程序,处理器从存储介质上读取并运行计算机程序以执行如权利要求1至7中任一项所述的PPG与ECG自动转换的人工智能算法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210688131.8A CN115270848B (zh) | 2022-06-17 | 2022-06-17 | 一种ppg与ecg自动转换智能算法、存储介质和计算机系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210688131.8A CN115270848B (zh) | 2022-06-17 | 2022-06-17 | 一种ppg与ecg自动转换智能算法、存储介质和计算机系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN115270848A true CN115270848A (zh) | 2022-11-01 |
CN115270848B CN115270848B (zh) | 2023-09-29 |
Family
ID=83761661
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210688131.8A Active CN115270848B (zh) | 2022-06-17 | 2022-06-17 | 一种ppg与ecg自动转换智能算法、存储介质和计算机系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115270848B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116584952A (zh) * | 2023-07-18 | 2023-08-15 | 牛津大学(苏州)科技有限公司 | 心电图异常监测方法、装置、设备和存储介质 |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20200196897A1 (en) * | 2018-12-20 | 2020-06-25 | Imec Vzw | Method for generating a model for generating a synthetic ecg and a method and system for analysis of heart activity |
CN111839495A (zh) * | 2020-07-30 | 2020-10-30 | 深圳前海微众银行股份有限公司 | 检测方法、设备和存储介质 |
WO2021087985A1 (zh) * | 2019-11-08 | 2021-05-14 | 深圳市欢太科技有限公司 | 模型训练方法、装置、存储介质及电子设备 |
CN113057648A (zh) * | 2021-03-22 | 2021-07-02 | 山西三友和智慧信息技术股份有限公司 | 一种基于复合lstm结构的ecg信号分类方法 |
WO2021151345A1 (zh) * | 2020-07-09 | 2021-08-05 | 平安科技(深圳)有限公司 | 识别模型的参数获取方法、装置、电子设备及存储介质 |
WO2021208490A1 (zh) * | 2020-04-13 | 2021-10-21 | 中国科学院深圳先进技术研究院 | 一种基于深度神经网络的血压测量方法及装置 |
CN114119968A (zh) * | 2021-12-16 | 2022-03-01 | 安徽心之声医疗科技有限公司 | 基于深度神经网络的纸质心电图信号区域提取方法及系统 |
CN114391824A (zh) * | 2022-01-18 | 2022-04-26 | 安徽大学 | 一种基于ppg信号的心率估计方法、电子设备及存储介质 |
WO2022106146A1 (en) * | 2020-11-17 | 2022-05-27 | Siemens Aktiengesellschaft | Transforming a trained artificial intelligence model into a trustworthy artificial intelligence model |
-
2022
- 2022-06-17 CN CN202210688131.8A patent/CN115270848B/zh active Active
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20200196897A1 (en) * | 2018-12-20 | 2020-06-25 | Imec Vzw | Method for generating a model for generating a synthetic ecg and a method and system for analysis of heart activity |
WO2021087985A1 (zh) * | 2019-11-08 | 2021-05-14 | 深圳市欢太科技有限公司 | 模型训练方法、装置、存储介质及电子设备 |
WO2021208490A1 (zh) * | 2020-04-13 | 2021-10-21 | 中国科学院深圳先进技术研究院 | 一种基于深度神经网络的血压测量方法及装置 |
WO2021151345A1 (zh) * | 2020-07-09 | 2021-08-05 | 平安科技(深圳)有限公司 | 识别模型的参数获取方法、装置、电子设备及存储介质 |
CN111839495A (zh) * | 2020-07-30 | 2020-10-30 | 深圳前海微众银行股份有限公司 | 检测方法、设备和存储介质 |
WO2022106146A1 (en) * | 2020-11-17 | 2022-05-27 | Siemens Aktiengesellschaft | Transforming a trained artificial intelligence model into a trustworthy artificial intelligence model |
CN113057648A (zh) * | 2021-03-22 | 2021-07-02 | 山西三友和智慧信息技术股份有限公司 | 一种基于复合lstm结构的ecg信号分类方法 |
CN114119968A (zh) * | 2021-12-16 | 2022-03-01 | 安徽心之声医疗科技有限公司 | 基于深度神经网络的纸质心电图信号区域提取方法及系统 |
CN114391824A (zh) * | 2022-01-18 | 2022-04-26 | 安徽大学 | 一种基于ppg信号的心率估计方法、电子设备及存储介质 |
Non-Patent Citations (1)
Title |
---|
荣凡稳;郑伟;陈冉;高军峰;: "基于深度学习的运动心率测量系统", 电子测量与仪器学报, no. 12 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116584952A (zh) * | 2023-07-18 | 2023-08-15 | 牛津大学(苏州)科技有限公司 | 心电图异常监测方法、装置、设备和存储介质 |
CN116584952B (zh) * | 2023-07-18 | 2023-12-08 | 牛津大学(苏州)科技有限公司 | 心电图异常监测方法、装置、设备和存储介质 |
Also Published As
Publication number | Publication date |
---|---|
CN115270848B (zh) | 2023-09-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US9949714B2 (en) | Method, electronic apparatus, and computer readable medium of constructing classifier for disease detection | |
Nogueira et al. | Classifying heart sounds using images of motifs, MFCC and temporal features | |
Tang et al. | PCG classification using multidomain features and SVM classifier | |
CN109864736A (zh) | 心电信号的处理方法、装置、终端设备及介质 | |
Alqudah et al. | Classification of heart sound short records using bispectrum analysis approach images and deep learning | |
US20140276115A1 (en) | Heart sound simulator | |
Singh et al. | ECG signal feature extraction trends in methods and applications | |
Milani et al. | Abnormal heart sound classification using phonocardiography signals | |
Ismail et al. | PCG classification through spectrogram using transfer learning | |
JP2023528988A (ja) | 心拍リズムを分析するための方法およびシステム | |
CN115270848B (zh) | 一种ppg与ecg自动转换智能算法、存储介质和计算机系统 | |
CN112949639B (zh) | 一种先天性心脏病心音智能分割分类算法、装置及存储介质 | |
KR102304339B1 (ko) | 인공 신경망을 활용하여 혈관 초음파 영상에 대한 분석 정보를 제공하기 위한 방법 및 컴퓨팅 장치 | |
CN110946566A (zh) | 一种基于U-Net网络的心拍分类方法、装置、设备及存储介质 | |
Deevi et al. | HeartNetEC: a deep representation learning approach for ECG beat classification | |
Ma et al. | Application of time-frequency domain and deep learning fusion feature in non-invasive diagnosis of congenital heart disease-related pulmonary arterial hypertension | |
Li et al. | Heart sound classification based on improved mel-frequency spectral coefficients and deep residual learning | |
Edinburgh et al. | DeepClean: Self-supervised artefact rejection for intensive care waveform data using deep generative learning | |
Malik et al. | Accurate classification of heart sound signals for cardiovascular disease diagnosis by wavelet analysis and convolutional neural network: preliminary results | |
Mustafa et al. | Detection of heartbeat sounds arrhythmia using automatic spectral methods and cardiac auscultatory | |
Tanasković et al. | A new algorithm for fetal heart rate detection: Fractional order calculus approach | |
CN113598782A (zh) | 预测室性心律失常起源位置的系统、电子设备和存储介质 | |
CN110786847B (zh) | 心电信号的建库方法和分析方法 | |
CN115828168A (zh) | 一种基于隐马尔可夫心动周期的心音分割方法及其装置 | |
Priyasad et al. | Detecting heart failure through voice analysis using self-supervised mode-based memory fusion |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |