CN110013275A - 一种基于超声多普勒的彩色频谱图及频谱血流检测方法 - Google Patents
一种基于超声多普勒的彩色频谱图及频谱血流检测方法 Download PDFInfo
- Publication number
- CN110013275A CN110013275A CN201910420762.XA CN201910420762A CN110013275A CN 110013275 A CN110013275 A CN 110013275A CN 201910420762 A CN201910420762 A CN 201910420762A CN 110013275 A CN110013275 A CN 110013275A
- Authority
- CN
- China
- Prior art keywords
- blood flow
- frequency spectrum
- color
- detection method
- blood
- 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
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/06—Measuring blood flow
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/48—Diagnostic techniques
- A61B8/488—Diagnostic techniques involving Doppler signals
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/52—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/5215—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data
- A61B8/5223—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data for extracting a diagnostic or physiological parameter from medical diagnostic data
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Medical Informatics (AREA)
- Surgery (AREA)
- Pathology (AREA)
- Radiology & Medical Imaging (AREA)
- Biophysics (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Physics & Mathematics (AREA)
- Molecular Biology (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Hematology (AREA)
- Physiology (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
Abstract
本发明涉及一种基于超声多普勒的彩色频谱图及频谱血流检测方法,实现方法如下:将取样容积在心脏大血管同一部位使用多个门;获得多普勒血流信号随时间变化的各功率谱,选择频带区域;获取血流流速,使用色彩映射表对血流速度进行色彩编码显示,并校正基线到一个频谱图像不发生正负频混叠的位置,获得该部位的血流彩色频谱参考曲线;动态计算并在血流彩色频谱参考曲线中实时显示出血流参数;从频谱曲线上可以了解血流性质、方向、流速,不仅可以更好地观测血流在空间位置的分布情况,而且还能进行更为精确的容积血流速度测量,同时动态计算并在血流彩色频谱参考曲线中实时显示出血流参数,解决了目前对于血流数据测试不全面、不直观的现状问题。
Description
技术领域
本发明涉及彩色频谱图及频谱血流检测技术领域,更具体地说,涉及一种基于超声多普勒的彩色频谱图及频谱血流检测方法。
背景技术
彩色血流成像是把体内血流速度分布进行彩色编码以实时显示,通过提供血流的方向、速度和紊动信息,可以决定血流的流速及其速度分布。传统的血流成像一般采用时域相关法在完成血流速度、方差和能量等表征血流状态的参数,然后经过彩色编码提供给使用者。传统方法的不足在于得到的血流信息是一定区域内的平均信息,对一些关键部位很难进行细致深入的定量分析。
频谱多普勒方法通过向被称为取样门的同一空间位置上聚焦发射重复的散射冲击信号来计算在取样门内的速度分布,并以频谱图形式显示速度分布信息,频谱图是一个随着时间改变的图表,在一个方向上表示时间,另一个方向上表示速度。传统的频谱多普勒是将属于同一个门内的基复信号进行平均从而计算得到时间序列上的频谱或时间维度。
频谱多普勒原理如下所示:
Gd=(2vGsc)cosθ,(v<<c)]]>
其中,v是流速,c为声波的传播速度,Gd表示被侦测出得多普勒频移,Gs为源频率,θ表示血流与超声波束的夹角。
实际的频谱多普勒采用如下公式:
x(d,n)=exp(j2πGcGsd)Σvαvexp(j2π2GcvcGftn)]]>
其中,x()为实际的接收信号,d和n跟观察窗口有关,Gs为快时采样频率,Gft为慢时采样频率,αv为流速分布。
但是,传统频谱多普勒方法也有缺陷:流速的空间信息与速度分布具有同等的重要性,而当取样门中的数据被求和时,血流的空间信息会丢失,血流中紊动信息常常通过Gate(门)中的频谱宽度或流速传播来推断。然而,由于距离血管壁的距离不同,血管中的层流血流的速度剖面类似于一个抛物线形式,传统多普勒显示方法不易于观测。
发明内容
本发明要解决的技术问题在于,针对现有技术的上述缺陷,提供一种基于超声多普勒的彩色频谱图及频谱血流检测方法。
本发明解决其技术问题所采用的技术方案是:
构造一种基于超声多普勒的彩色频谱图及频谱血流检测方法,其实现方法如下:
第一步:将取样容积在心脏大血管同一部位使用多个门;
第二步:获得多普勒血流信号随时间变化的各功率谱,选择频带区域;
第三步:获取血流流速,使用色彩映射表对血流速度进行色彩编码显示,并校正基线到一个频谱图像不发生正负频混叠的位置,获得该部位的血流彩色频谱参考曲线;
第四步:动态计算并在血流彩色频谱参考曲线中实时显示出血流参数。
本发明所述的基于超声多普勒的彩色频谱图及频谱血流检测方法,其中,所述第二步中,功率谱的计算方法包括快速傅里叶变换法。
本发明所述的基于超声多普勒的彩色频谱图及频谱血流检测方法,其中,所述第一步中,通过将取样容积置于心脏大血管的同一部位,沿着血管的横切面在不同的位置使用多个门。
本发明所述的基于超声多普勒的彩色频谱图及频谱血流检测方法,其中,所述第二步中,具体表现在超声多普勒探头经过射频超声回波信号解调和滤波及模数转换后获得连续多普勒血流信号随时间变化的各功率谱,选择频带区域。
本发明所述的基于超声多普勒的彩色频谱图及频谱血流检测方法,其中,所述第三步中,计算流速公式为:
多普勒频移(fd)与发射频率(fo)、血流速度(V)、超声束与血流间夹角(θ)的余弦成正比,与声速(C)成反比,公式为:
fd=±2v/λ=±2v/Cfo
fd=±2v·cosθ/Cfo
V=fdC/2fo·cosθ
式中fd、cosθ仪器均可显示,fo及C为已知,可以计算出V;
由此可得:
ΣH=HsG=HevH×HΣH=HsH=HevH=Μ;
式中,Hs为有效频带的起始点坐标,He为有效频带的终止点坐标,H为有效频带的频点坐标,vH为频点H处的像素值即速度的强度,Μ表示有效频带内各点的像素值。
本发明所述的基于超声多普勒的彩色频谱图及频谱血流检测方法,其中,所述血流参数包括血流速度、血液黏稠度、管径比值、搏动指数和阻力指数中一种或多种。
本发明所述的基于超声多普勒的彩色频谱图及频谱血流检测方法,其中,所述第三步中,所述血流彩色频谱参考曲线中,曲线横轴代表时间,纵轴代表血流速度。
本发明所述的基于超声多普勒的彩色频谱图及频谱血流检测方法,其中,所述第三步中,进行色彩编码,具体反映血液血管内流状态情况。
本发明的有益效果在于:将取样容积在心脏大血管同一部位使用多个门,获得多普勒血流信号随时间变化的各功率谱,选择频带区域;获取血流流速,使用色彩映射表对血流速度进行色彩编码显示,并校正基线到一个频谱图像不发生正负频混叠的位置,获得该部位的血流彩色频谱参考曲线,曲线横轴代表时间,纵轴代表血流速度,从频谱曲线上可以了解血流性质、方向、流速,通过血流彩色频谱曲线变化进而表达血流的改变,不仅可以更好地观测血流在空间位置的分布情况,而且还能进行更为精确的容积血流速度测量;同时动态计算并在血流彩色频谱参考曲线中实时显示出血流参数,血流参数包括血流速度、血液黏稠度、管径比值、搏动指数和阻力指数等,解决了目前对于血流数据测试不全面、不直观的现状问题。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将结合附图及实施例对本发明作进一步说明,下面描述中的附图仅仅是本发明的部分实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他附图:
图1是本发明较佳实施例的基于超声多普勒的彩色频谱图及频谱血流检测方法流程图;
图2是本发明较佳实施例的基于超声多普勒的彩色频谱图及频谱血流检测方法多普勒原理示意图;
图3是本发明较佳实施例的基于超声多普勒的彩色频谱图及频谱血流检测方法的门划分示意图;
图4是本发明较佳实施例的基于超声多普勒的彩色频谱图及频谱血流检测方法频谱计算结果示意图。
具体实施方式
为了使本发明实施例的目的、技术方案和优点更加清楚,下面将结合本发明实施例中的技术方案进行清楚、完整的描述,显然,所描述的实施例是本发明的部分实施例,而不是全部实施例。基于本发明的实施例,本领域普通技术人员在没有付出创造性劳动的前提下所获得的所有其他实施例,都属于本发明的保护范围。
本发明较佳实施例的基于超声多普勒的彩色频谱图及频谱血流检测方法,如图1所示,同时参阅图2-4,其实现方法如下:
S01:将取样容积在心脏大血管同一部位使用多个门;
S02:获得多普勒血流信号随时间变化的各功率谱,选择频带区域;
S03:获取血流流速,使用色彩映射表对血流速度进行色彩编码显示,并校正基线到一个频谱图像不发生正负频混叠的位置,获得该部位的血流彩色频谱参考曲线;
S04:动态计算并在血流彩色频谱参考曲线中实时显示出血流参数;
将取样容积在心脏大血管同一部位使用多个门,获得多普勒血流信号随时间变化的各功率谱,选择频带区域;获取血流流速,使用色彩映射表对血流速度进行色彩编码显示,并校正基线到一个频谱图像不发生正负频混叠的位置,获得该部位的血流彩色频谱参考曲线,曲线横轴代表时间,纵轴代表血流速度,从频谱曲线上可以了解血流性质、方向、流速,通过血流彩色频谱曲线变化进而表达血流的改变,不仅可以更好地观测血流在空间位置的分布情况,而且还能进行更为精确的容积血流速度测量;同时动态计算并在血流彩色频谱参考曲线中实时显示出血流参数,血流参数包括血流速度、血液黏稠度、管径比值、搏动指数和阻力指数等,解决了目前对于血流数据测试不全面、不直观的现状问题。
优选的,第二步中,功率谱的计算方法包括快速傅里叶变换法。
优选的,第一步中,通过将取样容积置于心脏大血管的同一部位,沿着血管的横切面在不同的位置使用多个门。
优选的,第二步中,具体表现在超声多普勒探头经过射频超声回波信号解调和滤波及模数转换后获得连续多普勒血流信号随时间变化的各功率谱,选择频带区域。
优选的,第三步中,计算流速公式为:
多普勒频移(fd)与发射频率(fo)、血流速度(V)、超声束与血流间夹角(θ)的余弦成正比,与声速(C)成反比,公式为:
fd=±2v/λ=±2v/Cfo
fd=±2v·cosθ/Cfo
V=fdC/2fo·cosθ
式中fd、cosθ仪器均可显示,fo及C为已知,可以计算出V;
由此可得:
ΣH=HsG=HevH×HΣH=HsH=HevH=Μ;
式中,Hs为有效频带的起始点坐标,He为有效频带的终止点坐标,H为有效频带的频点坐标,vH为频点H处的像素值即速度的强度,Μ表示有效频带内各点的像素值;根据本公式能够清楚的获得血流参数。
优选的,血流参数包括血流速度、血液黏稠度、管径比值、搏动指数和阻力指数中一种或多种。
优选的,第三步中,血流彩色频谱参考曲线中,曲线横轴代表时间,纵轴代表血流速度;
单位时间内流过血管某一截面的血量称为血流量,也称容积速度,血液在血管流动时,其血流速度与血流量成正比,与血管的截面成反比。根据泊肃叶(Poiseuilli)定律Q=π(P1-P2)r4/8ηL,由Re=VDσ/η式中的V为血液在血管内的平均流速(单位为cm/s),D为管腔直径(单位为cm),σ为血液密度(单位为g/cm3),η为血液沾滞度(单位为泊),Re为Reynolds数,没有单位。一般当Re数超过2000时,就可发生湍流。由上式可知,在血流速度快,血管口径大,血液粘滞度低的情况下,容易产生湍流;
从频谱曲线上可以了解血流性质、方向、流速,通过血流彩色频谱曲线变化进而表达血流的改变,可直观的看出血液在血管里的具体表现,其中层流体现出血液在血管中以相同的方向作规则的分层运动而无横向交换,其横截面上各点的流速出现速度梯度小,轴心快而靠管壁慢、频带窄、频谱光点密集、包络光滑、频谱与基线之间有明显的空窗,分布规律符合泊肃叶方程,多普勒血流声平滑且具有乐感;其中湍流体现出当血流遇到阻碍时,流体成分以多种方向和速度作无规则的运动,流体成分间相互错杂交换,湍流的峰值流速增加、速度分布分散、红细胞运动加速,速度梯度大、频带增宽、频谱光点疏散、包络毛刺、频谱与基线之间无空窗、血流的多普勒声粗糙刺耳。
优选的,第三步中,进行色彩编码,具体反映血液血管内流状态情况,血液血管内流呈规律层流状态处于血管央红细胞流快向周边逐渐减慢所TCD频谱表现红色集周边并蓝色"频窗"规律层流频谱血管现严重狭窄:①狭窄部位血流速度增快处于高流速红细胞数量减少呈现频谱紊乱湍流状态;②由于狭窄血管内径复原或代偿性扩张使处于边缘红细胞形种涡漩反流状态或量处于低流速红细胞血流表现向性典形狭窄血流频谱表现周边蓝色基底部"频窗"消失双向红色涡流或湍流替代。
应当理解的是,对本领域普通技术人员来说,可以根据上述说明加以改进或变换,而所有这些改进和变换都应属于本发明所附权利要求的保护范围。
Claims (8)
1.一种基于超声多普勒的彩色频谱图及频谱血流检测方法,其特征在于,实现方法如下:
第一步:将取样容积在心脏大血管同一部位使用多个门;
第二步:获得多普勒血流信号随时间变化的各功率谱,选择频带区域;
第三步:获取血流流速,使用色彩映射表对血流速度进行色彩编码显示,并校正基线到一个频谱图像不发生正负频混叠的位置,获得该部位的血流彩色频谱参考曲线;
第四步:动态计算并在血流彩色频谱参考曲线中实时显示出血流参数。
2.根据权利要求1所述的基于超声多普勒的彩色频谱图及频谱血流检测方法,其特征在于,所述第二步中,功率谱的计算方法包括快速傅里叶变换法。
3.根据权利要求1所述的基于超声多普勒的彩色频谱图及频谱血流检测方法,其特征在于,所述第一步中,通过将取样容积置于心脏大血管的同一部位,沿着血管的横切面在不同的位置使用多个门。
4.根据权利要求1所述的基于超声多普勒的彩色频谱图及频谱血流检测方法,其特征在于,所述第二步中,具体表现在超声多普勒探头经过射频超声回波信号解调和滤波及模数转换后获得连续多普勒血流信号随时间变化的各功率谱,选择频带区域。
5.根据权利要求1所述的基于超声多普勒的彩色频谱图及频谱血流检测方法,其特征在于,所述第三步中,计算流速公式为:
多普勒频移(fd)与发射频率(fo)、血流速度(V)、超声束与血流间夹角(θ)的余弦成正比,与声速(C)成反比,公式为:
fd=±2v/λ=±2v/Cfo
fd=±2v.cosθ/Cfo
V=fdC/2fo.cosθ
式中fd、cosθ仪器均可显示,fo及C为已知,可以计算出V;
由此可得:
ΣH=HsG=HevH×HΣH=HsH=HevH=Μ;
式中,Hs为有效频带的起始点坐标,He为有效频带的终止点坐标,H为有效频带的频点坐标,vH为频点H处的像素值即速度的强度,Μ表示有效频带内各点的像素值。
6.根据权利要求1所述的基于超声多普勒的彩色频谱图及频谱血流检测方法,其特征在于,所述血流参数包括血流速度、血液黏稠度、管径比值、搏动指数和阻力指数中一种或多种。
7.根据权利要求1所述的基于超声多普勒的彩色频谱图及频谱血流检测方法,其特征在于,所述第三步中,所述血流彩色频谱参考曲线中,曲线横轴代表时间,纵轴代表血流速度。
8.根据权利要求1所述的基于超声多普勒的彩色频谱图及频谱血流检测方法,其特征在于,所述第三步中,进行色彩编码,具体反映血液血管内流状态情况。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910420762.XA CN110013275A (zh) | 2019-05-20 | 2019-05-20 | 一种基于超声多普勒的彩色频谱图及频谱血流检测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910420762.XA CN110013275A (zh) | 2019-05-20 | 2019-05-20 | 一种基于超声多普勒的彩色频谱图及频谱血流检测方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN110013275A true CN110013275A (zh) | 2019-07-16 |
Family
ID=67194090
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910420762.XA Pending CN110013275A (zh) | 2019-05-20 | 2019-05-20 | 一种基于超声多普勒的彩色频谱图及频谱血流检测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110013275A (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111325202A (zh) * | 2020-01-19 | 2020-06-23 | 华中科技大学同济医学院附属协和医院 | 一种基于超声多普勒的彩色频谱图的频谱血流检测方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5564424A (en) * | 1995-06-30 | 1996-10-15 | Siemens Medical Systems, Inc. | Method and apparatus for pulsed doppler ultrasound beam-forming |
CN104095656A (zh) * | 2014-07-25 | 2014-10-15 | 声泰特(成都)科技有限公司 | 一种基于超声多普勒频谱的彩色血流成像及其显示方法 |
CN104146731A (zh) * | 2014-07-21 | 2014-11-19 | 深圳先进技术研究院 | 一种估算流体速度的系统和方法 |
CN105476661A (zh) * | 2016-01-13 | 2016-04-13 | 北京悦琦创通科技有限公司 | 调整多普勒频谱图的方法和设备 |
CN105476665A (zh) * | 2016-01-27 | 2016-04-13 | 成都欣声科技有限公司 | 一种基于超声的血流速度测量及血流流量测量方法 |
CN106388862A (zh) * | 2015-07-27 | 2017-02-15 | 美国西门子医疗解决公司 | 改善的多普勒成像 |
-
2019
- 2019-05-20 CN CN201910420762.XA patent/CN110013275A/zh active Pending
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5564424A (en) * | 1995-06-30 | 1996-10-15 | Siemens Medical Systems, Inc. | Method and apparatus for pulsed doppler ultrasound beam-forming |
CN104146731A (zh) * | 2014-07-21 | 2014-11-19 | 深圳先进技术研究院 | 一种估算流体速度的系统和方法 |
CN104095656A (zh) * | 2014-07-25 | 2014-10-15 | 声泰特(成都)科技有限公司 | 一种基于超声多普勒频谱的彩色血流成像及其显示方法 |
CN106388862A (zh) * | 2015-07-27 | 2017-02-15 | 美国西门子医疗解决公司 | 改善的多普勒成像 |
CN105476661A (zh) * | 2016-01-13 | 2016-04-13 | 北京悦琦创通科技有限公司 | 调整多普勒频谱图的方法和设备 |
CN105476665A (zh) * | 2016-01-27 | 2016-04-13 | 成都欣声科技有限公司 | 一种基于超声的血流速度测量及血流流量测量方法 |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111325202A (zh) * | 2020-01-19 | 2020-06-23 | 华中科技大学同济医学院附属协和医院 | 一种基于超声多普勒的彩色频谱图的频谱血流检测方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Shung et al. | The effects of hematocrit, shear rate, and turbulence on ultrasonic Doppler spectrum from blood | |
Tan et al. | Ultrasonic Doppler technique for application to multiphase flows: A review | |
Sahn | Instrumentation and physical factors related to visualization of stenotic and regurgitant jets by Doppler color flow mapping | |
David et al. | Modern spectral analysis techniques for blood flow velocity and spectral measurements with pulsed Doppler ultrasound | |
CN104095656B (zh) | 一种基于超声多普勒频谱的彩色血流成像及其显示方法 | |
Goddi et al. | High-frame rate vector flow imaging of the carotid bifurcation | |
Guo et al. | Three-dimensional power Doppler imaging: a phantom study to quantify vessel stenosis | |
CN105476665B (zh) | 一种基于超声的血流速度测量及血流流量测量方法 | |
US4932415A (en) | Method of color coding two dimensional ulltrasonic doppler velocity images of blood flow on a display | |
Li et al. | Measurement of mean velocity during pulsatile flow using time-averaged maximum frequency of Doppler ultrasound waveforms | |
Rubin et al. | Volume flow measurement using Doppler and grey-scale decorrelation | |
Boote et al. | Performance tests of Doppler ultrasound equipment with a tissue and blood‐mimicking phantom. | |
Wells | Ultrasonic colour flow imaging | |
Jeronimo et al. | Echo-Lagrangian particle tracking: an ultrasound-based method for extracting path-dependent flow quantities | |
CN110013275A (zh) | 一种基于超声多普勒的彩色频谱图及频谱血流检测方法 | |
Yu et al. | Ultrasonic parametric imaging of erythrocyte aggregation using the structure factor size estimator | |
Walter et al. | Absolute flow measurements using pulsed Doppler US. Work in progress. | |
Tamura et al. | Quantitative study of steady flow using color Doppler ultrasound | |
Browne et al. | Validation of a sensitivity performance index test protocol and evaluation of colour Doppler sensitivity for a range of ultrasound scanners | |
Karabetsos et al. | Design and development of a new ultrasonic doppler technique for estimation of the aggregation of red blood cells | |
Diebold et al. | In vitro flow mapping of regurgitant jets: systematic description of free jet with laser Doppler velocimetry | |
D'Luna et al. | In vitro Doppler detection of axisymmetric stenoses from transverse velocity measurements | |
Jeronimo et al. | Extraction of Particle Residence Time using echo-Lagrangian particle tracking | |
Wendt III et al. | Nuclear magnetic resonance velocity spectra of pulsatile flow in a rigid tube | |
Giarre et al. | Velocity vector reconstruction for color flow Doppler: Experimental evaluation of a new geometrical method |
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 |