CN108143414A - It is a kind of it is adaptive filter out surface electromyogram signal electro photoluminescence artefact in line method - Google Patents

It is a kind of it is adaptive filter out surface electromyogram signal electro photoluminescence artefact in line method Download PDF

Info

Publication number
CN108143414A
CN108143414A CN201810027796.8A CN201810027796A CN108143414A CN 108143414 A CN108143414 A CN 108143414A CN 201810027796 A CN201810027796 A CN 201810027796A CN 108143414 A CN108143414 A CN 108143414A
Authority
CN
China
Prior art keywords
artefact
electro photoluminescence
signal
pulse
discharge portion
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
Application number
CN201810027796.8A
Other languages
Chinese (zh)
Other versions
CN108143414B (en
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.)
Fuzhou University
Original Assignee
Fuzhou University
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 Fuzhou University filed Critical Fuzhou University
Priority to CN201810027796.8A priority Critical patent/CN108143414B/en
Publication of CN108143414A publication Critical patent/CN108143414A/en
Application granted granted Critical
Publication of CN108143414B publication Critical patent/CN108143414B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/24Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
    • A61B5/316Modalities, i.e. specific diagnostic methods
    • A61B5/389Electromyography [EMG]
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/24Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
    • A61B5/316Modalities, i.e. specific diagnostic methods
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/72Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7203Signal processing specially adapted for physiological signals or for diagnostic purposes for noise prevention, reduction or removal
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/72Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7235Details of waveform analysis
    • A61B5/725Details of waveform analysis using specific filters therefor, e.g. Kalman or adaptive filters

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • Veterinary Medicine (AREA)
  • Public Health (AREA)
  • Physics & Mathematics (AREA)
  • General Health & Medical Sciences (AREA)
  • Biophysics (AREA)
  • Pathology (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Medical Informatics (AREA)
  • Molecular Biology (AREA)
  • Signal Processing (AREA)
  • Artificial Intelligence (AREA)
  • Physiology (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Psychiatry (AREA)
  • Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
  • Electrotherapy Devices (AREA)

Abstract

The invention discloses it is a kind of it is adaptive filter out surface electromyogram signal electro photoluminescence artefact in line method, including step:1)Pure electro photoluminescence artefact is acquired using subthreshold value stimulation, line modeling is carried out to the discharge portion of pure electro photoluminescence artefact, obtains filtering mould;2)To the electromyography signal that each electro photoluminescence generates, using the pulse signal of electro photoluminescence as the initial signal of wave filter, mask filter is used in the stage pulse of electro photoluminescence artefact, is shielded to the negative-going pulse spike of electro photoluminescence artefact;3)After shielding, artefact estimated value is calculated, and the artefact estimated value is deducted from the electromyography signal of acquisition using Filtering Template, artefact discharge portion is filtered out, until next electric stimulation pulse starts.The present invention can effectively filter out electro photoluminescence artefact using mixed filtering method, remove the electro photoluminescence artefact in electromyography signal in real time, and will not be distorted muscle electric signal.

Description

