CN112494008A - 基于ppg信号的呼吸率测量方法及装置 - Google Patents
基于ppg信号的呼吸率测量方法及装置 Download PDFInfo
- Publication number
- CN112494008A CN112494008A CN202011180742.9A CN202011180742A CN112494008A CN 112494008 A CN112494008 A CN 112494008A CN 202011180742 A CN202011180742 A CN 202011180742A CN 112494008 A CN112494008 A CN 112494008A
- Authority
- CN
- China
- Prior art keywords
- ppg signal
- respiration rate
- signal
- ppg
- processing
- 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.)
- Pending
Links
- 230000029058 respiratory gaseous exchange Effects 0.000 title claims abstract description 45
- 238000000034 method Methods 0.000 title claims abstract description 26
- 238000012545 processing Methods 0.000 claims abstract description 32
- 238000005259 measurement Methods 0.000 claims abstract description 15
- 238000001914 filtration Methods 0.000 claims abstract description 11
- 238000000691 measurement method Methods 0.000 claims abstract description 11
- 230000036387 respiratory rate Effects 0.000 claims abstract description 9
- 230000003068 static effect Effects 0.000 claims description 7
- 239000000126 substance Substances 0.000 claims description 7
- 238000012935 Averaging Methods 0.000 claims description 5
- 238000012886 linear function Methods 0.000 claims description 3
- 238000007619 statistical method Methods 0.000 claims description 3
- 238000010586 diagram Methods 0.000 description 8
- 238000001514 detection method Methods 0.000 description 3
- CURLTUGMZLYLDI-UHFFFAOYSA-N Carbon dioxide Chemical compound O=C=O CURLTUGMZLYLDI-UHFFFAOYSA-N 0.000 description 2
- 210000002565 arteriole Anatomy 0.000 description 1
- 210000003403 autonomic nervous system Anatomy 0.000 description 1
- 239000008280 blood Substances 0.000 description 1
- 210000004369 blood Anatomy 0.000 description 1
- 230000036772 blood pressure Effects 0.000 description 1
- 230000036760 body temperature Effects 0.000 description 1
- 238000009529 body temperature measurement Methods 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 229910002092 carbon dioxide Inorganic materials 0.000 description 1
- 239000001569 carbon dioxide Substances 0.000 description 1
- 239000002131 composite material Substances 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000001939 inductive effect Effects 0.000 description 1
- 238000013186 photoplethysmography Methods 0.000 description 1
- 230000000241 respiratory effect Effects 0.000 description 1
- 238000009531 respiratory rate measurement Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 210000000115 thoracic cavity Anatomy 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/0205—Simultaneously evaluating both cardiovascular conditions and different types of body conditions, e.g. heart and respiratory condition
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/103—Detecting, measuring or recording devices for testing the shape, pattern, colour, size or movement of the body or parts thereof, for diagnostic purposes
- A61B5/11—Measuring movement of the entire body or parts thereof, e.g. head or hand tremor, mobility of a limb
-
- 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/7225—Details of analog processing, e.g. isolation amplifier, gain or sensitivity adjustment, filtering, baseline or drift compensation
-
- 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
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Physiology (AREA)
- Molecular Biology (AREA)
- General Health & Medical Sciences (AREA)
- Biophysics (AREA)
- Pathology (AREA)
- Signal Processing (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Medical Informatics (AREA)
- Veterinary Medicine (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- Physics & Mathematics (AREA)
- Public Health (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Psychiatry (AREA)
- Cardiology (AREA)
- Dentistry (AREA)
- Oral & Maxillofacial Surgery (AREA)
- Power Engineering (AREA)
- Pulmonology (AREA)
- Measurement Of The Respiration, Hearing Ability, Form, And Blood Characteristics Of Living Organisms (AREA)
Abstract
本发明提供一种基于PPG信号的呼吸率测量方法及装置，该方法包括：实时采集PPG信号；对所述PPG信号进行第一处理，以得到第一处理信号；所述第一处理包括滤波处理和差分处理；对所述第一处理信号进行第二处理，以得到上包络线；根据所述上包络线计算得到最终的呼吸率。此外，还提供了基于PPG信号的呼吸率测量装置。本发明提供的一种基于PPG信号的呼吸率测量方法及装置，通过基于PPG的信号实现呼吸率测量，相对于现有的呼吸测量方法受外界的影响小，使用场景广，操作简单易携带，并且测量更快速。
Description
技术领域
本发明涉及PPG信号的呼吸率测量技术领域，尤其涉及一种基于PPG信号的呼吸率测量方法及装置。
背景技术
现有测量呼吸的方法有胸廓阻抗、电感体积描记、温度测量、气流检测、湿度检测、二氧化碳水平测量等方法，但现有的测试方法中均存在测试设备复杂、测量受外界条件影响大、设备不能携带、测量慢等缺点。
发明内容
本发明主要目的是提供一种基于PPG信号的呼吸率测量方法及装置，旨在解决现有呼吸测量的问题。
为实现上述目的，本发明提供了一种基于PPG信号的呼吸率测量方法，所述基于PPG信号的呼吸率测量方法包括：
步骤S10：实时采集PPG信号；
步骤S20：对所述PPG信号进行第一处理，以得到第一处理信号；所述第一处理包括滤波处理和差分处理；
步骤S30：对所述第一处理信号进行第二处理，以得到上包络线；
步骤S40：根据所述上包络线计算得到最终的呼吸率。
进一步地，在步骤S10之前还包括：
步骤S1：收集用户的加速度计信号以得到用户运动状态；
步骤S2：判断所述用户运动状态是否为静止状态，如果不是静止状态则执行步骤S3，如果是则执行步骤S4；
步骤S3：提醒用户平静并保持至静止状态，返回步骤S1；
步骤S4：在用户处于静止状态时实时采集PPG信号。
进一步地，所述PPG信号为16Hz的PPG绿光信号。
进一步地，所述步骤S20包括：
步骤S210：对所述PPG信号使用一个四阶的平均滤波器和一个四阶的巴特沃兹带通滤波器进行滤波，得到预处理信号；
步骤S220：对所述预处理信号进行一阶差分处理，得到第一处理信号。
进一步地，所述平均滤波器为：
其中，
X(n)为输入信号，即所述PPG信号；
Y(n)为输出信号；
K和N为自然数。
进一步地，所述巴特沃兹带通滤波器为：
其中，
aK和bK为选择所需要频率信号计算得到的滤波器的系数。
进一步地，所述一阶差分为：
Y(n)＝X(n)-X(n-1)。
进一步地，所述步骤S30包括：
步骤S310：根据预设的第一阈值找出所述第一处理信号中的极大值点；
步骤S320：基于上包络线函数计算三次样条插值得到上包络线；所述上包络线函数在所求区间内[x1,xn]是二次连续可微函数，其中的每个小区间[xk,xk+1](k＝1,2,…,n-1)是一个不超过三次的多项式，其二阶导数存在且为线性函数，一阶导数连续。
进一步地，所述步骤S40包括：
步骤S410：根据预设的第二阈值找出所述上包络线指定区间的最大值点；
步骤S420：根据所述最大值点计算所述上包络线每一点的斜率；
步骤S430：使用所述斜率辅助对所述指定区间的所有呼吸率进行统计分析，计算得到最终的呼吸率。
同时，本发明还提供一种基于PPG信号的呼吸率测量装置，所述基于PPG信号的呼吸率测量装置包括如上所述的基于PPG信号的呼吸率测量方法。
本发明提供的一种基于PPG信号的呼吸率测量方法及装置，通过基于PPG的信号实现呼吸率测量，相对于现有的呼吸测量方法受外界的影响小，使用场景广，操作简单易携带，并且测量更快速。
附图说明
图1为本发明一实施例提供的基于PPG信号的呼吸率测量方法的流程示意图；
图2为本发明一实施例提供的用户运动状态检测方法的流程示意图；
图3为图1中步骤S20的流程示意图；
图4为图1中步骤S30的流程示意图；
图5为图1中步骤S40的流程示意图；
图6为本发明一实施例提供的PPG信号的示意图；
图7为本发明一实施例提供的预处理信号的示意图；
图8为本发明一实施例提供的第一处理信号的示意图；
图9为本发明一实施例提供的上包络线的示意图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅用以解释本发明，并不用于限定本发明。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
本发明一实施例提供一种基于PPG信号的呼吸率测量方法，所述PPG是光电容积脉搏波描记法(PhothPlethysmoGraphy，PPG)，基于PPG信号的呼吸率测量方法是一种基于PPG信号进行呼吸率检测的方法，通过检测微动脉、毛细血管中血液容积的变化来获得容积脉搏波，PPG信号是一种复合信号，主要由分布在0.007Hz-1.5Hz的5个频率组成，PPG信号和呼吸、血压、心率、体温和自主神经系统有关，因此，本发明利用PPG信号找到所蕴含的呼吸信号来进行呼吸率的计算。
请参阅图1，本发明的一实施例提供一种基于PPG信号的呼吸率测量方法，所述基于PPG信号的呼吸率测量方法包括：
步骤S10：实时采集PPG信号；
步骤S20：对所述PPG信号进行第一处理，以得到第一处理信号；所述第一处理包括滤波处理和差分处理；
步骤S30：对所述第一处理信号进行第二处理，以得到上包络线；
步骤S40：根据所述上包络线计算得到最终的呼吸率。
请参阅图2，具体在测量时，需要让用户保持静止的状态，因此，在进行所述步骤S10进行实时采集PPG信号前，先进行用户运动状态检测，其步骤包括：
步骤S1：收集用户的加速度计信号以得到用户运动状态；实际测量时，用户可佩戴智能穿戴设备如智能手表实现用户的加速度计信号的收集；
步骤S2：判断所述用户运动状态是否为静止状态，如果不是静止状态则执行步骤S3，如果是则执行步骤S4；
步骤S3：提醒用户平静并保持至静止状态，返回步骤S1；
步骤S4：在用户处于静止状态时实时采集PPG信号。
确保用户处于静止状态后再进行PPG信号的采集，所述PPG信号为16Hz的PPG绿光信号。
请参阅图3，对实时采集到的PPG信号进行第一处理，具体地，所述步骤S20包括：
步骤S210：对所述PPG信号使用一个四阶的平均滤波器和一个四阶的巴特沃兹带通滤波器进行滤波，得到预处理信号；其中，
所述平均滤波器为：
其中，
X(n)为输入信号，即所述PPG信号；
Y(n)为输出信号；
K和N为自然数。
用某一点邻域内几个值的均值来代替该点的值，计算量小，处理速度快，对高频噪声的滤除效果很好。
所述巴特沃兹带通滤波器为：
其中，
aK和bK为选择所需要频率信号计算得到的滤波器的系数。
巴特沃兹滤波器的时域等效差分方程，滤除高频信号和基线漂移干扰。
步骤S220：对所述预处理信号进行一阶差分处理，得到第一处理信号。
所述一阶差分为：
Y(n)＝X(n)-X(n-1)。
对实时采集的PPG信号用一个四阶的平均滤波器和一个四阶的巴特沃兹带通滤波器进行滤波，滤除较高频率的信号干扰和基线漂移产生的影响。由于呼吸信号频率较低，巴特沃兹带通滤波器不能完全消除基线漂移的干扰，因此对得到的滤波信号再进行一阶的差分就完成了对PPG信号的原始数据的第一处理，得到第一处理信号。
进一步地，对所述第一处理信号进行第二处理，以得到上包络线。
请参阅图4，所述步骤S30包括：
步骤S310：根据预设的第一阈值找出所述第一处理信号中的极大值点；
步骤S320：基于上包络线函数计算三次样条插值得到上包络线；所述上包络线函数在所求区间内[x1,xn]是二次连续可微函数，其中的每个小区间[xk,xk+1](k＝1,2,…,n-1)是一个不超过三次的多项式，其二阶导数存在且为线性函数，一阶导数连续。
进一步地，根据所述上包络线计算得到最终的呼吸率。
请参阅图5，所述步骤S40包括：
步骤S410：根据预设的第二阈值找出所述上包络线指定区间的最大值点；
步骤S420：根据所述最大值点计算所述上包络线每一点的斜率；
步骤S430：使用所述斜率辅助对所述指定区间的所有呼吸率进行统计分析，计算得到最终的呼吸率。
请一并参阅图6、图7、图8和图9，其中，图6为本发明一实施例提供的PPG信号的示意图，图7为经过预处理后的PPG信号的示意图，即经过了平均滤波器和巴特沃兹带通滤波器的滤波处理；图8为经过第一处理后的第一处理信号，图9为本发明一实施例提供的上包络线的示意图。原始的PPG信号经过处理后得到的上包络线和呼吸信号有较强的相关性，因此运用该上包络线来进行呼吸率的计算。
同时，本发明还提供一种基于PPG信号的呼吸率测量装置，所述基于PPG信号的呼吸率测量装置包括如上所述的基于PPG信号的呼吸率测量方法。
与现有技术相比，本发明提供的一种基于PPG信号的呼吸率测量方法及装置，通过基于PPG的信号实现呼吸率测量，相对于现有的呼吸测量方法受外界的影响小，使用场景广，操作简单易携带，并且测量更快速。
以上所述的仅是本发明的实施方式，在此应当指出，对于本领域的普通技术人员来说，在不脱离本发明创造构思的前提下，还可以做出改进，但这些均属于本发明的保护范围。
Claims (10)
1.一种基于PPG信号的呼吸率测量方法，其特征在于，包括：
步骤S10：实时采集PPG信号；
步骤S20：对所述PPG信号进行第一处理，以得到第一处理信号；所述第一处理包括滤波处理和差分处理；
步骤S30：对所述第一处理信号进行第二处理，以得到上包络线；
步骤S40：根据所述上包络线计算得到最终的呼吸率。
2.根据权利要求1所述的基于PPG信号的呼吸率测量方法，其特征在于，在步骤S10之前还包括：
步骤S1：收集用户的加速度计信号以得到用户运动状态；
步骤S2：判断所述用户运动状态是否为静止状态，如果不是静止状态则执行步骤S3，如果是则执行步骤S4；
步骤S3：提醒用户平静并保持至静止状态，返回步骤S1；
步骤S4：在用户处于静止状态时实时采集PPG信号。
3.根据权利要求1所述的基于PPG信号的呼吸率测量方法，其特征在于，所述PPG信号为16Hz的PPG绿光信号。
4.根据权利要求1所述的基于PPG信号的呼吸率测量方法，其特征在于，所述步骤S20包括：
步骤S210：对所述PPG信号使用一个四阶的平均滤波器和一个四阶的巴特沃兹带通滤波器进行滤波，得到预处理信号；
步骤S220：对所述预处理信号进行一阶差分处理，得到第一处理信号。
5.根据权利要求4所述的基于PPG信号的呼吸率测量方法，其特征在于，所述平均滤波器为：
其中，
X(n)为输入信号，即所述PPG信号；
Y(n)为输出信号；
K和N为自然数。
6.根据权利要求4所述的基于PPG信号的呼吸率测量方法，其特征在于，所述巴特沃兹带通滤波器为：
其中，
aK和bK为选择所需要频率信号计算得到的滤波器的系数。
7.根据权利要求4所述的基于PPG信号的呼吸率测量方法，其特征在于，所述一阶差分为：
Y(n)＝X(n)-X(n-1)。
8.根据权利要求1所述的基于PPG信号的呼吸率测量方法，其特征在于，所述步骤S30包括：
步骤S310：根据预设的第一阈值找出所述第一处理信号中的极大值点；
步骤S320：基于上包络线函数计算三次样条插值得到上包络线；所述上包络线函数在所求区间内[x1,xn]是二次连续可微函数，其中的每个小区间[xk,xk+1](k＝1,2,…,n-1)是一个不超过三次的多项式，其二阶导数存在且为线性函数，一阶导数连续。
9.根据权利要求1所述的基于PPG信号的呼吸率测量方法，其特征在于，所述步骤S40包括：
步骤S410：根据预设的第二阈值找出所述上包络线指定区间的最大值点；
步骤S420：根据所述最大值点计算所述上包络线每一点的斜率；
步骤S430：使用所述斜率辅助对所述指定区间的所有呼吸率进行统计分析，计算得到最终的呼吸率。
10.一种基于PPG信号的呼吸率测量装置，其特征在于，所述基于PPG信号的呼吸率测量装置包括如权利要求1至9任一项所述的基于PPG信号的呼吸率测量方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011180742.9A CN112494008A (zh) | 2020-10-29 | 2020-10-29 | 基于ppg信号的呼吸率测量方法及装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011180742.9A CN112494008A (zh) | 2020-10-29 | 2020-10-29 | 基于ppg信号的呼吸率测量方法及装置 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN112494008A true CN112494008A (zh) | 2021-03-16 |
Family
ID=74954623
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202011180742.9A Pending CN112494008A (zh) | 2020-10-29 | 2020-10-29 | 基于ppg信号的呼吸率测量方法及装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112494008A (zh) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114027822A (zh) * | 2021-04-19 | 2022-02-11 | 北京超思电子技术有限责任公司 | 一种基于ppg信号的呼吸率测量方法及装置 |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106037743A (zh) * | 2016-05-19 | 2016-10-26 | 广东宝莱特医用科技股份有限公司 | 一种提取呼吸率的方法及设备 |
| CN107087862A (zh) * | 2017-01-19 | 2017-08-25 | 上海大学 | 基于msp430f5529穿戴式多功能健康手环系统 |
-
2020
- 2020-10-29 CN CN202011180742.9A patent/CN112494008A/zh active Pending
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106037743A (zh) * | 2016-05-19 | 2016-10-26 | 广东宝莱特医用科技股份有限公司 | 一种提取呼吸率的方法及设备 |
| CN107087862A (zh) * | 2017-01-19 | 2017-08-25 | 上海大学 | 基于msp430f5529穿戴式多功能健康手环系统 |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114027822A (zh) * | 2021-04-19 | 2022-02-11 | 北京超思电子技术有限责任公司 | 一种基于ppg信号的呼吸率测量方法及装置 |
| WO2022222472A1 (zh) * | 2021-04-19 | 2022-10-27 | 北京超思电子技术有限责任公司 | 呼吸率测量方法及装置、电子设备、可读介质 |
| CN114027822B (zh) * | 2021-04-19 | 2022-11-25 | 北京超思电子技术有限责任公司 | 一种基于ppg信号的呼吸率测量方法及装置 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104783780B (zh) | 心电信号除噪方法及装置 | |
| US9986951B1 (en) | Low-distortion ECG denoising | |
| CN103099615A (zh) | 一种消除运动心电信号干扰的方法和装置 | |
| CN106618542A (zh) | 一种去噪心率信号检测装置及方法 | |
| CN106618560B (zh) | 脑电波信号的处理方法和装置 | |
| CN104622440B (zh) | 一种提取脉搏波时标点的方法及装置 | |
| CN117033911B (zh) | 基于智能眼镜数据的计步分析方法 | |
| CN110974189A (zh) | 脉搏波的信号质量检测方法及装置、设备、系统 | |
| CN105662345A (zh) | 心跳信号处理方法、装置和系统 | |
| CN112494001B (zh) | Ppg信号质量评估方法、装置、穿戴式设备及存储介质 | |
| CN110680307A (zh) | 一种运动环境下基于脉搏波传导时间的动态血压监测方法 | |
| CN108042107A (zh) | 一种ppg信号伪差校正方法 | |
| CN112587133A (zh) | 一种血氧饱和度的测量方法 | |
| CN112494008A (zh) | 基于ppg信号的呼吸率测量方法及装置 | |
| CN105816165B (zh) | 一种实时动态心率监测装置及监测方法 | |
| CN106236041B (zh) | 一种实时且准确的测量心率及呼吸率的算法及系统 | |
| TW201811261A (zh) | 訊號偵測方法 | |
| TW201332512A (zh) | 量測心率方法及裝置 | |
| CN107224284A (zh) | 一种全数字心电信号的噪声检测方法及系统 | |
| CN113375065B (zh) | 管道泄漏监测中趋势信号的消除方法及装置 | |
| CN106175755B (zh) | 用于睡眠状态分析的睡眠状态检测器 | |
| CN114463784A (zh) | 一种基于视频-音频多模态深度学习的多人跳绳分析方法 | |
| CN103705242B (zh) | 一种人体呼吸监测设备中吸气波的提取方法 | |
| CN106236067A (zh) | 一种心率信号处理器及信号处理方法 | |
| CN112353377A (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 |