CN114145759A - 肌电信号的补偿方法、装置、肌电检测设备及存储介质 - Google Patents
肌电信号的补偿方法、装置、肌电检测设备及存储介质 Download PDFInfo
- Publication number
- CN114145759A CN114145759A CN202111337718.6A CN202111337718A CN114145759A CN 114145759 A CN114145759 A CN 114145759A CN 202111337718 A CN202111337718 A CN 202111337718A CN 114145759 A CN114145759 A CN 114145759A
- Authority
- CN
- China
- Prior art keywords
- electromyographic
- electrode
- data
- electrodes
- sampling time
- 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
- 238000000034 method Methods 0.000 title claims abstract description 35
- 238000001514 detection method Methods 0.000 title claims abstract description 34
- 230000003183 myoelectrical effect Effects 0.000 title description 14
- 238000005070 sampling Methods 0.000 claims description 152
- 238000002567 electromyography Methods 0.000 claims description 20
- 238000013507 mapping Methods 0.000 claims description 20
- 238000004590 computer program Methods 0.000 claims description 13
- 238000010586 diagram Methods 0.000 description 14
- 210000003205 muscle Anatomy 0.000 description 14
- 238000012545 processing Methods 0.000 description 10
- 238000004891 communication Methods 0.000 description 7
- 230000006870 function Effects 0.000 description 7
- 230000000694 effects Effects 0.000 description 5
- 230000009471 action Effects 0.000 description 4
- 230000005540 biological transmission Effects 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 2
- 239000000835 fiber Substances 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 230000001902 propagating effect Effects 0.000 description 2
- 230000005236 sound signal Effects 0.000 description 2
- 239000000126 substance Substances 0.000 description 2
- RYGMFSIKBFXOCR-UHFFFAOYSA-N Copper Chemical compound [Cu] RYGMFSIKBFXOCR-UHFFFAOYSA-N 0.000 description 1
- 230000002159 abnormal effect Effects 0.000 description 1
- 229910052802 copper Inorganic materials 0.000 description 1
- 239000010949 copper Substances 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 230000014509 gene expression Effects 0.000 description 1
- 238000012417 linear regression Methods 0.000 description 1
- 239000004973 liquid crystal related substance Substances 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000013307 optical fiber Substances 0.000 description 1
- 230000002093 peripheral effect Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 230000003238 somatosensory effect Effects 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 238000012360 testing method 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/24—Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
- A61B5/316—Modalities, i.e. specific diagnostic methods
- A61B5/389—Electromyography [EMG]
-
- 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/7225—Details of analog processing, e.g. isolation amplifier, gain or sensitivity adjustment, filtering, baseline or drift compensation
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Surgery (AREA)
- General Health & Medical Sciences (AREA)
- Biophysics (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Medical Informatics (AREA)
- Molecular Biology (AREA)
- Physics & Mathematics (AREA)
- Animal Behavior & Ethology (AREA)
- Pathology (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Signal Processing (AREA)
- Power Engineering (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Physiology (AREA)
- Psychiatry (AREA)
- Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
Abstract
本公开提供了一种肌电信号的补偿方法、装置、肌电检测设备及存储介质，该方法包括：获取肌电检测设备上设置的至少两个电极所采集的肌电信号；在第一电极所采集的肌电信号出现数据缺失的情况下，确定其他电极中与所述第一电极相关的第二电极；其中，所述第一电极为所述至少两个电极中的任一个电极，所述其他电极为所述至少两个电极中除所述第一电极以外的电极；根据所述第二电极所采集的肌电信号，确定所述第一电极所采集的肌电信号中缺失的肌电数据，作为目标肌电数据；根据所述目标肌电数据对所述第一电极所采集的肌电信号进行补偿。
Description
技术领域
本公开涉及信号处理技术领域，更具体地，涉及一种肌电信号的补偿方法、一种肌电信号的补偿装置、一种肌电检测设备及一种可读存储介质。
背景技术
目前的肌电检测设备，通常设置有多个电极，并通过多个电极来同时采集测试者不同部位的肌电信号。
但在实际测试中，由于测试者的肌肉比较发达，在做某些动作时可能会使部分电极短暂地与人体脱离接触，这样就会导致瞬时其采集的肌电信号异常。或者，由于系统出现阻塞或者延迟，导致肌电信号中的部分数据的丢失。
如果肌电信号中出现数据丢失的问题，可能会导致无法根据肌电信号准确地识别测试者的动作，因此，提出一种对肌电信号进行补偿的方法是十分有必要的。
发明内容
本公开的一个目的是提供一种对肌电信号进行补偿的新技术方案。
根据本公开的第一方面，提供了一种肌电信号的补偿方法，包括：
获取肌电检测设备上设置的至少两个电极所采集的肌电信号；
在第一电极所采集的肌电信号出现数据缺失的情况下，确定其他电极中与所述第一电极相关的第二电极；其中，所述第一电极为所述至少两个电极中的任一个电极，所述其他电极为所述至少两个电极中除所述第一电极以外的电极；
根据所述第二电极所采集的肌电信号，确定所述第一电极所采集的肌电信号中缺失的肌电数据，作为目标肌电数据；
根据所述目标肌电数据对所述第一电极所采集的肌电信号进行补偿。
可选的，所述确定其他电极中与所述第一电极相关的第二电极，包括：
根据所述至少两个电极所采集的肌电信号，确定所述第一电极与所述其他电极之间的目标相关系数；
确定与所述第一电极之间的目标相关系数最大的其他电极，作为所述第二电极。
可选的，所述根据所述至少两个电极所采集的肌电信号，确定所述第一电极与所述其他电极中之间的目标相关系数，包括：
获取所述第一电极所采集的肌电信号出现数据缺失前的n个第一采样时刻的第一肌电数据；其中，n为正整数；
获取所述第一电极与所述其他电极的采样时刻差；其中，所述采样时刻差指示所述第一电极所采集的肌电信号与所述其他电极所采集的肌电信号中，相匹配的肌电数据的采样时刻之间的差异；
获取每一所述采样时刻差下，所述其他电极所采集的肌电信号中与所述第一肌电数据匹配的其他肌电数据；
根据所述第一肌电数据和每一所述采样时刻差下的所述其他肌电数据，确定所述第一电极与所述其他电极之间的目标相关系数。
可选的，所述获取每一所述采样时刻差下，所述其他电极所采集的肌电信号中与所述第一肌电数据匹配的其他肌电数据，包括：
遍历所述采样时刻差；
根据当前遍历的采样时刻差，得到与所述第一采样时刻对应的第二采样时刻；
获取所述其他电极所采集的肌电信号中所述第二采样时刻的肌电数据，作为所述当前遍历的采样时刻差下，与所述第一肌电数据匹配的所述其他肌电数据。
可选的，所述根据所述第一肌电数据和每一所述采样时刻差下的所述其他肌电数据，确定所述第一电极与所述其他电极之间的目标相关系数，包括：
根据所述第一肌电数据和每一所述采样时刻差下的所述其他肌电数据，确定所述第一电极与所述其他电极在每一所述采样时刻差下的相关系数；
选取最大的相关系数，作为所述第一电极与所述其他电极之间的目标相关系数。
可选的，所述根据所述第二电极所采集的肌电信号，确定所述第一电极所采集的肌电信号中缺失的肌电数据，作为目标肌电数据，包括：
获取所述第一电极采集的肌电信号和所述第二电极所采集的肌电信号中，相匹配的肌电数据之间的映射关系；
获取缺失的肌电数据所对应的所述第一电极的采样时刻，作为第三采样时刻；
根据目标采样时刻差，确定与所述第三采样时刻对应的第四采样时刻；其中，所述目标采样时刻差为所述目标相关系数所对应的采样时刻差；
获取所述第二电极所采集的肌电信号中所述第四采样时刻的肌电数据，作为参考肌电数据；
根据所述参考肌电数据和所述映射关系，得到所述目标肌电数据。
可选的，所述获取所述第一电极采集的肌电信号和所述第二电极所采集的肌电信号中，相匹配的肌电数据之间的映射关系，包括：
获取所述目标采样时刻差下，所述第二电极所采集的肌电信号中与所述第一肌电数据匹配的第二肌电数据；
根据所述第一肌电数据和所述第二肌电数据，得到所述映射关系。
根据本公开的第二方面，提供了一种肌电信号的补偿装置，包括：
信号获取模块，用于获取肌电检测设备上设置的至少两个电极所采集的肌电信号；
相关电极确定模块，用于在第一电极所采集的肌电信号出现数据缺失的情况下，确定其他电极中与所述第一电极相关的第二电极；其中，所述第一电极为所述至少两个电极中的任一个电极，所述其他电极为所述至少两个电极中除所述第一电极以外的电极；
缺失数据确定模块，用于根据所述第二电极所采集的肌电信号，确定所述第一电极所采集的肌电信号中缺失的肌电数据，作为目标肌电数据；
信号补偿模块，用于根据所述目标肌电数据对所述第一电极所采集的肌电信号进行补偿。
根据本公开的第三方面，提供了一种肌电检测设备，包括：
如本公开第二方面所述的装置；或者，
处理器和存储器，所述存储器用于存储指令，所述指令用于控制所述处理器执行根据本公开第一方面所述的方法。
根据本公开的第四方面，提供了一种可读存储介质，其上存储有计算机程序，所述计算机程序在被处理器执行时实现如本公开第一方面所述的方法。
通过本实施例的方法，在第一电极所采集的肌电信号出现数据缺失的情况下，确定与第一电极相关的第二电极；根据第二电极所采集的肌电信号来确定第一电极所采集的肌电信号中所缺失的肌电数据，作为目标肌电数据，根据目标肌电数据对第一电极所采集的肌电信号进行补偿。本实施例基于人体在执行动作时，肌肉之间具有相关性的特点，在多个电极中利用相关系数找出与出现数据缺失的第一电极最为相似的第二电极，利用第一电极所采集的肌电信号中丢失数据前的一段数据、以及第二电极所采集的肌电信号，来求出用于替代丢失的数据的目标肌电数据，对第一电极所采集的肌电信号进行补偿，可以使得补偿后的肌电信号更加接近真实数据。而且，在肌电信号中缺失的肌电数据较多的情况下，根据本实施例的方法所补偿得到的肌电信号的准确度也较高。
通过以下参照附图对本公开的示例性实施例的详细描述，本公开的其它特征及其优点将会变得清楚。
附图说明
被结合在说明书中并构成说明书的一部分的附图示出了本公开的实施例，并且连同其说明一起用于解释本公开的原理。
图1是显示可用于实现本公开的实施例的肌电检测设备的硬件配置的一个例子的框图。
图2示出了本公开的实施例的肌电信号的补偿方法的流程图。
图3示出了本公开的实施例的肌电信号的补偿装置的框图。
图4示出了本公开的实施例的肌电检测设备的框图。
具体实施方式
现在将参照附图来详细描述本公开的各种示例性实施例。应注意到：除非另外具体说明，否则在这些实施例中阐述的部件和步骤的相对布置、数字表达式和数值不限制本公开的范围。
以下对至少一个示例性实施例的描述实际上仅仅是说明性的，决不作为对本公开及其应用或使用的任何限制。
对于相关领域普通技术人员已知的技术、方法和设备可能不作详细讨论，但在适当情况下，所述技术、方法和设备应当被视为说明书的一部分。
在这里示出和讨论的所有例子中，任何具体值应被解释为仅仅是示例性的，而不是作为限制。因此，示例性实施例的其它例子可以具有不同的值。
应注意到：相似的标号和字母在下面的附图中表示类似项，因此，一旦某一项在一个附图中被定义，则在随后的附图中不需要对其进行进一步讨论。
<硬件配置>
图1是可用于实现本公开实施例的肌电检测设备的硬件配置的一个例子的框图。
该肌电检测设备1000可以包括但不限于处理器1100、存储器1200、接口装置1300、通信装置1400、显示装置1500、输入装置1600、扬声器1700、麦克风1800等等。其中，处理器1100可以是中央处理器CPU、图形处理器GPU、微处理器MCU等，用于执行计算机程序，该计算机程序可以采用比如x86、Arm、RISC、MIPS、SSE等架构的指令集编写。存储器1200例如包括ROM(只读存储器)、RAM(随机存取存储器)、诸如硬盘的非易失性存储器等。接口装置1300例如包括USB接口、串行接口、并行接口等。通信装置1400例如能够利用光纤或电缆进行有线通信，或者进行无线通信，具体地可以包括WiFi通信、蓝牙通信、2G/3G/4G/5G通信等。显示装置1500例如是液晶显示屏、触摸显示屏等。输入装置1600例如可以包括触摸屏、键盘、体感输入等。扬声器1700用于输出音频信号。麦克风1800用于采集音频信号。
应用于本公开实施例中，肌电检测设备1000的存储器1200用于存储计算机程序，该计算机程序用于控制所述处理器1100进行操作以实现根据本公开实施例的方法。技术人员可以根据本公开所公开方案设计该计算机程序。该计算机程序如何控制处理器进行操作，这是本领域公知，故在此不再详细描述。该肌电检测设备1000可以安装有智能操作系统(例如Windows、Linux、安卓、IOS等系统)和应用软件。
本领域技术人员应当理解，尽管在图1中示出了肌电检测设备1000的多个装置，但是，本公开实施例的肌电检测设备1000可以仅涉及其中的部分装置，例如，只涉及处理器1100和存储器1200等。
下面，参照附图描述根据本发明的各个实施例和例子。
<方法实施例>
在本实施例中，提供一种肌电信号的补偿方法。该方法由肌电检测设备实施。在一个例子中，该肌电检测设备可以是图1所示的肌电检测设备1000。
根据图2所示，本实施例的肌电信号的补偿方法可以包括如下步骤S2100～S2400。
步骤S2100，获取肌电检测设备上设置的至少两个电极所采集的肌电信号。
在本实施例中，肌电检测设备上可以设置有至少两个电极，该电极可以是用于采集使用该肌电检测设备的用户的肌电信号。具体的，每个电极可以是按照对应的频率，来采集肌电数据，得到肌电信号。
假设肌电检测设备上设置有四个电极，电极1所采集的肌电信号中，可以包括如下肌电数据：31856.0、31986.0、32103.0、32157.0、32093.0、……、32104.0、32149.0、32145.0、32079.0；电极2所采集的肌电信号中，可以包括如下肌电数据：31298.0、31406.0、31553.0、31646.0、31663.0、……、31551.0、31580.0、31582.0、31538.0；电极3所采集的肌电信号中，可以包括如下肌电数据：32501.0、32534.0、32527.0、32541.0、32513.0、……、32492.0、32472.0、32455.0、32429.0；电极4所采集的肌电信号中，可以包括如下肌电数据：33382.0、33326.0、33344.0、33356.0、33395.0、……、33380.0、33397.0、33360.0、33350.0。
步骤S2200，在第一电极所采集的肌电信号出现数据缺失的情况下，确定其他电极中与第一电极相关的第二电极。
其中，第一电极为至少两个电极中的任一个电极，其他电极为至少两个电极中除第一电极以外的电极。
在第一电极所采集的肌电信号出现数据缺失的情况下，第一电极在其中一个采样时刻所采集的肌电数据为零。
在本公开的一个实施例中，确定其他电极中与第一电极相关的第二电极，可以包括如下所示的步骤S2210～S2220：
步骤S2210，根据至少两个电极所采集的肌电信号，确定第一电极与其他电极之间的目标相关系数。
本实施例中的相关系数，可以是反映两个电极所采集的肌电信号的相关程度。
具体的，用户在做手势1的情况下，可能需要肌肉1和肌肉2同时工作，而用户在做手势2的情况下，可能需要肌肉2和肌肉3同时工作。在肌电检测设备上，可以是设置有分别用于采集反映肌肉1、肌肉2和肌肉3的活动情况的肌电信号的电极。如果电极1用于采集反映肌肉1的活动情况的肌电信号，电极2用于采集反映肌肉2的活动情况的肌电信号，电极3用于采集反映肌肉3的活动情况的肌电信号，那么，用户在做手势1的情况下，电极1和电极2之间相关；用户在做手势2的情况下，电极2和电极3之间相关。
在本公开的一个实施例中，根据至少两个电极所采集的肌电信号，确定第一电极与其他电极之间的目标相关系数，可以包括如下所示的步骤S2211～S2214：
步骤S2211，获取第一电极所采集的肌电信号出现数据缺失前的n个第一采样时刻的第一肌电数据。
在本实施例中，第一电极所采集的肌电信号中出现数据缺失的采样时刻可以是第三采样时刻，即第一电极所采集的肌电信号中，缺失第三采样时刻的肌电数据。第三采样时刻的前n个采样时刻，即为第一采样时刻。其中，n为预先根据应用场景或具体需求所设定的正整数。例如，n可以是20。
第一肌电数据，可以是第一电极在第一采样时刻所采集的肌电数据。具体的，第一肌电数据可以是第一电极所采集的肌电信号中，缺失的肌电数据之前采集的n个肌电数据。
步骤S2212，获取第一电极与其他电极的采样时刻差。
其中，采样时刻差指示第一电极所采集的肌电信号与其他电极所采集的肌电信号中，相匹配的肌电数据的采样时刻之间的差异。
在本实施例中，由于不同电极的采样可能会有时间差，而且，用户在做不同手势时，肌肉的活动情况也可能会有时间差，因此，不同电极之间可以是具有采样时刻差。
进一步地，第一电极与不同电极之间的采样时刻差，或者，用户在做不同手势时第一电极与任一个电极之间的采样时刻差，可能是不同的。因此，本实施例中，可以是预先根据应用场景或具体需求设置多个采样时刻差，以从中选取最能够反映第一电极与其他电极之间的采样延时情况的采样时刻差。
在一个例子中，本实施例中所获取的采样时刻差，可以是-m～m个时间单位，1个时间单位为相邻两个采样时刻差之间的时间间隔。其中，m为预先根据应用场景或具体需求所设定的正整数。
步骤S2213，获取每一采样时刻差下，其他电极所采集的肌电信号中与第一肌电数据匹配的其他肌电数据。
在本公开的一个实施例中，获取每一采样时刻差下，其他电极所采集的肌电信号中与第一肌电数据匹配的其他肌电数据，可以包括如下所示的步骤S22131～S22133：
步骤S22131，遍历采样时刻差。
步骤S22132，根据当前遍历的采样时刻差，得到与第一采样时刻对应的第二采样时刻。
在当前遍历的采样时刻差为δ(其中，δ为大于等于-m且小于等于m的整数)个时间单位，第一采样时刻为第i(其中，i为大于等于1且小于等于n的整数)个采样时刻的情况下，与该第一采样时刻对应的第二采样时刻可以是第i+δ个采样时刻。
步骤S22133，获取其他电极所采集的肌电信号中第二采样时刻的肌电数据，作为当前遍历的采样时刻差下，与第一肌电数据匹配的其他肌电数据。
在当前遍历的采样时刻差为δ个时间单位，第一采样时刻为第i个采样时刻，第一采样时刻的第一肌电数据可以为xi的情况下，当前遍历的采样时刻差下，与第一肌电数据匹配的其他肌电数据可以表示为yi+δ。
步骤S2214，根据第一肌电数据和每一采样时刻差下的其他肌电数据，确定第一电极与其他电极之间的目标相关系数。
在本公开的一个实施例中，根据第一肌电数据和每一采样时刻差下的其他肌电数据，确定第一电极与其他电极之间的目标相关系数，可以包括如下所示的步骤S22141～S22142：
步骤S22141，根据第一肌电数据和每一采样时刻差下的其他肌电数据，确定第一电极与其他电极在每一采样时刻差下的相关系数。
第一电极与其他电极在采样时刻差为采样时刻差为δ个时间单位的情况下的相关系数rδ，可以是通过如下公式确定：
其中，
为n个第一采样时刻的第一肌电数据的平均值，
为n个第二采样时刻的其他肌电数据的平均值。
步骤S22142，选取最大的相关系数，作为第一电极与其他电极之间的目标相关系数。
在本实施例中，可以是选取第一电极与任一个其他电极在所有采样时刻差下最大的相关系数，作为第一电极与个其他电极之间的目标相关系数。
步骤S2220，确定与第一电极之间的目标相关系数最大的其他电极，作为第二电极。
例如，在第一电极为电极1，电极1与电极2之间的目标相关系数可以是0.747781，电极1与电极3之间的目标相关系数可以是0.126012，电极1与电极4之间的目标相关系数可以是0.174313的情况下，电极1与电极2之间的目标相关系数最大，因此，可以是将电极2作为与第一电极相关的第二电极。
步骤S2300，根据第二电极所采集的肌电信号，确定第一电极所采集的肌电信号中缺失的肌电数据，作为目标肌电数据。
在本公开的一个实施例中，根据第二电极所采集的肌电信号，确定第一电极所采集的肌电信号中缺失的肌电数据，作为目标肌电数据，可以包括如下所示的步骤S2310～S2350：
步骤S2310，获取第一电极采集的肌电信号和第二电极所采集的肌电信号中，相匹配的肌电数据之间的映射关系。
本实施例中，第一电极采集的肌电信号和第二电极所采集的肌电信号中，相匹配的肌电数据之间的映射关系，可以是用一元线性回归方程表示。
在一个例子中，该映射关系可以是通过如下公式表示：
z＝a*x+b
其中，x和z分别为第一电极采集的肌电信号和第二电极所采集的肌电信号中相匹配的肌电数据，a和b为系数。该系数a和b可以是通过第一肌电数据、及第二电极所采集的肌电信号中与第一肌电数据匹配的第二肌电数据所确定。
在本公开的一个实施例中，获取第一电极采集的肌电信号和第二电极所采集的肌电信号中，相匹配的肌电数据之间的映射关系，可以包括如下所示的步骤S2311～S2312：
步骤S2311，获取目标采样时刻差下，第二电极所采集的肌电信号中与第一肌电数据匹配的第二肌电数据。
在本实施例中，目标采样时刻差可以是目标相关系数所对应的采样时刻差。也就是说，第一电极和第二电极在目标采样时刻差下的相关系数最大。
在目标采样时刻差为ε个时间单位，第一肌电数据为xi的情况下，与该第一肌电数据匹配的第二肌电数据可以表示为zi+∈。
步骤S2312，根据第一肌电数据和第二肌电数据，得到该映射关系。
在本实施例中，可以通过如下公式来确定系数a：
可以通过如下公式来确定系数b：
其中，
为n个第一采样时刻的第一肌电数据的平均值，
为n个第二采样时刻的第二肌电数据的平均值。
步骤S2320，获取缺失的肌电数据所对应的第一电极的采样时刻，作为第三采样时刻。
步骤S2330，根据目标采样时刻差，确定与第三采样时刻对应的第四采样时刻。
其中，目标采样时刻差为目标相关系数所对应的采样时刻差。
在目标采样时刻差为ε个时间单位，第三采样时刻为第n+1个采样时刻的情况下，与第三采样时刻对应的第四采样时刻可以是第n+ε+1个采样时刻。
步骤S2340，获取第二电极所采集的肌电信号中第四采样时刻的肌电数据，作为参考肌电数据。
在本实施例中，第二电极所采集的肌电信号中第四采样时刻的参考肌电数据，可以是zn+∈+1。
步骤S2350，根据参考肌电数据和映射关系，得到目标肌电数据。
具体的，可以是将参考肌电数据zn+∈+1代入映射关系z＝a*x+b中的z，得到的x即为第一电极所采集的肌电信号中第三采样时刻的所缺失的肌电数据，作为目标肌电数据。
步骤S2400，根据目标肌电数据对第一电极所采集的肌电信号进行补偿。
在本实施例中，可以是将通过步骤S2300所得到的目标肌电数据，补偿至第一电极所采集的肌电信号中对应第三采样时刻的位置上，得到补偿后的肌电信号。
通过本公开的实施例，在第一电极所采集的肌电信号出现数据缺失的情况下，确定与第一电极相关的第二电极；根据第二电极所采集的肌电信号来确定第一电极所采集的肌电信号中所缺失的肌电数据，作为目标肌电数据，根据目标肌电数据对第一电极所采集的肌电信号进行补偿。本实施例基于人体在执行动作时，肌肉之间具有相关性的特点，在多个电极中利用相关系数找出与出现数据缺失的第一电极最为相似的第二电极，利用第一电极所采集的肌电信号中丢失数据前的一段数据、以及第二电极所采集的肌电信号，来求出用于替代丢失的数据的目标肌电数据，对第一电极所采集的肌电信号进行补偿，可以使得补偿后的肌电信号更加接近真实数据。而且，在肌电信号中缺失的肌电数据较多的情况下，根据本实施例的方法所补偿得到的肌电信号的准确度也较高。
<装置实施例>
在本实施例中，提供一种肌电信号的补偿装置3000，如图3所示，包括信号获取模块3100、相关电极确定模块3200、缺失数据确定模块3300和信号补偿模块3400。该信号获取模块3100用于获取肌电检测设备上设置的至少两个电极所采集的肌电信号；该相关电极确定模块3200用于在第一电极所采集的肌电信号出现数据缺失的情况下，确定其他电极中与第一电极相关的第二电极；其中，第一电极为至少两个电极中的任一个电极，其他电极为至少两个电极中除第一电极以外的电极；该缺失数据确定模块3300用于根据第二电极所采集的肌电信号，确定第一电极所采集的肌电信号中缺失的肌电数据，作为目标肌电数据；该信号补偿模块3400用于根据目标肌电数据对第一电极所采集的肌电信号进行补偿。
在本公开的一个实施例中，该相关电极确定模块3200还可以用于：
根据至少两个电极所采集的肌电信号，确定第一电极与其他电极之间的目标相关系数；
确定与第一电极之间的目标相关系数最大的其他电极，作为第二电极。
在本公开的一个实施例中，根据至少两个电极所采集的肌电信号，确定第一电极与其他电极中之间的目标相关系数，包括：
获取第一电极所采集的肌电信号出现数据缺失前的n个第一采样时刻的第一肌电数据；其中，n为正整数；
获取第一电极与其他电极的采样时刻差；其中，采样时刻差指示第一电极所采集的肌电信号与其他电极所采集的肌电信号中，相匹配的肌电数据的采样时刻之间的差异；
获取每一采样时刻差下，其他电极所采集的肌电信号中与第一肌电数据匹配的其他肌电数据；
根据第一肌电数据和每一采样时刻差下的其他肌电数据，确定第一电极与其他电极之间的目标相关系数。
在本公开的一个实施例中，获取每一采样时刻差下，其他电极所采集的肌电信号中与第一肌电数据匹配的其他肌电数据，包括：
遍历采样时刻差；
根据当前遍历的采样时刻差，得到与第一采样时刻对应的第二采样时刻；
获取其他电极所采集的肌电信号中第二采样时刻的肌电数据，作为当前遍历的采样时刻差下，与第一肌电数据匹配的其他肌电数据。
在本公开的一个实施例中，根据第一肌电数据和每一采样时刻差下的其他肌电数据，确定第一电极与其他电极之间的目标相关系数，包括：
根据第一肌电数据和每一采样时刻差下的其他肌电数据，确定第一电极与其他电极在每一采样时刻差下的相关系数；
选取最大的相关系数，作为第一电极与其他电极之间的目标相关系数。
在本公开的一个实施例中，该缺失数据确定模块3300还可以用于：
获取第一电极采集的肌电信号和第二电极所采集的肌电信号中，相匹配的肌电数据之间的映射关系；
获取缺失的肌电数据所对应的第一电极的采样时刻，作为第三采样时刻；
根据目标采样时刻差，确定与第三采样时刻对应的第四采样时刻；其中，目标采样时刻差为目标相关系数所对应的采样时刻差；
获取第二电极所采集的肌电信号中第四采样时刻的肌电数据，作为参考肌电数据；
根据参考肌电数据和映射关系，得到目标肌电数据。
在本公开的一个实施例中，获取第一电极采集的肌电信号和第二电极所采集的肌电信号中，相匹配的肌电数据之间的映射关系，包括：
获取目标采样时刻差下，第二电极所采集的肌电信号中与第一肌电数据匹配的第二肌电数据；
根据第一肌电数据和第二肌电数据，得到映射关系。
本领域技术人员应当明白，可以通过各种方式来实现肌电信号的补偿装置3000。例如，可以通过指令配置处理器来实现肌电信号的补偿装置3000。例如，可以将指令存储在ROM中，并且当启动设备时，将指令从ROM读取到可编程器件中来实现肌电信号的补偿装置3000。例如，可以将肌电信号的补偿装置3000固化到专用器件(例如ASIC)中。可以将肌电信号的补偿装置3000分成相互独立的单元，或者可以将它们合并在一起实现。肌电信号的补偿装置3000可以通过上述各种实现方式中的一种来实现，或者可以通过上述各种实现方式中的两种或更多种方式的组合来实现。
在本实施例中，肌电信号的补偿装置3000可以具有多种实现形式，例如，肌电信号的补偿装置3000可以是任何的提供肌电信号的补偿服务的软件产品或者应用程序中运行的功能模块，或者是这些软件产品或者应用程序的外设嵌入件、插件、补丁件等，还可以是这些软件产品或者应用程序本身。
<肌电检测设备>
在本实施例中，还提供一种肌电检测设备4000。
在一个例子中，该肌电检测设备4000可以包括前述的肌电信号的补偿装置3000。
在另一个例子中，如图4所示，该肌电检测设备4000可以包括：
存储器4100和处理器4200，该存储器4100用于存储可执行的指令；该指令用于控制处理器4200执行前述的肌电信号的补偿方法。
本实施例中的肌电检测设备4000，可以是能够通过设置的多个电极来采集肌电信号的设备。在一个例子中，该肌电检测设备4000可以是可穿戴设备。
<可读存储介质>
在本实施例中，还提供一种可读存储介质，其上存储有计算机程序，计算机程序在被处理器执行时实现如本公开任意实施例的肌电信号的补偿方法。
本公开可以是系统、方法和/或计算机程序产品。计算机程序产品可以包括计算机可读存储介质，其上载有用于使处理器实现本公开的各个方面的计算机可读程序指令。
计算机可读存储介质可以是可以保持和存储由指令执行设备使用的指令的有形设备。计算机可读存储介质例如可以是――但不限于――电存储设备、磁存储设备、光存储设备、电磁存储设备、半导体存储设备或者上述的任意合适的组合。计算机可读存储介质的更具体的例子(非穷举的列表)包括：便携式计算机盘、硬盘、随机存取存储器(RAM)、只读存储器(ROM)、可擦式可编程只读存储器(EPROM或闪存)、静态随机存取存储器(SRAM)、便携式压缩盘只读存储器(CD-ROM)、数字多功能盘(DVD)、记忆棒、软盘、机械编码设备、例如其上存储有指令的打孔卡或凹槽内凸起结构、以及上述的任意合适的组合。这里所使用的计算机可读存储介质不被解释为瞬时信号本身，诸如无线电波或者其他自由传播的电磁波、通过波导或其他传输媒介传播的电磁波(例如，通过光纤电缆的光脉冲)、或者通过电线传输的电信号。
这里所描述的计算机可读程序指令可以从计算机可读存储介质下载到各个计算/处理设备，或者通过网络、例如因特网、局域网、广域网和/或无线网下载到外部计算机或外部存储设备。网络可以包括铜传输电缆、光纤传输、无线传输、路由器、防火墙、交换机、网关计算机和/或边缘服务器。每个计算/处理设备中的网络适配卡或者网络接口从网络接收计算机可读程序指令，并转发该计算机可读程序指令，以供存储在各个计算/处理设备中的计算机可读存储介质中。
用于执行本公开操作的计算机程序指令可以是汇编指令、指令集架构(ISA)指令、机器指令、机器相关指令、微代码、固件指令、状态设置数据、或者以一种或多种编程语言的任意组合编写的源代码或目标代码，所述编程语言包括面向对象的编程语言—诸如Smalltalk、C++等，以及常规的过程式编程语言—诸如“C”语言或类似的编程语言。计算机可读程序指令可以完全地在用户计算机上执行、部分地在用户计算机上执行、作为一个独立的软件包执行、部分在用户计算机上部分在远程计算机上执行、或者完全在远程计算机或服务器上执行。在涉及远程计算机的情形中，远程计算机可以通过任意种类的网络—包括局域网(LAN)或广域网(WAN)—连接到用户计算机，或者，可以连接到外部计算机(例如利用因特网服务提供商来通过因特网连接)。在一些实施例中，通过利用计算机可读程序指令的状态信息来个性化定制电子电路，例如可编程逻辑电路、现场可编程门阵列(FPGA)或可编程逻辑阵列(PLA)，该电子电路可以执行计算机可读程序指令，从而实现本公开的各个方面。
这里参照根据本公开实施例的方法、装置(系统)和计算机程序产品的流程图和/或框图描述了本公开的各个方面。应当理解，流程图和/或框图的每个方框以及流程图和/或框图中各方框的组合，都可以由计算机可读程序指令实现。
这些计算机可读程序指令可以提供给通用计算机、专用计算机或其它可编程数据处理装置的处理器，从而生产出一种机器，使得这些指令在通过计算机或其它可编程数据处理装置的处理器执行时，产生了实现流程图和/或框图中的一个或多个方框中规定的功能/动作的装置。也可以把这些计算机可读程序指令存储在计算机可读存储介质中，这些指令使得计算机、可编程数据处理装置和/或其他设备以特定方式工作，从而，存储有指令的计算机可读介质则包括一个制造品，其包括实现流程图和/或框图中的一个或多个方框中规定的功能/动作的各个方面的指令。
也可以把计算机可读程序指令加载到计算机、其它可编程数据处理装置、或其它设备上，使得在计算机、其它可编程数据处理装置或其它设备上执行一系列操作步骤，以产生计算机实现的过程，从而使得在计算机、其它可编程数据处理装置、或其它设备上执行的指令实现流程图和/或框图中的一个或多个方框中规定的功能/动作。
附图中的流程图和框图显示了根据本公开的多个实施例的系统、方法和计算机程序产品的可能实现的体系架构、功能和操作。在这点上，流程图或框图中的每个方框可以代表一个模块、程序段或指令的一部分，所述模块、程序段或指令的一部分包含一个或多个用于实现规定的逻辑功能的可执行指令。在有些作为替换的实现中，方框中所标注的功能也可以以不同于附图中所标注的顺序发生。例如，两个连续的方框实际上可以基本并行地执行，它们有时也可以按相反的顺序执行，这依所涉及的功能而定。也要注意的是，框图和/或流程图中的每个方框、以及框图和/或流程图中的方框的组合，可以用执行规定的功能或动作的专用的基于硬件的系统来实现，或者可以用专用硬件与计算机指令的组合来实现。对于本领域技术人员来说公知的是，通过硬件方式实现、通过软件方式实现以及通过软件和硬件结合的方式实现都是等价的。
以上已经描述了本公开的各实施例，上述说明是示例性的，并非穷尽性的，并且也不限于所披露的各实施例。在不偏离所说明的各实施例的范围和精神的情况下，对于本技术领域的普通技术人员来说许多修改和变更都是显而易见的。本文中所用术语的选择，旨在最好地解释各实施例的原理、实际应用或对市场中的技术改进，或者使本技术领域的其它普通技术人员能理解本文披露的各实施例。本公开的范围由所附权利要求来限定。
Claims (10)
1.一种肌电信号的补偿方法，其特征在于，包括：
获取肌电检测设备上设置的至少两个电极所采集的肌电信号；
在第一电极所采集的肌电信号出现数据缺失的情况下，确定其他电极中与所述第一电极相关的第二电极；其中，所述第一电极为所述至少两个电极中的任一个电极，所述其他电极为所述至少两个电极中除所述第一电极以外的电极；
根据所述第二电极所采集的肌电信号，确定所述第一电极所采集的肌电信号中缺失的肌电数据，作为目标肌电数据；
根据所述目标肌电数据对所述第一电极所采集的肌电信号进行补偿。
2.根据权利要求1所述的方法，其特征在于，所述确定其他电极中与所述第一电极相关的第二电极，包括：
根据所述至少两个电极所采集的肌电信号，确定所述第一电极与所述其他电极之间的目标相关系数；
确定与所述第一电极之间的目标相关系数最大的其他电极，作为所述第二电极。
3.根据权利要求2所述的方法，其特征在于，所述根据所述至少两个电极所采集的肌电信号，确定所述第一电极与所述其他电极中之间的目标相关系数，包括：
获取所述第一电极所采集的肌电信号出现数据缺失前的n个第一采样时刻的第一肌电数据；其中，n为正整数；
获取所述第一电极与所述其他电极的采样时刻差；其中，所述采样时刻差指示所述第一电极所采集的肌电信号与所述其他电极所采集的肌电信号中，相匹配的肌电数据的采样时刻之间的差异；
获取每一所述采样时刻差下，所述其他电极所采集的肌电信号中与所述第一肌电数据匹配的其他肌电数据；
根据所述第一肌电数据和每一所述采样时刻差下的所述其他肌电数据，确定所述第一电极与所述其他电极之间的目标相关系数。
4.根据权利要求3所述的方法，其特征在于，所述获取每一所述采样时刻差下，所述其他电极所采集的肌电信号中与所述第一肌电数据匹配的其他肌电数据，包括：
遍历所述采样时刻差；
根据当前遍历的采样时刻差，得到与所述第一采样时刻对应的第二采样时刻；
获取所述其他电极所采集的肌电信号中所述第二采样时刻的肌电数据，作为所述当前遍历的采样时刻差下，与所述第一肌电数据匹配的所述其他肌电数据。
5.根据权利要求3所述的方法，其特征在于，所述根据所述第一肌电数据和每一所述采样时刻差下的所述其他肌电数据，确定所述第一电极与所述其他电极之间的目标相关系数，包括：
根据所述第一肌电数据和每一所述采样时刻差下的所述其他肌电数据，确定所述第一电极与所述其他电极在每一所述采样时刻差下的相关系数；
选取最大的相关系数，作为所述第一电极与所述其他电极之间的目标相关系数。
6.根据权利要求5所述的方法，其特征在于，所述根据所述第二电极所采集的肌电信号，确定所述第一电极所采集的肌电信号中缺失的肌电数据，作为目标肌电数据，包括：
获取所述第一电极采集的肌电信号和所述第二电极所采集的肌电信号中，相匹配的肌电数据之间的映射关系；
获取缺失的肌电数据所对应的所述第一电极的采样时刻，作为第三采样时刻；
根据目标采样时刻差，确定与所述第三采样时刻对应的第四采样时刻；其中，所述目标采样时刻差为所述目标相关系数所对应的采样时刻差；
获取所述第二电极所采集的肌电信号中所述第四采样时刻的肌电数据，作为参考肌电数据；
根据所述参考肌电数据和所述映射关系，得到所述目标肌电数据。
7.根据权利要求6所述的方法，其特征在于，所述获取所述第一电极采集的肌电信号和所述第二电极所采集的肌电信号中，相匹配的肌电数据之间的映射关系，包括：
获取所述目标采样时刻差下，所述第二电极所采集的肌电信号中与所述第一肌电数据匹配的第二肌电数据；
根据所述第一肌电数据和所述第二肌电数据，得到所述映射关系。
8.一种肌电信号的补偿装置，其特征在于，包括：
信号获取模块，用于获取肌电检测设备上设置的至少两个电极所采集的肌电信号；
相关电极确定模块，用于在第一电极所采集的肌电信号出现数据缺失的情况下，确定其他电极中与所述第一电极相关的第二电极；其中，所述第一电极为所述至少两个电极中的任一个电极，所述其他电极为所述至少两个电极中除所述第一电极以外的电极；
缺失数据确定模块，用于根据所述第二电极所采集的肌电信号，确定所述第一电极所采集的肌电信号中缺失的肌电数据，作为目标肌电数据；
信号补偿模块，用于根据所述目标肌电数据对所述第一电极所采集的肌电信号进行补偿。
9.一种肌电检测设备，其特征在于，包括：
如权利要求8所述的装置；或者，
处理器和存储器，所述存储器用于存储指令，所述指令用于控制所述处理器执行根据权利要求1-7中任一项所述的方法。
10.一种可读存储介质，其上存储有计算机程序，所述计算机程序在被处理器执行时实现如权利要求1至7中任一项所述的方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111337718.6A CN114145759B (zh) | 2021-11-11 | 2021-11-11 | 肌电信号的补偿方法、装置、肌电检测设备及存储介质 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111337718.6A CN114145759B (zh) | 2021-11-11 | 2021-11-11 | 肌电信号的补偿方法、装置、肌电检测设备及存储介质 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN114145759A true CN114145759A (zh) | 2022-03-08 |
| CN114145759B CN114145759B (zh) | 2024-05-14 |
Family
ID=80460061
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202111337718.6A Active CN114145759B (zh) | 2021-11-11 | 2021-11-11 | 肌电信号的补偿方法、装置、肌电检测设备及存储介质 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114145759B (zh) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115114962A (zh) * | 2022-07-19 | 2022-09-27 | 歌尔股份有限公司 | 基于表面肌电信号的控制方法、装置及可穿戴设备 |
| CN116196020A (zh) * | 2023-03-24 | 2023-06-02 | 苏州海臻医疗器械有限公司 | 一种基于表面肌电信号的患者肌肉疲劳度评估系统及方法 |
| CN116269450A (zh) * | 2023-03-21 | 2023-06-23 | 苏州海臻医疗器械有限公司 | 一种基于肌电信号的患者肢体康复状态评估系统及方法 |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2003103484A2 (en) * | 2002-06-05 | 2003-12-18 | Nervetrack Ltd. | Method and apparatus for measuring nerve signals in nerve fibers |
| JP2010264168A (ja) * | 2009-05-18 | 2010-11-25 | Seiko Epson Corp | 生体状態測定装置 |
| JP2016054887A (ja) * | 2014-09-09 | 2016-04-21 | 日本電信電話株式会社 | 欠損生体信号推定方法 |
| CN112912003A (zh) * | 2018-10-22 | 2021-06-04 | 日本光电工业株式会社 | 心电图的缺失导联的补偿 |
| CN113555132A (zh) * | 2020-04-24 | 2021-10-26 | 华为技术有限公司 | 多源数据处理方法、电子设备和计算机可读存储介质 |
-
2021
- 2021-11-11 CN CN202111337718.6A patent/CN114145759B/zh active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2003103484A2 (en) * | 2002-06-05 | 2003-12-18 | Nervetrack Ltd. | Method and apparatus for measuring nerve signals in nerve fibers |
| JP2010264168A (ja) * | 2009-05-18 | 2010-11-25 | Seiko Epson Corp | 生体状態測定装置 |
| JP2016054887A (ja) * | 2014-09-09 | 2016-04-21 | 日本電信電話株式会社 | 欠損生体信号推定方法 |
| CN112912003A (zh) * | 2018-10-22 | 2021-06-04 | 日本光电工业株式会社 | 心电图的缺失导联的补偿 |
| CN113555132A (zh) * | 2020-04-24 | 2021-10-26 | 华为技术有限公司 | 多源数据处理方法、电子设备和计算机可读存储介质 |
| WO2021213165A1 (zh) * | 2020-04-24 | 2021-10-28 | 华为技术有限公司 | 多源数据处理方法、电子设备和计算机可读存储介质 |
Non-Patent Citations (2)
| Title |
|---|
| 丁其川;赵新刚;韩建达;: "基于肌电信号容错分类的手部动作识别", 机器人, no. 01, 15 January 2015 (2015-01-15) * |
| 高峥;徐震;: "基于多元回归KNN的油田缺失数据填充方法", 信息技术, no. 04, 16 April 2020 (2020-04-16) * |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115114962A (zh) * | 2022-07-19 | 2022-09-27 | 歌尔股份有限公司 | 基于表面肌电信号的控制方法、装置及可穿戴设备 |
| CN116269450A (zh) * | 2023-03-21 | 2023-06-23 | 苏州海臻医疗器械有限公司 | 一种基于肌电信号的患者肢体康复状态评估系统及方法 |
| CN116269450B (zh) * | 2023-03-21 | 2023-12-19 | 苏州海臻医疗器械有限公司 | 一种基于肌电信号的患者肢体康复状态评估系统及方法 |
| CN116196020A (zh) * | 2023-03-24 | 2023-06-02 | 苏州海臻医疗器械有限公司 | 一种基于表面肌电信号的患者肌肉疲劳度评估系统及方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN114145759B (zh) | 2024-05-14 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN114145759A (zh) | 肌电信号的补偿方法、装置、肌电检测设备及存储介质 | |
| CN106897688B (zh) | 交互式投影装置、控制交互式投影的方法和可读存储介质 | |
| CN108352459B (zh) | 显示器和包括显示器的电子设备 | |
| CN109491586B (zh) | 虚拟对象控制方法及装置、电子设备、存储介质 | |
| CN103519794A (zh) | 测量装置、测量方法、程序、存储介质和测量系统 | |
| EP2947543A1 (en) | Method for displaying battery information of electronic device and electronic device performing the same | |
| US8253706B2 (en) | Apparatus using a differential analog-to-digital converter | |
| KR20180136242A (ko) | 생체 센서의 전극을 이용한 충전 방법 및 이를 적용한 전자 장치 | |
| CN107810466B (zh) | 触摸传感器中的多个相关 | |
| CN109379678A (zh) | 非线性补偿方法、装置、存储介质及终端设备 | |
| US10585529B2 (en) | Selecting correlation reference based on noise estimation | |
| WO2020038453A1 (zh) | 用于脉搏检测的可穿戴设备以及脉搏检测方法 | |
| EP3971694B1 (en) | Noise measurement circuit, self-capacitance measurement method, touch chip and electronic device | |
| US11941175B2 (en) | Skin-to-skin contact detection | |
| CN111225603A (zh) | 提供与用户活动相对应的压力指数的电子设备和方法 | |
| US20130305088A1 (en) | Electronic device and test card thereof | |
| US20170177144A1 (en) | Touch display device and touch display method | |
| Vizitiu et al. | Exhaustive description of the system architecture and prototype implementation of an IoT-based eHealth biometric monitoring system for elders in independent living | |
| US20200409470A1 (en) | Wearable gesture recognition device and associated operation method and system | |
| US20170082695A1 (en) | Method and apparatus for determining battery parameters | |
| US9607536B2 (en) | Testing apparatus for electronic device | |
| CN113467647A (zh) | 皮肤与皮肤接触检测 | |
| CN105829998A (zh) | 将装置绑定到计算设备 | |
| CN111035393A (zh) | 三维步态数据处理方法、系统、服务器及存储介质 | |
| CN107450774B (zh) | 触控检测方法、电路、存储介质、处理器和终端 |
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 |