It is a kind of it is adaptive filter out surface electromyogram signal electro photoluminescence artefact in line method
Technical field
The present invention relates to electrical stimulation technology fields, and in particular to one kind adaptively filters out surface electromyogram signal electro photoluminescence artefact In line method.
Background technology
When electric stimulation and electromyographic signal collection system are simultaneously in use, electro photoluminescence can seriously pollute collected surface Electromyography signal (electromyography, EMG).The electromyography signal acquired under electro photoluminescence mainly includes muscle electric signal (M Wave), independent desire electromyography signal and electro photoluminescence artefact.The Amplitude Ration surface electromyogram signal amplitude of stimulus artifact is much higher, and With M waves and independent desire electromyography signal aliasing over the frequency domain, electromyography signal seriously polluted.Therefore there are electro photoluminescence feelings Under condition, analysis and use to electromyography signal should carry out filtering out for electro photoluminescence artefact first.
At present in the literature research delivered, the method that filters out electro photoluminescence artefact in surface electromyogram signal mainly has following Several method.(1)The shielding mode kept using sampling, i.e., when exporting electro photoluminescence, stop the acquisition of EMG signal.And when thorn When sharp frequency is higher, shielded signal is more, has seriously affected the quality of EMG signal.The method of shielding can delete it is excessive or A part of stimulus artifact is left, this depends on the length of shielding window.When stimulus is constant-current source, the back segment of electro photoluminescence artefact Response is slow, and electro photoluminescence can cover longer time region for artefact waveform.(2)Use experience mode decomposition, wavelet decomposition etc. Method decomposes signal, and the signal difference component after decomposition is filtered.Method based on signal decomposition, such as Empirical mode decomposition can decomposite intrinsic mode functions component, although being to decompose different frequency contents from trend, Can illustrate that some eigencomponent is electro photoluminescence artefact there is no index, wavelet decomposition there are it is similary the problem of.And based on decomposition Method is not suitable for application on site, it is impossible to be used in real-time system.(3)Join using the method for artefact template or using having Examine the adaptive filter method of signal.The validity of this method depends on the accuracy of template and the adaptability of template, needs Study the relationship of position of electrical stimulation waveforms feature, stimulating electrode and EMG sensors etc. and electro photoluminescence artefact.Have been reported explanation The position of stimulating electrode and EMG sensors is affected to electrical stimulation waveforms feature, and institute is in this way for testing installation position It puts more demanding.And when using offline template method, the change of template can be led to when experiment condition changes so that occur Noise can not filter out or useful signal is also weakened while noise filtering.
Invention content
In view of the deficiencies of the prior art, the present invention propose it is a kind of it is adaptive filter out surface electromyogram signal electro photoluminescence artefact Line method can effectively filter out electro photoluminescence artefact, and not destroy muscle electric signal;Electro photoluminescence artefact and muscle electric signal are mixed Folded degree does not have particular/special requirement;The template used is adaptable, can adapt to the change of different experimental subjects and electrode position Change;The online filtering synchronous with electro photoluminescence is realized, available for being used in real-time system;For having other physiology of similar form Electro photoluminescence artefact in electric signal can be filtered using the same method.
To achieve the above object, the technical scheme is that:It is a kind of adaptively to filter out surface electromyogram signal electro photoluminescence puppet Mark in line method, including:
Step S1:Pure electro photoluminescence artefact is acquired using subthreshold value stimulation, the discharge portion of pure electro photoluminescence artefact is built online Mould obtains Filtering Template;
Step S2:To the electromyography signal that each electro photoluminescence generates, believe using the pulse signal of electro photoluminescence as the starting of wave filter Number, mask filter is used in the stage pulse of electro photoluminescence artefact, is shielded to the negative-going pulse spike of electro photoluminescence artefact;
Step S3:After shielding, artefact estimated value is calculated, and the puppet is deducted from the electromyography signal of acquisition using Filtering Template Mark estimated value filters out artefact discharge portion, until next electric stimulation pulse starts.
Further, the modeling process in the step S1 is specially:
The discharge portion of pure electro photoluminescence artefact is modeled using autoregressive process model, is distinguished using recurrent least square method Know,dRank autoregressive process model is as follows:
Wherein,dIt is the exponent number of autoregressive process model,y i For the observation of time series,w i To estimate weights,e k For white Gaussian Noise process.
Compared with prior art, the present invention has advantageous effect:
(1)Realize the online filtering synchronous with electro photoluminescence, can be used for real-time system, to the electro photoluminescence artefact in each period into Row filtering, mixed filtering method can effectively filter out electro photoluminescence artefact, and will not be distorted muscle electric signal;
(2)Template modeling is carried out using the surface electromyogram signal of subthreshold value electro photoluminescence in the pre-stimulation stage, enables the template of foundation It enough adapts to experimental subjects and electrode puts up the otherness of position;
(3)The filtering shielding used has adaptive length of window, without manually setting parameter;
(4)The present invention does not have particular/special requirement to the aliasing degree of electro photoluminescence artefact and muscle electric signal, suitable for more generally feelings Condition.
Description of the drawings
Fig. 1 is that a kind of adaptive flow in line method for filtering out surface electromyogram signal electro photoluminescence artefact of the present invention is illustrated Figure;
Fig. 2 is the oscillogram of one embodiment of the invention electro photoluminescence artefact.
Specific embodiment
The present invention will be further described with reference to the accompanying drawings and embodiments.
As shown in Figure 1, the present invention it is a kind of it is adaptive filter out surface electromyogram signal electro photoluminescence artefact in line method, first Step is pre-stimulation, and pure electro photoluminescence artefact is acquired using subthreshold value stimulation, the discharge portion of electro photoluminescence artefact is carried out online Modeling, second step is in filtering stage, and the spike stage of electro photoluminescence artefact is filtered using screen method, to electro photoluminescence puppet The discharge portion of mark uses template convolution.
First step is pre-stimulation, in the pre-stimulation stage, acquires pure electro photoluminescence artefact using subthreshold value electro photoluminescence and believes Number for establishing artefact template.Subthreshold value stimulation refers to the electro photoluminescence less than human motion unit activating threshold value.Subthreshold value stimulates not Action potential can be generated, collected surface electromyogram signal is pure electro photoluminescence artefact.The waveform of electro photoluminescence artefact such as Fig. 2 institutes Show, the dotted line left side is known as the stage pulse of electro photoluminescence artefact, and the discharge regime of electro photoluminescence artefact is known as on the right side of dotted line.Utilize acquisition Purified signal as observation, establish the model of electro photoluminescence artefact discharge portion as follow-up Filtering Template.Utilize autoregression Model models electro photoluminescence artefact discharge portion using autoregression (autoregressive, AR) process model, using passing Least square method (recursive least squares, RLS) is returned to recognize.dRank autoregression AR models have following form:
WhereindIt is the exponent number of AR models,y i For the observation of time series,w i To estimate weights,e k For white Gaussian noise process.
Second step is the filtering of electro photoluminescence artefact.Filtering is divided into two parts, i.e. mask filter and template convolution, such as Shown in Fig. 1.Because under different intensity of electric stimulus, the stage pulse of the electro photoluminescence artefact signal in each period differs greatly and discharges Stage has similitude, so the stage pulse to electro photoluminescence artefact uses mask filter, to the discharge portion of electro photoluminescence artefact Use template convolution.In filtering stage, filtering is, rising using the pulse signal of electro photoluminescence as wave filter synchronous with electro photoluminescence Beginning signal is first begin to shield, and shields to the negative-going pulse spike of electro photoluminescence artefact.After shielding, artefact electric discharge rank is utilized The Filtering Template of section calculates, and is deducted from acquisition signal and count counted artefact estimated value, artefact discharge portion is filtered out, under One electric stimulation pulse starts, and completes the electro photoluminescence artefact filtering of a cycle.
Particular embodiments described above elaborates the purpose of the present invention, technical solution and achievement, is answered Understand, the above is only a specific embodiment of the present invention, is not intended to restrict the invention, all essences in the present invention God and any modification, equivalent substitution, improvement and etc. within principle, done, should all be included in the protection scope of the present invention.

