CN102940490B - 基于非线性动力学的运动想象脑电信号特征提取方法 - Google Patents
基于非线性动力学的运动想象脑电信号特征提取方法 Download PDFInfo
- Publication number
- CN102940490B CN102940490B CN201210434254.5A CN201210434254A CN102940490B CN 102940490 B CN102940490 B CN 102940490B CN 201210434254 A CN201210434254 A CN 201210434254A CN 102940490 B CN102940490 B CN 102940490B
- Authority
- CN
- China
- Prior art keywords
- sample point
- coarse
- vectors
- eeg signals
- subsignal
- 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.)
- Expired - Fee Related
Links
Images
Landscapes
- Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
Abstract
本发明公开了一种基于非线性动力学的运动想象脑电信号特征提取方法。该方法解决了运动想象脑电信号在实际应用环境中可能出现信噪比低从而导致分类正确率低和不稳定的问题。本发明具体步骤包括:(1)采集数据:(2)空间滤波;(3)基线校正;(4)带通滤波;(5)特征提取;(6)分类。本发明具有能够有效保证稳定的脑电信号分类正确率、有效的提供更好实时性的优点。
Description
技术领域
本发明属于信息技术领域,更一步涉及在生命科学领域中应用脑-机接口(Brain-Computer Interface,BCI)系统的基于非线性动力学的运动想象脑电信号特征提取方法。本发明用于对运动想象脑电信号特征的提取,通过分类器对特征进行分类,实现单侧手指运动想象的判别,最终可应用于轮椅、鼠标等运动想象BCI在线系统,从而提高有缺陷人士与外界的交流能力。
背景技术
在准备和执行单侧手指运动想象时,人们大脑皮层的功能性连接发生改变,从而导致其对侧大脑运动感觉区mu和beta节律的脑电信号能量减弱,而其同侧大脑运动感觉区mu和beta节律的脑电信号能量增强。这种想象单侧手指运动时特定脑区特定频率脑电信号的能量变化,被称为事件相关去同步现象。该现象是判别左右手指运动想象脑电信号最根本的特征。因此,通过受试者运动想象时脑电信号的分析,判别受试者运动想象的方向,从而实现对外界装置的控制。
北京工业大学拥有的专利技术“想象单侧肢体运动的脑电特征的提取方法”(申请号200810056839.1,申请日2008.01.25,授权号CN101219048B,授权日2010.06.23)中提出一种想象单侧肢体运动的脑电特征的提取方法。该专利技术将共同空间模式(Common SpatialPattern,CSP)方法和线性判别分析(FDA)相结合,降低了输入矢量的维数,提高了分类器的推广性,在一定程度上提高了分类正确率。共同空间模式方法对于信噪比高的运动想象脑电信号的分类虽然能获得很好的正确率,但是该专利技术仍然存在的不足是,在实际应用中,通常不具备实验室较为理想的数据采集环境,所以实际应用采集的运动想象脑电信号往往与实验室采集的数据相比包含大量噪声,信噪比低。用共同空间模式方法分析处理实际应用采集的运动想象脑电信号的正确率不高,稳定性差,算法运算量大,不利于实际应用的推广。
发明内容
本发明的目的在于克服上述已有脑电信号特征提取方法的不足,提出一种基于非线性动力学的运动想象脑电信号特征的提取方法。该方法充分考虑运动想象时脑电信号在实际环境中信噪比低的特征,保证了一定分类正确率。
实现本发明方法的主要思路是:将采集的多通道脑电信号做空间滤波、基线校正、带通滤波处理,用非线性动力学的多尺度熵方法提取脑电信号特征向量;利用朴素贝叶斯分类器分类,采取交叉验证方法,得到分类正确率。
根据上述主要思路,本发明方法的具体步骤如下:
(1)采集数据:
脑电信号采集系统通过受试者佩戴的电极帽,采集受试者想象单侧手指运动的脑电信号,受试者手指不能有实际上的动作。
(2)空间滤波:
采用共同平均参考的方法,将受试者电极帽上每个电极采集的脑电信号减去所有电极采集的脑电信号的均值,得到共同平均参考后的脑电信号。
(3)基线校正:
将共同平均参考后的脑电信号减去基线,得到基线校正后的脑电信号。
(4)带通滤波:
利用有限脉冲响应滤波器,对基线校正后的脑电信号进行带通滤波,得到不同频带的子信号。
(5)特征提取:
对所有不同频带的子信号在不同的尺度下进行粗粒化处理,对经过粗粒化之后的所有子信号进行样本熵估计,将所有得到的样本熵估计值组成特征向量。
(6)分类:
利用朴素贝叶斯分类器,对得到的特征向量进行分类,采取交叉验证方法对分类结果进行验证,得到分类正确率。
本发明与现有技术相比具有如下优点:
第一,本发明由于考虑了运动想象时脑电信号在实际环境中信噪比低的情况,利用多尺度熵作为特征向量,克服了现有技术中要求运动想象脑电信号信噪比高的、稳定性差的不足,使得本发明具有能够有效保证稳定的脑电信号分类正确率的优点。
第二,本发明由于考虑了运动想象时脑电信号处理的实时性,利用多尺度熵作为特征向量,算法运算量小,克服了现有技术中算法运算量大的问题,使得本发明具有能够有效的提供更好实时性的优点。
附图说明
图1为本发明的流程图;
图2为本发明采集数据步骤中屏幕提示符示意图;
图3为本发明采集数据步骤的示意图。
具体实施方式
下面结合附图1对本发明做进一步的描述。
步骤1,采集数据:
脑电信号采集系统通过受试者佩戴的电极帽,采集想象单侧手指运动的脑电信号。脑电信号由佩戴在受试者头上的电极帽获取,并通过脑电放大器放大以及模/数转换器转换,输入计算机,以信号电压幅值的形式存储并显示。
受试者佩戴电极帽,坐在椅子上平视距其1m左右的显示器。脑电信号采集系统的采样频率为250Hz,测试电极分别为C3和C4,脑电信号的波动值域为±100μV。
在本发明采集数据步骤中需要给受试者运动想象提示,具体的提示种类参照图2。在图2中显示器的提示符有三种,显示器中标有十字叉的代表准备提示,显示器中标有向左的箭头代表想象左手运动提示,向右的箭头代表想象右手运动提示。
在本发明采集数据实施例中图2所示的运动想象提示按设定的时间次序出现,参考图3。在图3中开始0秒至3秒时屏幕出现准备提示,并在第2秒出现短暂的提示音(1kHz,70ms)。随后,屏幕出现想象左手运动提示或想象右手提示,并持续1.25秒。在第4秒,受试者开始想象相应的手指运动,并持续3秒。每次实施例后休息1.5秒至2.5秒。左右手指运动想象实验各120次,顺序上随机排列。步骤2,空间滤波:
采用共同平均参考的方法,对得到的脑电信号计算每个电极信号所有样本点的均值,将所有电极脑电信号的样本点的均值做平均得到总体样本点的均值,将所有脑电信号的每个样本点减去总体样本点的均值,从而得到共同平均参考之后的脑电信号。
步骤3,基线校正:
以受试者想象单侧手指运动之前的200ms脑电信号为基线,将共同平均参考空间滤波后的脑电信号减去基线,得到基线校正后的脑电信号。
步骤4,带通滤波:
利用有限脉冲响应滤波器,对基线校正后的脑电信号进行带通滤波,得到特定频带的子信号,重复执行带通滤波,得到频带为8-13Hz的子信号和18-26Hz子信号。
步骤5,特征提取:
对采集的多个电极中的C3和C4电极,分别通过步骤4的带通滤波,得到两个电极、两个频带的4个子信号。将这4个子信号分别在尺度1、尺度2、尺度3下进行粗粒化处理,得到12个粗粒化序列。本发明实施例中,尺度1的数值取1,尺度2的数值取2,尺度3的数值取3。对得到的12个粗粒化序列分别计算在阈值r=0.15×SD下的样本熵,其中SD为求样本熵的信号的标准差,在本发明的实施例中先通过matlab中的var函数取得信号的方差,将方差开平方得到该标准差。将得到的样本熵作为特征值,从而得到12个特征值,将这12个特征值组成一个特征向量。
对两个不同电极、两个不同频带的子信号在不同尺度下进行粗粒化处理的具体步骤是:
第一步:按照下式计算粗粒化样本点;
其中,yj (τ)为粗粒化样本点;τ为尺度的值;j为粗粒化样本点下标,代表第j个粗粒化样本点,其范围是1≤j≤Q/τ;Q为子信号的样本点总数;i为子信号的第i个样本点,其范围是1≤i≤Q,其初值为(j-1)τ+1;xi为子信号的样本点;
第二步:按照粗粒化样本点下标从小到大的顺序,将得到的粗粒化样本点构建成粗粒化序列。
对粗粒化处理得到的粗粒化序列进行样本熵估计方法如下:
第一步:按照下式将粗粒化序列的样本点组成一个m维矢量组:
um(p)={yp,yp+1,...yp+m-1},1≤p≤N-m
其中:um(p)为一个m维的矢量组;m为矢量组的维数;p为粗粒化序列的第p个样本点,其范围是1≤p≤N-m;N为粗粒化序列的样本点总数;yp为粗粒化序列的样本点;
第二步:按照下式计算两个m维矢量组之间的距离:
d=max{|y(k+a)-y(l+a)|:0≤a≤m-1,k≠l}
其中:d为两个m维矢量组之间的距离;m为矢量组的维数;max{}为求最大值的函数;y(k+a)为矢量组um(k)中第a个元素;y(l+a)为矢量组um(l)中第a个元素;um(k)为第k个矢量组;um(l)为第l个矢量组;
第三步:根据阈值r,对每个p≤N-m的值,统计距离d小于阈值r的数目(称为模板匹配数)及此数目与距离总数N-m-1的比值,该比值记作
第四步:按照下式求模板匹配数与距离总个数的比值对矢量组总数的平均值
其中:Bm(r)为模板匹配总数与距离总个数的比值对矢量组总数的平均值;r为阈值;m为矢量组的维数;N为粗粒化序列的样本点总数;为模板匹配数与其距离总个数的比值;c为第c个模板匹配数与距离总个数的比值,其范围是1≤c≤N-m;
第五步:增大维数m重复第一步、第二步、第三步、第四步,得到m+1维的模板匹配数与距离总个数的比值对矢量组总数的平均值Bm+1(r);
其中:Bm+1(r)为模板匹配总数与距离总个数的比值对矢量组总数的平均值;r为阈值;m为矢量组的维数;N为粗粒化序列的样本点总数;为模板匹配数与其距离总个数的比值;c为第c个模板匹配数与距离总个数的比值,其范围是1≤c≤N-m;
第六步:熵估计值为
SampEn(m,r,N)=-ln[Bm+1(r)/Bm(r)]
其中:SampEn(m,r,N)为熵估计值的表达式;m为矢量组的维数;r为阈值;N为粗粒化序列的样本点总数;ln()为对数函数;Bm+1(r)为m+1维的模板匹配数与距离总数的比值;Bm(r)为m维的模板匹配数与距离总数的比值。
步骤6,特征分类:
利用朴素贝叶斯分类器,对特征向量分类,采取交叉验证方法,对分类结果进行交叉验证,从而得到分类正确率。
选择步骤5得到的任意一个脑电信号的特征向量作为测试数据,其余脑电信号的特征向量作为训练数据;朴素贝叶斯分类器利用训练数据建立分类模型,将测试数据代入分类模型得到分类类别,比较分类类别与实际任务类别,得到正确分类或错误分类的结果;依次将每个脑电信号的特征向量作为一次测试数据,统计所有测试数据的分类结果,得到分类正确率。
Claims (6)
1.基于非线性动力学的运动想象脑电信号特征提取方法,具体步骤如下:
(1)采集数据:
脑电信号采集系统通过受试者佩戴的电极帽,采集受试者想象单侧手指运动时的电极帽中多个电极的脑电信号;
(2)空间滤波:
采用共同平均参考的方法,对采集到的多个电极的脑电信号进行处理,得到共同平均参考后的脑电信号;
(3)基线校正:
将共同平均参考后的脑电信号减去基线,得到基线校正后的脑电信号;
(4)带通滤波:
利用有限脉冲响应滤波器,对基线校正后的脑电信号进行带通滤波,得到特定频带的子信号,重复执行带通滤波,得到两个不同频带的子信号;
(5)特征提取:
对两个不同频带的子信号在不同尺度下进行粗粒化处理,对粗粒化处理得到的粗粒化序列进行样本熵估计,将两个不同频带、不同尺度下得到的样本熵估计值组成特征向量;
所述的粗粒化的具体步骤是:
第一步:按照下式计算粗粒化样本点:
其中,yj (τ)为粗粒化样本点;τ为尺度的值;j为粗粒化样本点的序号,代表第j个粗粒化样本点,其范围是1≤j≤Q/τ;Q为子信号的样本点总数;i为子信号的第i个样本点,其范围是1≤i≤Q,其初值为(j-1)τ+1;xi为子信号的样本点;
第二步:按照粗粒化样本点的序号从小到大的顺序,将得到的粗粒化样本点构建成粗粒化序列;
所述的信号的样本熵估计方法如下;
第一步:按照下式将粗粒化序列的样本点组成一个矢量组:
um(p)={yp,yp+1,...yp+m-1}
其中:um(p)为一个m维的矢量组;m为矢量组的维数;p为粗粒化序列的第p个样本点,其范围是1≤p≤N-m;N为粗粒化序列的样本点总数;yp为粗粒化序列的样本点;
第二步:按照下式计算两个m维矢量组之间的距离:
d=max{|y(k+a)-y(l+a)|}
其中:d为两个m维矢量组之间的距离;m为矢量组的维数;max{}为求最大值的函数;y(k+a)为矢量组um(k)中第a个元素;y(l+a)为矢量组um(l)中第a个元素;um(k)为第k个矢量组;um(l)为第l个矢量组,k和l互不相等;
第三步:设定阈值r,r=0.15×SD,其中SD为子信号的标准方差,统计距离d小于阈值r的矢量组的数目并记为模板匹配数,计算模板匹配数与距离总个数的比值;
第四步:按照下式求模板匹配数与距离总个数的比值对矢量组总数的平均值:
其中:Bm(r)为模板匹配总数与距离总个数的比值对矢量组总数的平均值;r为阈值;m为矢量组的维数;N为粗粒化序列的样本点总数;为模板匹配数与其距离总个数的比值;c为第c个模板匹配数与距离总个数的比值,其范围是1≤c≤N-m;
第五步:增大维数m重复第一步、第二步、第三步、第四步,得到m+1维的模板匹配数与距离总个数的比值对矢量组总数的平均值Bm+1(r);
第六步:根据下式计算熵估计值:
SampEn(m,r,N)=-ln[Bm+1(r)/Bm(r)]
其中:SampEn(m,r,N)为熵估计值的表达式;m为矢量组的维数;r为阈值;N为粗粒化序列的样本点总数;ln()为对数函数;Bm+1(r)为m+1维的模板匹配数与距离总数的比值;Bm(r)为m维的模板匹配数与距离总数的比值;
(6)分类:
利用朴素贝叶斯分类器,对得到的特征向量进行分类,采取交叉验证方法对分类结果进行验证,得到分类正确率。
2.根据权利要求1所述的基于非线性动力学的运动想象脑电信号特征提取方法,其特征在于:步骤(2)中所述的共同平均参考方法如下:
第一步:对步骤(1)得到的每个电极的脑电信号的所有样本点取均值,得到每个电极的均值;
第二步:将第一步得到的每个电极的均值求和再除以电极数目,得到总体样本点的均值;
第三步:将所有电极的脑电信号的每个样本点分别减去总体样本点的均值。
3.根据权利要求1所述的基于非线性动力学的运动想象脑电信号特征提取方法,其特征在于:步骤(3)所述的基线为受试者想象单侧手指运动之前的200ms脑电信号。
4.根据权利要求1所述的基于非线性动力学的运动想象脑电信号特征提取方法,其特征在于:步骤(4)所述的不同频带的子信号包括8-13Hz的子信号以及18-26Hz子信号。
5.根据权利要求1所述的基于非线性动力学的运动想象脑电信号特征提取方法,其特征在于:步骤(5)中所述不同的尺度是指尺度的数量不少于3个,且相互不重复。
6.根据权利要求1所述的运动想象脑电信号特征的提取方法,其特征在于:步骤(6)中所述的交叉验证方法是指,选择步骤(5)得到的任意一个脑电信号的特征向量作为测试数据,其余脑电信号的特征向量作为训练数据;分类器利用训练数据建立分类模型,将测试数据代入分类模型得到分类类别,比较分类类别与实际任务类别,得到正确分类或错误分类的结果;依次将每个脑电信号的特征向量作为一次测试数据,统计所有测试数据的分类结果,得到分类正确率。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210434254.5A CN102940490B (zh) | 2012-10-19 | 2012-10-19 | 基于非线性动力学的运动想象脑电信号特征提取方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210434254.5A CN102940490B (zh) | 2012-10-19 | 2012-10-19 | 基于非线性动力学的运动想象脑电信号特征提取方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102940490A CN102940490A (zh) | 2013-02-27 |
CN102940490B true CN102940490B (zh) | 2014-06-18 |
Family
ID=47723563
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210434254.5A Expired - Fee Related CN102940490B (zh) | 2012-10-19 | 2012-10-19 | 基于非线性动力学的运动想象脑电信号特征提取方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102940490B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105956624A (zh) * | 2016-05-06 | 2016-09-21 | 东南大学 | 基于空时频优化特征稀疏表示的运动想象脑电分类方法 |
Families Citing this family (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103258215A (zh) * | 2013-05-10 | 2013-08-21 | 杭州电子科技大学 | 一种多导联间相关性分析的脑电特征提取方法 |
CN103340624B (zh) * | 2013-07-22 | 2015-01-14 | 上海交通大学 | 少通道下运动想象脑电特征的提取方法 |
CN104083163A (zh) * | 2014-07-16 | 2014-10-08 | 南京大学 | 一种非线性参数脑电地形图的方法 |
CN106108893B (zh) * | 2016-06-20 | 2019-03-01 | 杭州电子科技大学 | 基于眼电、脑电的运动想象训练人机交互系统设计方法 |
CN106156737A (zh) * | 2016-07-05 | 2016-11-23 | 北京工业大学 | 基于二维空间的时间序列样本熵的计算方法及系统 |
CN106473736B (zh) * | 2016-10-11 | 2019-05-21 | 天津大学 | 基于复杂网络的脑电信号分析方法及应用 |
CN107095668B (zh) * | 2017-04-14 | 2019-10-01 | 陇东学院 | 一种基于时域分析的自适应脑电波信号异常的检测方法 |
WO2020042511A1 (zh) * | 2018-08-28 | 2020-03-05 | 天津大学 | 基于空间滤波与模版匹配的运动电位脑机接口编解码方法 |
CN109192007A (zh) * | 2018-09-21 | 2019-01-11 | 杭州电子科技大学 | 一种基于肌电运动感知的ar手语识别方法及教学方法 |
CN109620215A (zh) * | 2018-12-14 | 2019-04-16 | 中航华东光电(上海)有限公司 | 一种脑电特征的提取方法 |
CN109871882A (zh) * | 2019-01-24 | 2019-06-11 | 重庆邮电大学 | 基于高斯伯努利卷积深度置信网络的脑电信号分类方法 |
CN110652293A (zh) * | 2019-10-22 | 2020-01-07 | 燕山大学 | 一种用于运动想象脑电信号的自适应预处理优化方法 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101980106A (zh) * | 2010-10-15 | 2011-02-23 | 华南理工大学 | 一种脑机接口的二维光标控制方法及装置 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20050085744A1 (en) * | 2003-10-20 | 2005-04-21 | Stmicroelectronics S.R.I. | Man-machine interfaces system and method, for instance applications in the area of rehabilitation |
EP2126785A2 (en) * | 2007-01-25 | 2009-12-02 | NeuroVista Corporation | Systems and methods for identifying a contra-ictal condition in a subject |
-
2012
- 2012-10-19 CN CN201210434254.5A patent/CN102940490B/zh not_active Expired - Fee Related
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101980106A (zh) * | 2010-10-15 | 2011-02-23 | 华南理工大学 | 一种脑机接口的二维光标控制方法及装置 |
Non-Patent Citations (3)
Title |
---|
基于运动想象的脑电信号特征提取与分类;李丽君等;《医疗卫生装备 》;20110115;第333-336页 * |
张毅等.基于AR和SVM的运动想象脑电信号识别.《华中科技大学学报(自然科学版)》.2011,(第39期),第103-106页. * |
李丽君等.基于运动想象的脑电信号特征提取与分类.《医疗卫生装备 》.2011,第333-336页. |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105956624A (zh) * | 2016-05-06 | 2016-09-21 | 东南大学 | 基于空时频优化特征稀疏表示的运动想象脑电分类方法 |
CN105956624B (zh) * | 2016-05-06 | 2019-05-21 | 东南大学 | 基于空时频优化特征稀疏表示的运动想象脑电分类方法 |
Also Published As
Publication number | Publication date |
---|---|
CN102940490A (zh) | 2013-02-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102940490B (zh) | 基于非线性动力学的运动想象脑电信号特征提取方法 | |
CN101596101B (zh) | 依据脑电信号判定疲劳状态的方法 | |
CN108143409B (zh) | 睡眠阶段分期方法及装置 | |
CN100482155C (zh) | 基于脑机交互的注意力状态即时检测系统及检测方法 | |
CN110353702A (zh) | 一种基于浅层卷积神经网络的情感识别方法及系统 | |
CN103610447B (zh) | 一种基于前额脑电信号的脑力负荷在线检测方法 | |
CN104771163A (zh) | 基于csp和r-csp算法的脑电信号特征提取方法 | |
CN109359619A (zh) | 一种基于卷积盲源分离的高密度表面肌电信号分解方法 | |
CN103699230A (zh) | 一种基于图标脑电控制的数字界面交互方法 | |
CN107595247A (zh) | 一种基于脑电信号的麻醉深度的监测方法及系统 | |
CN101515200A (zh) | 基于瞬态视觉诱发脑电的目标选择方法 | |
CN103300852A (zh) | 基于运动想象脑电控制遥控小车的训练方法 | |
CN105640500A (zh) | 基于独立分量分析的扫视信号特征提取方法和识别方法 | |
CN103019382A (zh) | 一种利用诱发电位反映脑部主观动机信号的脑机接口方法 | |
CN113017645B (zh) | 一种基于空洞卷积神经网络的p300信号检测方法 | |
CN104571504A (zh) | 一种基于想象动作的在线脑-机接口方法 | |
CN113017627A (zh) | 一种基于双通道相位同步特征融合的抑郁症和双相障碍脑网络分析方法 | |
CN106648087A (zh) | 一种基于意识任务的特征eeg处理方法 | |
Barban et al. | Yet another artefact rejection study: an exploration of cleaning methods for biological and neuromodulatory noise | |
CN109009098A (zh) | 一种运动想象状态下的脑电信号特征识别方法 | |
CN106333676B (zh) | 清醒状态下的脑电信号数据类型的标注装置 | |
CN106073767B (zh) | Eeg信号的相位同步度量、耦合特征提取及信号识别方法 | |
CN103892797B (zh) | 一种用于睡眠结构分析的信号处理方法和装置 | |
Dzitac et al. | Identification of ERD using fuzzy inference systems for brain-computer interface | |
Jia | Neural network in the application of EEG signal classification method |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20140618 Termination date: 20191019 |
|
CF01 | Termination of patent right due to non-payment of annual fee |