CN112386227A - 一种多尺度时间序列生理信号的因果网络分析方法 - Google Patents
一种多尺度时间序列生理信号的因果网络分析方法 Download PDFInfo
- Publication number
- CN112386227A CN112386227A CN202011239121.3A CN202011239121A CN112386227A CN 112386227 A CN112386227 A CN 112386227A CN 202011239121 A CN202011239121 A CN 202011239121A CN 112386227 A CN112386227 A CN 112386227A
- Authority
- CN
- China
- Prior art keywords
- imf
- signals
- signal
- causal
- instantaneous phase
- 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
Images
Classifications
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
- G16H40/00—ICT specially adapted for the management or administration of healthcare resources or facilities; ICT specially adapted for the management or operation of medical equipment or devices
- G16H40/60—ICT specially adapted for the management or administration of healthcare resources or facilities; ICT specially adapted for the management or operation of medical equipment or devices for the operation of medical equipment or devices
- G16H40/63—ICT specially adapted for the management or administration of healthcare resources or facilities; ICT specially adapted for the management or operation of medical equipment or devices for the operation of medical equipment or devices for local operation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N5/00—Computing arrangements using knowledge-based models
- G06N5/04—Inference or reasoning models
-
- 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
Abstract
本发明公开了一种多尺度时间序列生理信号的因果网络分析方法,包括以下步骤:S1.利用NA‑MEMD算法对待分析的生理信号u1,u2,...,um进行分解;S2.对于两个不同的生理信号ui,uj,进行因果分析,i=1,2,...m,j=1,2...,m,且i≠j,得到两个信号之间的因果关系;S3.对于u1,u2,...,um中任意两个信号,重复执行步骤S2,直至得到u1,u2,...,um中各个信号两两之间的因果关系,形成因果网络。本发明能够有效分析出生理信号的因果网络,为生理信号的应用提供了便利。
Description
技术领域
本发明涉及生理信号处理,特别是涉及一种多尺度时间序列生理信号的因果网络分析方法。
背景技术
多尺度时间序列生理信号的因果网络分析方法可在脑机接口技术、脑功能和结构机制研究、脑器交互技术(大脑与眼、心、肺、肌等人体器官构成的生理网络)中得以应用和推广。
基于Gall提出的神经科学功能分离和功能整合表明,不同的大脑功能在专门的皮质区域是存在局部局限性的,因此通过功能连接(FC)来评估大脑感知、认知和行为的科学标准一直由格兰杰因果关系分析主导。FC被用作为一种视觉工具,以显示大脑区域之间空间距离的时间相关性。
在大多数真实情况下,因果关系可能是依赖时间的、同时的和相互的。然而,现有的时间序列因果关系分析方法大多基于预测的概念,可能无法描述瞬时事件之间的相互因果关系运动。在这种情况下,格兰杰因果关系是基于以下假设:因果关系是可分离的,这在许多线性随机系统中很有用,但可能不适用于复杂的动态过程(例如与大脑相关的生理网络)。此外,还开发了收敛交叉映射(CCM)方法,以适应因果效应的不可分离性。在大脑网络的因果建模中,动态因果建模和转移熵也很普遍,但基本上基于贝叶斯预测。杨等人提出了一种通过希尔伯特-黄变换(HHT)评估因果相互作用的方法。然而,它的理论框架不能直接在大脑相关的复杂信息系统上进行推广应用,而且在先验因果相互作用方面也缺乏充分必要的实践。
协变量方法所依据的是伽利略和休谟的因果归纳,由于因果关系以某些方式从感官输入中编码:即,候选因果与结果之间的共变量可以定义为在候选原因存在时效果概率与在候选原因不存在的情况下,该概率之间的差值。Kantian的因果强度表示解释了因果信息的先验知识的存在。与通用的基于HHT的因果分解方法相比,经验模式分解(EMD)不能保证跨多个时间序列的IMF数相等,并导致模式(比例)对齐(跨通道相互依赖),从而导致多元时间序列中时间尺度的延伸。例如,在大脑-生理网络中,从信号采集设备收集到的振荡分量通常以不同的频段采样。另外,EMD方法还容易受到导致模式混合问题(单通道独立性)的噪声影响,导致不足以识别分解的固有因果分量(ICCs)。
到目前为止,多时间尺度的因果关系分析已基本应用于时间序列的ICC评估。多尺度格兰杰因果关系通过过滤和下采样步骤来解决,以获得双变量过程的再缩放表示。基于转移熵的多尺度分析主要基于傅里叶和小波变换得到的固定的基函数。但是,当应用于自然继承的非线性和非固定属性的复杂过程时,这些分解过程不会捕获不同时间尺度上的重要特征。
发明内容
本发明的目的在于克服现有技术的不足,提供一种多尺度时间序列生理信号的因果网络分析方法,能够有效分析出生理信号的因果网络,为生理信号的应用提供了便利。
本发明的目的是通过以下技术方案来实现的:一种多尺度时间序列生理信号的因果网络分析方法,包括以下步骤:
S1.输入待分析的生理信号:
u1={u1,1,u1,2,...,u1,t}
u2={u2,1,u2,2,...,u2,t}
......
um={um,1,um,2,...,um,t}
利用NA-MEMD算法对待分析的生理信号u1,u2,...,um进行分解:
其中,表示利用NA-MEMD算法对信号进行分解,NA-MEMD算法是指基于噪声辅助的多元经验模式分解算法,m表示输入的生理信号数目,并且满足m≥2,t∈N+,N+表示正整数;表示NA-MEMD算法选择的辅助噪声,互不相关且均为随机高斯噪声,选择的辅助噪声数目,n表示每个输入的生理信号分解后得到的IMF信号个数;
S2.对于两个不同的生理信号ui,uj,进行因果分析,i=1,2,...m,j=1,2...,m,且i≠j,分析过程如下:
S201.将生理信号ui分解得到的信号{IMFi,1,IMFi,2,...,IMFi,n}与生理信号uj分解得到的信号{IMFj,1,IMFj,2,...,IMFj,n}进行配对,得到n个IMF对:
(IMFi,1,IMFj,1),(IMFi,2,IMFj,2),....,(IMFi,n,IMFj,n);
每个IMF对中,两个IMF信号的时间长度相等;
S202.计算每一个IMF对的平均瞬时相位差,并将计算得到的平均瞬时相位差与预设的阈值进行比较,筛选出平均瞬时相位差小于预设阈值的IMF对,生成ICC集:
S204.分别计算ICC集中每一个IMF对的相位相关性;
S204.信号重分解:
分解后取出u′j对应的分解信号{IMF′j,1,IMF′j,2,...,IMF′j,n}
分解后取出u′i对应的分解信号{IMF′i,1,IMF′i,2,...,IMF′i,n};
获取绝对因果强度ACS:
S206.根据绝对因果强度ACS,计算如下比值:
若该比值大于1,则ui为因,uj为果;
若该比值小于1,则ui为果,uj为因;
若该比值等于1,则ui、uj无因果关系或者互为因果;
由此得到ui、uj的因果分析结果;
S3.对于u1,u2,...,um中任意两个信号,重复执行步骤S2,直至得到u1,u2,...,um中各个信号两两之间的因果关系,形成因果网络。
其中,所述步骤S202包括:
S2021.设定各个IMF对的平均瞬时相位差阈值δ1,δ2,...δn;
S2022.计算第h个IMF对(IMFi,h,IMFj,h)的平均瞬时相位差:
设mean(φi,h)表示IMFi,h的平均瞬时相位,即IMFi,h在时间长度内的瞬时相位平均值;mean(φj,h)表示IMFj,h的平均瞬时相位,即IMFj,h时间长度内的瞬时相位平均值;
则第h个IMF对(IMFi,h,IMFj,h)的平均瞬时相位差为:
|mean(φi,h)-mean(φj,h)|;
将|mean(φi,h)-mean(φj,h)|与对应阈值δh进行比较,判断是否满足:
|mean(φi,h)-mean(φj,h)|<δh;
如果满足,则将对应的(IMFi,h,IMFj,h)加入ICC集中;
若不满足,则抛弃(IMFi,h,IMFj,h);
S2023.在h=1,2,...n时,分别重复执行步骤S2022,最终得到ICC集为:
本发明的有益效果是:本发明提供了一种多尺度时间序列生理信号的因果网络分析方法,有效分析出生理信号的因果网络,相较于传统方法,本申请的分析结果能够更加完整地体现对时间序列的因果关系,为生理信号在脑机接口技术、脑功能和结构机制研究、脑器交互技术上的应用提供了技术条件。
附图说明
图1为本发明的方法流程图。
具体实施方式
下面结合附图进一步详细描述本发明的技术方案,但本发明的保护范围不局限于以下所述。
如图1所示,一种多尺度时间序列生理信号的因果网络分析方法,包括以下步骤:
S1.输入待分析的生理信号:
u1={u1,1,u1,2,...,u1,t}
u2={u2,1,u2,2,...,u2,t}
......
um={um,1,um,2,...,um,t}
利用NA-MEMD算法对待分析的生理信号u1,u2,...,um进行分解:
其中,表示利用NA-MEMD算法对信号进行分解,NA-MEMD算法是指基于噪声辅助的多元经验模式分解算法,m表示输入的生理信号数目,并且满足m≥2,t∈N+,N+表示正整数;表示NA-MEMD算法选择的辅助噪声,互不相关且均为随机高斯噪声,为选择的辅助噪声数目,n表示每个输入的生理信号分解后得到的IMF信号个数;
S2.对于两个不同的生理信号ui,uj,进行因果分析,i=1,2,...m,j=1,2...,m,且i≠j,分析过程如下:
S201.将生理信号ui分解得到的信号{IMFi,1,IMFi,2,...,IMFi,n}与生理信号uj分解得到的信号{IMFj,1,IMFj,2,...,IMFj,n}进行配对,得到n个IMF对:
(IMFi,1,IMFj,1),(IMFi,2,IMFj,2),....,(IMFi,n,IMFj,n);
每个IMF对中,两个IMF信号的时间长度相等;
S202.计算每一个IMF对的平均瞬时相位差,并将计算得到的平均瞬时相位差与预设的阈值进行比较,筛选出平均瞬时相位差小于预设阈值的IMF对,生成ICC集:
S204.分别计算ICC集中每一个IMF对的相位相关性;
S204.信号重分解:
分解后取出u′j对应的分解信号{IMF′j,1,IMF′j,2,...,IMF′j,n}
分解后取出ui′对应的分解信号{IMF′i,1,IMF′i,2,...,IMF′i,n};
获取绝对因果强度ACS:
S206.根据绝对因果强度ACS,计算如下比值:
若该比值大于1,则ui为因,uj为果;
若该比值小于1,则ui为果,uj为因;
若该比值等于1,则ui、uj无因果关系或者互为因果;
由此得到ui、uj的因果分析结果;
S3.对于u1,u2,...,um中任意两个信号,重复执行步骤S2,直至得到u1,u2,...,um中各个信号两两之间的因果关系,形成因果网络。
其中,所述步骤S202包括:
S2021.设定各个IMF对的平均瞬时相位差阈值δ1,δ2,...δn;
S2022.计算第h个IMF对(IMFi,h,IMFj,h)的平均瞬时相位差:
设mean(φi,h)表示IMFi,h的平均瞬时相位,即IMFi,h在时间长度内的瞬时相位平均值;mean(φj,h)表示IMFj,h的平均瞬时相位,即IMFj,h时间长度内的瞬时相位平均值;
则第h个IMF对(IMFi,h,IMFj,h)的平均瞬时相位差为:
|mean(φi,h)-mean(φj,h)|;
将|mean(φi,h)-mean(φj,h)|与对应阈值δh进行比较,判断是否满足:
|mean(φi,h)-mean(φj,h)|<δh;
如果满足,则将对应的(IMFi,h,IMFj,h)加入ICC集中;
若不满足,则抛弃(IMFi,h,IMFj,h);
S2023.在h=1,2,...n时,分别重复执行步骤S2022,最终得到ICC集为:
在本申请的实施例中,两个生理信号A、B之间的因果关系定义为:如果A和B中分解的IMF处于相同的时间尺度,当B的IMF从B本身删除时,若A和B中的IMF之间的瞬时相位依赖关系消失,则变量A会导致变量B,即A是B的因(A为因,B为果);若当B的IMF从B本身删除时,A和B中的IMF之间的瞬时相位依赖关系没有消失,则变量A不会导致B,即A不是B的因。
因果关系需要满足如下条件:
(1)任何因果关系都基于ICC在多个时间序列中的瞬时相位相关性;
(2)基于相位的因果本身是能够分离的。
以上虽然描述了本发明的具体实施方法,但是本领域的技术人员应当理解,这些仅是举例说明,在不背离本发明原理和实现的前提下,可以对这些实施方案做出多种变更或修改,因此,本发明的保护范围由所附权利要求书限定。
Claims (2)
1.一种多尺度时间序列生理信号的因果网络分析方法,其特征在于:包括以下步骤:
S1.输入待分析的生理信号:
u1={u1,1,u1,2,...,u1,t}
u2={u2,1,u2,2,...,u2,t}
......
um={um,1,um,2,...,um,t}
利用NA-MEMD算法对待分析的生理信号u1,u2,...,um进行分解:
......
......
其中,表示利用NA-MEMD算法对信号进行分解,NA-MEMD算法是指基于噪声辅助的多元经验模式分解算法,m表示输入的生理信号数目,并且满足m≥2,t∈N+,N+表示正整数;表示NA-MEMD算法选择的辅助噪声,互不相关且均为随机高斯噪声,为选择的辅助噪声数目,n表示每个输入的生理信号分解后得到的IMF信号个数;
S2.对于两个不同的生理信号ui,uj,进行因果分析,i=1,2,...m,j=1,2...,m,且i≠j,分析过程如下:
S201.将生理信号ui分解得到的信号{IMFi,1,IMFi,2,...,IMFi,n}与生理信号uj分解得到的信号{IMFj,1,IMFj,2,...,IMFj,n}进行配对,得到n个IMF对:
(IMFi,1,IMFj,1),(IMFi,2,IMFj,2),....,(IMFi,n,IMFj,n);
每个IMF对中,两个IMF信号的时间长度相等;
S202.计算每一个IMF对的平均瞬时相位差,并将计算得到的平均瞬时相位差与预设的阈值进行比较,筛选出平均瞬时相位差小于预设阈值的IMF对,生成ICC集:
S204.分别计算ICC集中每一个IMF对的相位相关性;
S204.信号重分解:
分解后取出u′j对应的分解信号{IMF′j,1,IMF′j,2,...,IMF′j,n}
分解后取出u′i对应的分解信号{IMF′i,1,IMF′i,2,...,IMF′i,n};
获取绝对因果强度ACS:
S206.根据绝对因果强度ACS,计算如下比值:
若该比值大于1,则ui为因,uj为果;
若该比值小于1,则ui为果,uj为因;
若该比值等于1,则ui、uj无因果关系或者互为因果;
由此得到ui、uj的因果分析结果;
S3.对于u1,u2,...,um中任意两个信号,重复执行步骤S2,直至得到u1,u2,...,um中各个信号两两之间的因果关系,形成因果网络。
2.根据权利要求1所述的一种多尺度时间序列生理信号的因果网络分析方法,其特征在于:所述步骤S202包括:
S2021.设定各个IMF对的平均瞬时相位差阈值δ1,δ2,...δn;
S2022.计算第h个IMF对(IMFi,h,IMFj,h)的平均瞬时相位差:
设mean(φi,h)表示IMFi,h的平均瞬时相位,即IMFi,h在时间长度内的瞬时相位平均值;mean(φj,h)表示IMFj,h的平均瞬时相位,即IMFj,h时间长度内的瞬时相位平均值;
则第h个IMF对(IMFi,h,IMFj,h)的平均瞬时相位差为:
|mean(φi,h)-mean(φj,h)|;
将|mean(φi,h)-mean(φj,h)|与对应阈值δh进行比较,判断是否满足:
|mean(φi,h)-mean(φj,h)|<δh;
如果满足,则将对应的(IMFi,h,IMFj,h)加入ICC集中;
若不满足,则抛弃(IMFi,h,IMFj,h);
S2023.在h=1,2,...n时,分别重复执行步骤S2022,最终得到ICC集为:
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011239121.3A CN112386227B (zh) | 2020-11-09 | 2020-11-09 | 一种多尺度时间序列生理信号的因果网络分析方法 |
US17/396,800 US20220148721A1 (en) | 2020-11-09 | 2021-08-09 | Analysis method for causal inference of physiological network in multiscale time series signals |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011239121.3A CN112386227B (zh) | 2020-11-09 | 2020-11-09 | 一种多尺度时间序列生理信号的因果网络分析方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112386227A true CN112386227A (zh) | 2021-02-23 |
CN112386227B CN112386227B (zh) | 2021-07-02 |
Family
ID=74597875
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011239121.3A Active CN112386227B (zh) | 2020-11-09 | 2020-11-09 | 一种多尺度时间序列生理信号的因果网络分析方法 |
Country Status (2)
Country | Link |
---|---|
US (1) | US20220148721A1 (zh) |
CN (1) | CN112386227B (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113143296A (zh) * | 2021-04-20 | 2021-07-23 | 上海外国语大学 | 沟通障碍的智能评估方法及评估系统 |
CN114154543A (zh) * | 2021-11-30 | 2022-03-08 | 电子科技大学 | 基于因果特征的人体上肢运动神经行为辨识模型构建方法 |
CN114460916A (zh) * | 2021-12-22 | 2022-05-10 | 浙江大学 | 一种基于因果分析的子系统波动信号分析方法 |
CN114510966A (zh) * | 2022-01-14 | 2022-05-17 | 电子科技大学 | 一种基于图神经网络的端到端大脑因果网络构建方法 |
CN116743602A (zh) * | 2023-08-03 | 2023-09-12 | 广东广宇科技发展有限公司 | 一种基于收敛交叉映射的通信线路安全分析方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109758145A (zh) * | 2018-12-15 | 2019-05-17 | 北京交通大学 | 基于脑电信号因果关系的自动睡眠分期方法 |
CN111859274A (zh) * | 2020-07-03 | 2020-10-30 | 自然资源部第一海洋研究所 | 一种用于度量复杂动力系统本征因果关系的方法 |
-
2020
- 2020-11-09 CN CN202011239121.3A patent/CN112386227B/zh active Active
-
2021
- 2021-08-09 US US17/396,800 patent/US20220148721A1/en not_active Abandoned
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109758145A (zh) * | 2018-12-15 | 2019-05-17 | 北京交通大学 | 基于脑电信号因果关系的自动睡眠分期方法 |
CN111859274A (zh) * | 2020-07-03 | 2020-10-30 | 自然资源部第一海洋研究所 | 一种用于度量复杂动力系统本征因果关系的方法 |
Non-Patent Citations (2)
Title |
---|
YANG, AC等: "Causal decomposition in the mutual causation", 《NATURE COMMUNICATIONS》 * |
杜鹏等: "非线性格兰杰因果关系在睡眠生理信号分析中的应用", 《数据采集与处理》 * |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113143296A (zh) * | 2021-04-20 | 2021-07-23 | 上海外国语大学 | 沟通障碍的智能评估方法及评估系统 |
CN114154543A (zh) * | 2021-11-30 | 2022-03-08 | 电子科技大学 | 基于因果特征的人体上肢运动神经行为辨识模型构建方法 |
CN114154543B (zh) * | 2021-11-30 | 2023-04-28 | 电子科技大学 | 基于因果特征的人体上肢运动神经行为辨识模型构建方法 |
CN114460916A (zh) * | 2021-12-22 | 2022-05-10 | 浙江大学 | 一种基于因果分析的子系统波动信号分析方法 |
CN114510966A (zh) * | 2022-01-14 | 2022-05-17 | 电子科技大学 | 一种基于图神经网络的端到端大脑因果网络构建方法 |
CN116743602A (zh) * | 2023-08-03 | 2023-09-12 | 广东广宇科技发展有限公司 | 一种基于收敛交叉映射的通信线路安全分析方法 |
CN116743602B (zh) * | 2023-08-03 | 2023-12-01 | 广东广宇科技发展有限公司 | 一种基于收敛交叉映射的通信线路安全分析方法 |
Also Published As
Publication number | Publication date |
---|---|
CN112386227B (zh) | 2021-07-02 |
US20220148721A1 (en) | 2022-05-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112386227B (zh) | 一种多尺度时间序列生理信号的因果网络分析方法 | |
Toderici et al. | Full resolution image compression with recurrent neural networks | |
CN101515200B (zh) | 基于瞬态视觉诱发脑电的目标选择方法 | |
CN113128552B (zh) | 一种基于深度可分离因果图卷积网络的脑电情绪识别方法 | |
Capizzi et al. | A wavelet based prediction of wind and solar energy for long-term simulation of integrated generation systems | |
Zhang et al. | Diffusion kernel attention network for brain disorder classification | |
JP6099032B2 (ja) | 信号処理装置、信号処理方法及びコンピュータプログラム | |
Xiang et al. | Estimation of the evolutionary spectra with application to stationarity test | |
Tabassum et al. | Comparison of filter banks of DWT in recovery of image using one dimensional signal vector | |
Samal et al. | Ensemble median empirical mode decomposition for emotion recognition using EEG signal | |
Al-Assaf et al. | Identification of fractional chaotic system parameters | |
Al Bassam et al. | Wavelet theory and application in communication and signal processing | |
Hwang et al. | Null space component analysis of one-shot single-channel source separation problem | |
Feng et al. | Singular spectral analysis‐based denoising without computing singular values via augmented Lagrange multiplier algorithm | |
Yarga et al. | Efficient spike encoding algorithms for neuromorphic speech recognition | |
Potter et al. | Competing ICA techniques in biomedical signal analysis | |
Kowalski et al. | Random models for sparse signals expansion on unions of bases with application to audio signals | |
Desai et al. | Compressive sensing in speech processing: A survey based on sparsity and sensing matrix | |
Zamani et al. | Separation of nonlinearly mixed sources using end-to-end deep neural networks | |
Zhang et al. | Signal processing using dictionaries, atoms, and deep learning: A common analysis-synthesis framework | |
CN115270847A (zh) | 基于小波包分解和卷积神经网络的设计决策脑电识别方法 | |
Horapong et al. | Progressive latent fingerprint enhancement using two-stage spectrum boosting with matched filter and sparse autoencoder | |
Bartler et al. | Training variational autoencoders with discrete latent variables using importance sampling | |
Chuma et al. | FPGA implementation of a de-noising using Haar level 5 wavelet transform | |
Flaherty et al. | AI to Judge AI-Generated Images: Both Semantics and Perception Matter |
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 |