Claims (2)

1. it is a kind of it is adaptive filter out surface electromyogram signal electro photoluminescence artefact in line method, which is characterized in that including:
Step S1:Pure electro photoluminescence artefact is acquired using subthreshold value stimulation, the discharge portion of pure electro photoluminescence artefact is built online Mould obtains Filtering Template;
Step S2:To the electromyography signal that each electro photoluminescence generates, believe using the pulse signal of electro photoluminescence as the starting of wave filter Number, mask filter is used in the stage pulse of electro photoluminescence artefact, is shielded to the negative-going pulse spike of electro photoluminescence artefact;
Step S3:After shielding, artefact estimated value is calculated, and the puppet is deducted from the electromyography signal of acquisition using Filtering Template Mark estimated value filters out artefact discharge portion, until next electric stimulation pulse starts.
2. it is according to claim 1 it is adaptive filter out surface electromyogram signal electro photoluminescence artefact in line method, feature exists In the modeling process in the step S1 is specially:
The discharge portion of pure electro photoluminescence artefact is modeled using autoregressive process model, is distinguished using recurrent least square method Know,dRank autoregressive process model is as follows:
Wherein,dIt is the exponent number of autoregressive process model,y i For the observation of time series,w i To estimate weights,e k For white Gaussian Noise process.
CN201810027796.8A 2018-01-11 2018-01-11 Online method for adaptively filtering surface electromyographic signal electrical stimulation artifact Active CN108143414B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810027796.8A CN108143414B (en) 2018-01-11 2018-01-11 Online method for adaptively filtering surface electromyographic signal electrical stimulation artifact

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810027796.8A CN108143414B (en) 2018-01-11 2018-01-11 Online method for adaptively filtering surface electromyographic signal electrical stimulation artifact

Publications (2)

Publication Number Publication Date
CN108143414A true CN108143414A (en) 2018-06-12
CN108143414B CN108143414B (en) 2021-05-18

Family

ID=62461377

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810027796.8A Active CN108143414B (en) 2018-01-11 2018-01-11 Online method for adaptively filtering surface electromyographic signal electrical stimulation artifact

Country Status (1)

Country Link
CN (1) CN108143414B (en)

Citations (18)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO1996025093A1 (en) * 1995-02-17 1996-08-22 Ep Technologies, Inc. Systems and methods for filtering signals derived from biological events
CN2530301Y (en) * 2001-03-09 2003-01-08 北京卡迪欧医疗设备有限责任公司 Analyser for electrocardiosignal
JP2004240518A (en) * 2003-02-03 2004-08-26 Nippon Telegr & Teleph Corp <Ntt> Signal analyzing method and device, and signal analyzing program and recording medium recording the program
US20070118047A1 (en) * 2005-10-20 2007-05-24 Brian Tracey Automated stimulus artifact removal for nerve conduction studies
CN101317794A (en) * 2008-03-11 2008-12-10 清华大学 Myoelectric control ability detecting and training method for hand-prosthesis with multiple fingers and multiple degrees of freedom
US20100068751A1 (en) * 2008-09-18 2010-03-18 Imec Method and System for Artifact Reduction
CN102164537A (en) * 2008-09-17 2011-08-24 Med-El电气医疗器械有限公司 Stimulus artifact removal for neuronal recordings
CN102697496A (en) * 2012-06-07 2012-10-03 天津大学 Filtering method for functional electrical stimulation surface electromyogram signal
US8521268B2 (en) * 2011-05-10 2013-08-27 Medtronic, Inc. Techniques for determining cardiac cycle morphology
CN203506713U (en) * 2013-09-25 2014-04-02 中国科学院昆明动物研究所 Digitalized multi-channel electric stimulus artifact eliminating device
US20150087934A1 (en) * 2000-04-05 2015-03-26 Neuropace, Inc. Devices and methods for monitoring physiological information relating to sleep with an implantable device
WO2015065599A1 (en) * 2013-10-28 2015-05-07 Medtronic, Inc. Devices for sensing physiological signals during stimulation therapy
US9116835B1 (en) * 2014-09-29 2015-08-25 The United States Of America As Represented By The Secretary Of The Army Method and apparatus for estimating cerebral cortical source activations from electroencephalograms
CN105031812A (en) * 2015-06-09 2015-11-11 电子科技大学 Functional electrostimulation closed-loop control system and method of electromyographic signal feedback
CN105122273A (en) * 2013-01-17 2015-12-02 科迪影技术股份有限公司 Non-local mean filtering for electrophysiological signals
CN105873506A (en) * 2013-11-07 2016-08-17 赛佛欧普手术有限公司 Systems and methods for detecting nerve function
CN105997064A (en) * 2016-05-17 2016-10-12 成都奥特为科技有限公司 Method for identifying human lower limb surface EMG signals (electromyographic signals)
CN107233094A (en) * 2016-03-29 2017-10-10 上海海神医疗电子仪器有限公司 A kind of electromyogram evoked potentuial measuring system

Patent Citations (18)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO1996025093A1 (en) * 1995-02-17 1996-08-22 Ep Technologies, Inc. Systems and methods for filtering signals derived from biological events
US20150087934A1 (en) * 2000-04-05 2015-03-26 Neuropace, Inc. Devices and methods for monitoring physiological information relating to sleep with an implantable device
CN2530301Y (en) * 2001-03-09 2003-01-08 北京卡迪欧医疗设备有限责任公司 Analyser for electrocardiosignal
JP2004240518A (en) * 2003-02-03 2004-08-26 Nippon Telegr & Teleph Corp <Ntt> Signal analyzing method and device, and signal analyzing program and recording medium recording the program
US20070118047A1 (en) * 2005-10-20 2007-05-24 Brian Tracey Automated stimulus artifact removal for nerve conduction studies
CN101317794A (en) * 2008-03-11 2008-12-10 清华大学 Myoelectric control ability detecting and training method for hand-prosthesis with multiple fingers and multiple degrees of freedom
CN102164537A (en) * 2008-09-17 2011-08-24 Med-El电气医疗器械有限公司 Stimulus artifact removal for neuronal recordings
US20100068751A1 (en) * 2008-09-18 2010-03-18 Imec Method and System for Artifact Reduction
US8521268B2 (en) * 2011-05-10 2013-08-27 Medtronic, Inc. Techniques for determining cardiac cycle morphology
CN102697496A (en) * 2012-06-07 2012-10-03 天津大学 Filtering method for functional electrical stimulation surface electromyogram signal
CN105122273A (en) * 2013-01-17 2015-12-02 科迪影技术股份有限公司 Non-local mean filtering for electrophysiological signals
CN203506713U (en) * 2013-09-25 2014-04-02 中国科学院昆明动物研究所 Digitalized multi-channel electric stimulus artifact eliminating device
WO2015065599A1 (en) * 2013-10-28 2015-05-07 Medtronic, Inc. Devices for sensing physiological signals during stimulation therapy
CN105873506A (en) * 2013-11-07 2016-08-17 赛佛欧普手术有限公司 Systems and methods for detecting nerve function
US9116835B1 (en) * 2014-09-29 2015-08-25 The United States Of America As Represented By The Secretary Of The Army Method and apparatus for estimating cerebral cortical source activations from electroencephalograms
CN105031812A (en) * 2015-06-09 2015-11-11 电子科技大学 Functional electrostimulation closed-loop control system and method of electromyographic signal feedback
CN107233094A (en) * 2016-03-29 2017-10-10 上海海神医疗电子仪器有限公司 A kind of electromyogram evoked potentuial measuring system
CN105997064A (en) * 2016-05-17 2016-10-12 成都奥特为科技有限公司 Method for identifying human lower limb surface EMG signals (electromyographic signals)

Non-Patent Citations (6)

* Cited by examiner, † Cited by third party
Title
AKIRA FURUI等: "An artificial EMG generation model based on signal-dependent noise and related application to motion classification", 《PLOS ONE》 *
DEREK T. O’KEEFFE等: "Stimulus artifact removal using a software-based two-stage peak detection algorithm", 《JOURNAL OF NEUROSCIENCE METHODS》 *
GERT PFURTSCHELLER等: "Separability of EEG Signals Recorded During Right and Left Motor Imagery Using Adaptive Autoregressive Parameters", 《IEEE TRANSACTIONS ON REHABILITATION ENGINEERING》 *
XIN, YI, JUN, ET AL.: "A Blink Restoration System With Contralateral EMG Triggered Stimulation and Real-Time Artifact Blanking", 《IEEE TRANSACTIONS ON BIOMEDICAL CIRCUITS AND SYSTEMS》 *
宫琴, 叶大田, 郭连生,等: "预测器—减法器—恢复器构成的滤波器去除刺激伪迹的研究", 《声学学报》 *
杨基海 等: "利用自组织竞争神经网络提取NMEG信号的MUAP模板", 《生物医学工程学杂志》 *

Also Published As

Publication number Publication date
CN108143414B (en) 2021-05-18

Similar Documents

Publication Publication Date Title
DE3732699C2 (en) Implantable pacemaker
CN103845137B (en) Based on the robot control method of stable state vision inducting brain-machine interface
CN108742611B (en) Autonomous myoelectricity extraction method under random dynamic parameter electrical stimulation condition
Yi et al. A blink restoration system with contralateral EMG triggered stimulation and real-time artifact blanking
Nimunkar et al. EMD-based 60-Hz noise filtering of the ECG
CN108294748A (en) A kind of eeg signal acquisition and sorting technique based on stable state vision inducting
DE69940071D1 (en) DEVICE FOR NOISE REDUCTION OF A FOLUTE ECG SIGNAL
CN107095672B (en) Method for removing artifact signals in electroencephalogram signals
CN108143414A (en) It is a kind of it is adaptive filter out surface electromyogram signal electro photoluminescence artefact in line method
Cline et al. Advanced artifact removal for automated TMS-EEG data processing
DE4106856A1 (en) METHOD AND DEVICE FOR FILTERING BASELINE FLUCTUATIONS FROM AN ELECTROCARDIOGRAM
Hu et al. ECG cancellation for surface electromyography measurement using independent component analysis
CN205850003U (en) A kind of SCM Based myoelectricity boost pulse instrument
Metcalfe et al. An enhancement to velocity selective discrimination of neural recordings: Extraction of neuronal firing rates
CN113208622A (en) Electroencephalogram EEG signal denoising method and system based on deep neural network technology
DE2717530A1 (en) Maternal heart signal noise removal for ECG signals - by storing shape of interference pulse for subtraction from subsequent pulses
McKenzie et al. Same-electrode stimulation and recording with dynamic hardware artefact suppression
CN114052752B (en) Method for filtering interference of multi-channel surface electromyographic signal center electric QRS complex
Andrews et al. Artifact Adaptive Ideal Filtering of EMG Signals Contaminated by Spinal Cord Transcutaneous Stimulation
Yochum et al. EMG artifacts removal during electrical stimulation, a CWT based technique
Medithe et al. Removal of ocular artifacts in EEG
CN115607160A (en) Method and system for detecting induced bioelectric signals under magnetoelectric stimulation
Gath et al. Advances in processing and pattern analysis of biological signals
CN107296607A (en) Muscular recuperation condition detection method
Choudhary et al. ECG Signal denoising using UWT and features extraction using LabVIEW

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