CN113925528A - 多普勒成像方法和超声设备 - Google Patents
多普勒成像方法和超声设备 Download PDFInfo
- Publication number
- CN113925528A CN113925528A CN202010609264.2A CN202010609264A CN113925528A CN 113925528 A CN113925528 A CN 113925528A CN 202010609264 A CN202010609264 A CN 202010609264A CN 113925528 A CN113925528 A CN 113925528A
- Authority
- CN
- China
- Prior art keywords
- sampling gate
- sampling
- gate
- echo signal
- fed back
- 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
- 238000003384 imaging method Methods 0.000 title claims abstract description 51
- 238000002604 ultrasonography Methods 0.000 title claims description 25
- 238000005070 sampling Methods 0.000 claims abstract description 295
- 230000005540 biological transmission Effects 0.000 claims abstract description 46
- 238000001514 detection method Methods 0.000 claims abstract description 45
- 238000000034 method Methods 0.000 claims abstract description 45
- 239000000523 sample Substances 0.000 claims description 39
- 238000012545 processing Methods 0.000 claims description 38
- 230000008569 process Effects 0.000 claims description 25
- 238000012163 sequencing technique Methods 0.000 claims description 6
- 230000015572 biosynthetic process Effects 0.000 abstract description 9
- 230000033001 locomotion Effects 0.000 abstract description 8
- 238000003786 synthesis reaction Methods 0.000 abstract description 8
- 230000035945 sensitivity Effects 0.000 abstract description 6
- 230000017531 blood circulation Effects 0.000 description 10
- 238000010586 diagram Methods 0.000 description 10
- 238000001228 spectrum Methods 0.000 description 8
- 238000012937 correction Methods 0.000 description 7
- 238000012986 modification Methods 0.000 description 5
- 230000004048 modification Effects 0.000 description 5
- 238000004590 computer program Methods 0.000 description 4
- 239000002131 composite material Substances 0.000 description 3
- 230000003287 optical effect Effects 0.000 description 3
- 230000003595 spectral effect Effects 0.000 description 3
- 230000000694 effects Effects 0.000 description 2
- 230000006870 function Effects 0.000 description 2
- 239000013307 optical fiber Substances 0.000 description 2
- 230000000644 propagated effect Effects 0.000 description 2
- 238000012935 Averaging Methods 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000000701 chemical imaging Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 238000005192 partition Methods 0.000 description 1
- 230000000737 periodic effect Effects 0.000 description 1
- 230000001902 propagating effect Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 238000012285 ultrasound imaging Methods 0.000 description 1
- 230000003313 weakening effect Effects 0.000 description 1
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/44—Constructional features of the ultrasonic, sonic or infrasonic diagnostic device
-
- 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
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/54—Control of the diagnostic device
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Biomedical Technology (AREA)
- Medical Informatics (AREA)
- Pathology (AREA)
- Radiology & Medical Imaging (AREA)
- Biophysics (AREA)
- Physics & Mathematics (AREA)
- Heart & Thoracic Surgery (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
Abstract
本申请公开了一种多普勒成像方法和超声设备,用于解决现有技术对多采样门成像时各采样门偏转角度一致且成像质量有待提高的问题。该方法逐一对每个取样门进行宽波束聚焦发送;并接收每个取样门反馈的回波信号;当对目标检测区域的所有取样门完成一轮聚焦处理之后,分别对同一取样门反馈的回波信号进行波束合成,得到扫描信息并进行多普勒成像。本申请实施例中,通过对同一目标检测区域的不同取样门采用各自的偏转角度,能够提高对运动检测的灵敏度,进一步通过对同一采样门的回波信号进行波束合成,能够提高信号强度,提高对多采样门同时进行多普勒成像时的成像质量。
Description
技术领域
本申请实施例涉及超声成像技术领域,尤其涉及一种多普勒成像方法和超声设备。
背景技术
在医学超声领域,频谱多普勒技术被广泛的用于血流信息定量检测。该技术能够准确的检测特定位置的血流/组织运动速度,并以频谱的形式将血流/组织运动速度直观的显示出来。频谱图是一个随时间变化的二维图表,在一个方向上表示时间,在另一个方向上表示血流/组织的运动速度。
现有的超声频谱多普勒成像方法大多选定一个取样门,按照设定的脉冲重复频率向同一取样门发射一定长度的脉冲,然后对取样门内的回波信号进行分析,得到取样门内随时间变化的血流速度和方向频谱。
相关技术中也出现了能够同时对多个取样门进行多普勒成像的技术。该技术能够根据用户设定的多个取样门位置从扫描线信息中提取对应的取样门位置的扫描信息。然而不同位置的扫描线能量不同,离聚焦中心越远,扫描线能量越弱,回波信号的信噪比越低,成像质量越差,从而限制了其进一步发展。同时,发明人发现,该技术中多个取样门的设置必须是同一偏转角度,不支持对不同偏转角度的取样门进行频谱成像。
发明内容
本申请实施例提供一种多普勒成像方法和超声设备,用以解决现有技术中成像质量差且取样门不支持任一偏转角度的问题。
根据示例性的实施方式中的一个方面,提供一种多普勒成像方法,其中,目标检测区域中包括预先设置的多个取样门,且各取样门的偏转角度支持任意角度,所述方法包括:
逐一对每个取样门执行聚焦处理,所述聚焦处理包括:以所述取样门为宽波束的聚焦位置,根据所述取样门的发射系数向所述目标检测区域发射宽波束;并接收每个取样门反馈的回波信号;
当对所述目标检测区域的所有取样门完成一轮聚焦处理之后,分别对同一取样门反馈的回波信号进行波束合成,得到扫描信息;
基于各取样门的扫描信息进行多普勒成像。
可选的,各取样门的设置信息中包括位置信息、尺寸信息、深度信息以及偏转角度;
所述方法还包括:
根据各取样门的位置信息、尺寸信息、深度信息以及偏转角度确定用于发送宽波束的发射振元的数量,以使所述发射振元发送的宽波束能够覆盖所述目标检测区域。
可选的,所述分别对同一取样门反馈的回波信号进行波束合成,得到扫描信息,包括:
对任一取样门,对所述取样门反馈的各辅回波信号分别进行相位补偿,所述辅回波信号为对所述取样门之外的取样门执行聚焦处理时所述取样门反馈的回波信号;
对所述取样门的辅回波信号以及主回波信号进行加权求和,得到所述扫描信息;
其中,所述主回波信号为对所述取样门执行聚焦处理时所述取样门反馈的回波信号,且各回波信号的权重与回波信号的幅值成正比。
可选的,所述逐一对每个取样门执行聚焦处理,包括:
根据预设规则确定各取样门的排序;
根据所述排序顺序,依序逐一对每个取样门执行聚焦处理。
可选的,所述发射系数中包括变迹系数、延迟系数以及宽波束类型;
其中,对任一取样门,所述宽波束类型为以下中的任一种:强聚焦波束、弱聚焦波束和平面波。
可选的,所述发射系数中包括宽波束发射频率,所述方法还包括:
根据以下方法确定所述宽波束发射频率:
获取每个取样门的发射周期;其中所述发射周期与发射振元至取样门的距离成正比;
以各取样门中最大发射周期的倒数为所述宽波束发射频率。
根据示例性的实施方式中的另一方面,提供一种超声设备,其中,目标检测区域中包括预先设置的多个取样门,且各取样门的偏转角度支持任意角度,所述超声设备包括:
探头,被配置为发射宽波束,并接收各取样门反馈的回波信号;
显示单元,被配置为显示超声图像;
处理器,分别与所述探头以及所述显示单元相连接,被配置为:
逐一对每个取样门执行聚焦处理,所述聚焦处理包括:以所述取样门为宽波束的聚焦位置,根据所述取样门的发射系数控制所述探头向所述目标检测区域发射宽波束;并通过所述探头接收每个取样门反馈的回波信号;
当对所述目标检测区域的所有取样门完成一轮聚焦处理之后,分别对同一取样门反馈的回波信号进行波束合成,得到扫描信息;
基于各取样门的扫描信息进行多普勒成像。
可选的,各取样门的设置信息中包括位置信息、尺寸信息、深度信息以及偏转角度;
所述处理器还被配置为:
根据各取样门的位置信息、尺寸信息、深度信息以及偏转角度确定所述探头中用于发送宽波束的发射振元的数量,以使所述探头发送的宽波束能够覆盖所述目标检测区域。
可选的,所述处理器执行分别对同一取样门反馈的回波信号进行波束合成,得到扫描信息时,被配置为:
对任一取样门,对所述取样门反馈的各辅回波信号分别进行相位补偿,所述辅回波信号为对所述取样门之外的取样门执行聚焦处理时所述取样门反馈的回波信号;
对所述取样门的辅回波信号以及主回波信号进行加权求和,得到所述扫描信息;
其中,所述主回波信号为对所述取样门执行聚焦处理时所述取样门反馈的回波信号,且各回波信号的权重与回波信号的幅值成正比。
可选的,所述处理器执行逐一对每个取样门执行聚焦处理时,被配置为:
根据预设规则确定各取样门的排序;
根据所述排序顺序,依序逐一对每个取样门执行聚焦处理。
可选的,所述发射系数中包括变迹系数、延迟系数以及宽波束类型;
其中,对任一取样门,所述宽波束类型为以下中的任一种:强聚焦波束、弱聚焦波束和平面波。
可选的,所述发射系数中包括宽波束发射频率,所述处理器还被配置为:
根据以下方法确定所述宽波束发射频率:
获取每个取样门的发射周期;其中所述发射周期与发射振元至取样门的距离成正比;
以各取样门中最大发射周期的倒数为所述宽波束发射频率。
本申请实施例中提供的一种多普勒成像方法及超声设备,本申请实施例中,通过对同一目标检测区域的不同取样门采用各自的偏转角度,能够提高对运动检测的灵敏度,进一步通过对同一采样门的回波信号进行波束合成,能够提高信号强度,提高对多采样门同时进行多普勒成像时的成像质量。
本申请的其它特征和优点将在随后的说明书中阐述,并且,部分地从说明书中变得显而易见,或者通过实施本申请而了解。本申请的目的和其他优点可通过在所写的说明书、权利要求书、以及附图中所特别指出的结构来实现和获得。
附图说明
为了更清楚地说明本申请实施例中的技术方案,下面将对实施例描述中所需要使用的附图作简要介绍,显而易见地,下面描述中的附图仅仅是本申请的一些实施例,对于本领域的普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为根据本申请一个实施例的超声设备的结构示意图;
图2为根据本申请实施例提供的一种多普勒成像方法的应用场景示意图;
图3为根据本申请的一个实施例的应用原理的示意图;
图4为根据本申请实施例提供的一种多普勒成像方法的流程示意图;
图5为根据本申请实施例提供的一种取样门和发送阵列序列的示意图;
图6为根据本申请实施例提供的一种多普勒成像方法的流程示意图。
具体实施方式
为了使本申请的目的、技术方案和优点更加清楚,下面将结合附图对本申请作进一步地详细描述,显然,所描述的实施例仅仅是本申请一部分实施例,而不是全部的实施例。基于本申请中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其它实施例,都属于本申请保护的范围。
图1示出了超声设备100的结构示意图。
下面以超声设备100为例对实施例进行具体说明。应该理解的是,图1所示超声设备100仅是一个范例,并且超声设备100可以具有比图1中所示的更多的或者更少的部件,可以组合两个或多个的部件,或者可以具有不同的部件配置。图中所示出的各种部件可以在包括一个或多个信号处理和/或专用集成电路在内的硬件、软件、或硬件和软件的组合中实现。
图1中示例性示出了根据示例性实施例中超声设备100的硬件配置框图。
如图1所示,超声设备100例如可以包括:处理器110、存储器120、显示单元130和探头140;其中,
探头140,被配置为发射宽波束,并接收各取样门反馈的回波信号;
显示单元130,被配置为显示超声图像;
存储器120被配置为存储用于超声图像所需的数据,可包括软件程序,应用界面数据等;
处理器110,分别与所述探头、所述显示单元以及所述处理器相连接,被配置为:
逐一对每个取样门执行聚焦处理,所述聚焦处理包括:以所述取样门为宽波束的聚焦位置,根据所述取样门的发射系数控制所述探头向所述目标检测区域发射宽波束;并通过所述探头接收每个取样门反馈的回波信号;
当对所述目标检测区域的所有取样门完成一轮聚焦处理之后,分别对同一取样门反馈的回波信号进行波束合成,得到扫描信息;
基于各取样门的扫描信息进行多普勒成像。
在一些可能的实施例中,各取样门的设置信息中包括位置信息、尺寸信息、深度信息以及偏转角度;
所述处理器还被配置为:
根据各取样门的位置信息、尺寸信息、深度信息以及偏转角度确定所述探头中用于发送宽波束的发射振元的数量,以使所述探头发送的宽波束能够覆盖所述目标检测区域。
在一些可能的实施例中,所述处理器执行分别对同一取样门反馈的回波信号进行波束合成,得到扫描信息时,被配置为:
对任一取样门,对所述取样门反馈的各辅回波信号分别进行相位补偿,所述辅回波信号为对所述取样门之外的取样门执行聚焦处理时所述取样门反馈的回波信号;
对所述取样门的辅回波信号以及主回波信号进行加权求和,得到所述扫描信息;
其中,所述主回波信号为对所述取样门执行聚焦处理时所述取样门反馈的回波信号,且各回波信号的权重与回波信号的幅值成正比。
在一些可能的实施例中,所述处理器执行逐一对每个取样门执行聚焦处理时,被配置为:
根据预设规则确定各取样门的排序;
根据所述排序顺序,依序逐一对每个取样门执行聚焦处理。
在一些可能的实施例中,所述发射系数中包括变迹系数、延迟系数以及宽波束类型;
其中,对任一取样门,所述宽波束类型为以下中的任一种:强聚焦波束、弱聚焦波束和平面波。
在一些可能的实施例中,所述发射系数中包括宽波束发射频率,所述处理器还被配置为:
根据以下方法确定所述宽波束发射频率:
获取每个取样门的发射周期;其中所述发射周期与发射振元至取样门的距离成正比;
以各取样门中最大发射周期的倒数为所述宽波束发射频率。
参阅图2,本申请实施例提供的一种多普勒成像方法的应用场景示意图,该场景包括用户200、超声设备201。
当超声波声束方向与血流流向垂直时,根据多普勒效应,此时的频移接近零,这种情况对运动检测的灵敏度较低;当超声声束方向与血流流向一致,根据多普勒效应,此时的频移最大,这种情况对运动检测的灵敏度较高。超声设备在测量血流速度时,通过取样门位置来设定检测深度,通过取样门的偏转角度来控制声束方向,通过取样容积偏转设置声束方向与血流方向的夹角,以便提高频谱多普勒成像的精度和灵敏度。当用户想观察不同位置的血流信息时,需要通过移动取样门的位置、调节取样门的偏转角度、调整取样容积角度来获得最佳的频谱图像。故此,合理设置不同取样门的偏转角度有利于提高成像质量。
本申请提供的方法中,首先,用户在超声设备201的用户界面对目标检测区域的多个取样门进行设置,其中,每个采样门的偏转角度可设置为任意角度,也即各取样门可以按照实际需求设置相应的偏转角度。超声设备201根据用户的设置对目标检测区域发送宽波束以便于获得超声图像。其中,在成像过程中以每个采样为聚焦位置向目标检测区域发射宽波束。同时,无论对那个采样门聚焦,由于宽波束是面向正规目标检测区域发射的,故此,各采样门均会反馈回波信号。在后续的信号处理过程中,对同一采样门反馈的回波信号进行波束合成,来提高同一个采样门的信号质量,提高信噪比,故此能够提升成像质量。由于是基于适用于采样门的偏转角度对每个采样门分别聚焦,所以能够很好的进行多普勒成像。
需要说明的是,上述图2所示的应用场景仅是一种示例,本申请实施例对此不做限定。
图3为根据本申请一个实施例的应用原理的示意图。其中,该部分可由图1所示超声设备的部分模块或功能组件实现,下面将仅针对主要的部件进行说明,而其它部件,如存储器、控制器、控制电路等,此处将不进行赘述。
如图3所示,应用环境中可以包括经由输入输出单元提供的待用户操作的用户界面310、用于显示所述用户界面的显示单元320以及处理器330。
显示单元320可以包括显示面板321、背光组件322。其中,显示面板321被配置为对超声图像进行显示,背光组件322位于显示面板321背面,背光组件322可以包括多个背光分区(图中未示出),各背光分区可以发光,以点亮显示面板321。
处理器330可以被配置为控制背光组件322中各背光分区的背光源亮度,以及控制探头发射宽波束和接收回波信号。
其中处理器330可以包括聚焦处理单元331、波束合成单元332、频谱生成单元333。其中聚焦处理单元331可以被配置为逐一对每个取样门执行聚焦处理,聚焦处理包括:以取样门为宽波束的聚焦位置,根据取样门的发射系数向目标检测区域发射宽波束;并接收每个取样门反馈的回波信号。波束合成单元332被配置为当对目标检测区域的所有取样门完成一轮聚焦处理之后,分别对同一取样门反馈的回波信号进行波束合成,得到扫描信息。频谱生成单元333被配置为基于各取样门的扫描信息进行多普勒成像。
为了能够提高成像质量,图4中详细的示出了本申请实施例提供的一种多普勒成像方法的流程示意图,该方法中目标检测区域中包括预先设置的多个取样门,且各取样门的偏转角度支持任意角度,当用户设置好各个取样门后,可执行以下步骤:
步骤401:逐一对每个取样门执行聚焦处理,聚焦处理包括:以取样门为宽波束的聚焦位置,根据取样门的发射系数向目标检测区域发射宽波束;并接收每个取样门反馈的回波信号。
传统的频谱多普勒成像扫查时通常对单个取样门进行频谱成像,发射的声束比较窄。本申请为了同时观察多个取样门,需要发射的超声声场足够大,并且能够覆盖所有的取样门,同时能量足够强。因此本申请采用宽波束发射方法,使得发射声场能够覆盖用户设置的所有取样门,同时考虑到取样门数量、位置以及偏转角度,自动计算并形成聚焦声场,实现声波能量在不同取样门位置的聚焦。
宽波束发射根据实际取样门的位置,大小以及偏转角度计算并形成聚焦曲线,不再是单一焦点,而是根据取样门位置和偏转角度设置的多个焦点(即每个取样门可视为一个焦点),从而实现不同取样门位置的声波发射能量聚焦。如图5所示:编号分别为0-128的发射振元用于发送宽波束。在发射宽波束时取样门A、B、C能够分别根据各自的发射系数进行聚焦,图中示出了各个取样门A、B、C的偏转角度是不同的。
实现过程可实施为:首先,响应用户操作,确定各个取样门的设置信息,该设置信息中可包括各取样门的位置信息、尺寸信息、深度信息以及偏转角度等;然后根据各取样门的设置信息确定发送宽波束的发射振元的数量,以使发射振元发送的宽波束能够覆盖目标检测区域。其中,每次聚焦一个具体的采样门时发射的宽波束可以不覆盖整个目标检测区域,但是所有发射振元能够对所有采样门进行聚焦并发送宽波束和接收回波信号。也即,如图5所示,在对A采样门聚焦时,可以不覆盖到A采样门之外的采样门,但是整个探头的发射振元和接收振元能够完成对所有采样门的多普勒成像。
在一些实施例中,发射宽波束的发射振元数量由用户设置的取样门位置和数量决定,例如不同的取样门位置较近,那么发射采用的振元数较少;如果不同的取样门位置较远,那么发射采用的振元数较多。所有的发射振元在换能器上是连续排布的,最远的两个取样门分别对应的发射振元之间的距离大于两个横向距离最远的取样门之间的间距。发射宽波束时,可以由部分发射振元发射宽波束,也即,不同取样门对应不同发射振元,对某一取样门聚焦处理时,由该取样门对应的发射振元发射宽波束。
由此,根据用户设置的多个采样门的设置信息确定的发射振元数量,能够确保同时完成对多个采样门的多普勒成像,保证形成一个声场覆盖目标检测区域。
在完成对采样门的设置之后,可以根据各采样门的设置信息生成发射系数,例如变迹系数、延迟系数,以便于能够更好的进行聚焦成像。此外,发射系数中还可以进一步包括宽波束类型,宽波束类型为以下中的任一种:强聚焦波束、弱聚焦波束和平面波。在同一目标检测区域中,不同取样门可采用不同的宽波束类型。例如,一般情况下,宽波束发射采用强聚焦的方式。具体根据用户预设的取样门位置和偏转角度,分别控制发射振元的延迟及变迹,使得发射声场覆盖所有的取样门,发射声场的焦点位于某一个取样门的中心位置,并在各个取样门处接收回波信号。特殊的情况下,当用户预设的取样门数量较多或者距离很远,此时发射的宽波束可采用弱聚焦或平面波,根据用户预设的所有取样门的偏转角度,使发射振元发射弱聚焦声场或者平面波,并按照用户预设的偏转角度依次调整声场偏角,并在各个取样门处接收每个偏转角度下的回波信号。
在另一个实施例中,每个取样门都有对应的发射周期,例如对一个取样门,以该取样门作为聚焦位置发送宽波束时,从发送到接收到回波信号所需的时长可以定义为该取样门的发射周期,故此,发射周期与发射振元至取样门的距离成正比。为了能够保证以一个稳定的频率发送宽波束,且能够保证所有取样门的回波信号都能够被准确的及时接收到,本申请实施例中以各取样门中最大发射周期的倒数为超声设备的宽波束发射频率。例如,宽波束声波在不同取样门位置的发射周期定义为PRT,宽波束发射频率PRF由声波从探头表面传播至最远的取样门的时间决定,例如,若取样门有3个,分别为A、B、C取样门,则宽波束发射频率PRF可定义如公式(1)所示:
其中,max()表示取最大值。
多普勒成像之前可以先对取样门进行排序,得到取样门序列,然后根据排序顺序,依序逐一对每个取样门执行聚焦处理,以便于有序的对各取样门进行多普勒成像。
综上,继续以A、B、C三个取样门为例,发射宽波束和接收回波信号的过程可简述为:
(1)对A取样门进行宽波束聚焦发射,接收A、B、C的回波信号;
(2)对B取样门进行宽波束聚焦发射,接收A、B、C的回波信号;
(3)对C取样门进行宽波束聚焦发射,接收A、B、C的回波信号;
(4)重复上述过程,分别接收A、B、C的回波信号;
可选的特殊情况下的发射-接收序列如下:
(1)根据A取样门的偏转角度发射弱聚焦波或平面波,接收A、B、C的回波信号;
(2)根据B取样门的偏转角度发射弱聚焦波或平面波,接收A、B、C的回波信号;
(3)根据C取样门的偏转角度发射弱聚焦波或平面波,接收A、B、C的回波信号;
(4)重复上述过程,分别接收A、B、C的回波信号。
当然,需要说明的是,以上仅为在三个取样门的情况下的发射-接收序列,当取样门为N时,以上述方式发射和接收,共发射宽波束N次。此外,对不同取样门采用的宽波束类型可不同,例如对A取样门采用强聚焦波束,然后对B取样门采用平面波。具体实施时,可根据实际情况确定不同取样门的宽波束类型。
步骤402:当对目标检测区域的所有取样门完成一轮聚焦处理之后,分别对同一取样门反馈的回波信号进行波束合成,得到扫描信息;
例如,当分别对A、B、C三个取样门完成一轮宽波束聚焦发射后,A取样门共反馈了3次回波信号,即对A取样门聚焦时反馈的回波信号、对B取样门聚焦时反馈的回波信号以及对C取样门聚焦时反馈的回波信号。将A取样门反馈的所有回波信号进行波束符合能够提高信号强度,以此来提供取样门A的成像质量。
在一些实施例中,为了能够很好的提高各取样门的信号强度,本申请实施例中采用相位校正后波束合成的方式,并对波数合成后的IQ信号进行复合加权处理,用于下一步生成频谱。
由于发射的声场是宽波束声场,必然会带来信号能量减弱的问题,本申请为了提升信号能量和信噪比,通过信号复合加权的方式进行信号处理。由于用户预设多个取样门,对于同一个取样门的回波信号在发射和接收环节均存在相位偏移,为了使波束合成后的回波信号能量最大,需要对不同发射周期的信号进行相位校正和复合加权处理。
相位校正的目的是使得同一取样门位置的信号相位一致,以获得加权信号能量最大化,本申请对发射和接收的相位改变均在波束合成环节实现。以图2中取样门A、B、C为例,当聚焦于取样门A时,接收A的回波信号YA→A,A位置的回波信号在t时刻为如公式(2)所示:
当聚焦于取样门B时,接收A的回波信号YB→A,A位置的回波信号在t+Δt时刻表示为如公式(3)所示:
UB→A表示聚焦于取样门B时A取样门的回波信号幅值。
同理,当聚焦于取样门C时,接收A的回波信号YC→A,A位置的回波信号在t+Δt时刻表示为如公式(4)所示:
UC→A表示聚焦于取样门C时A取样门的回波信号幅值。
采用校正值对相位校正,然后对相位校正后的回波信号进行滤波、解调和门内平均处理(如求取信号幅值或信号强度的平均值),得到每个偏转角度下的IQ(in-phasequadrature,同相正交)数据。对每个取样门在不同偏转角度下的IQ数据进行复合加权,得到加权信号,加权公式(5)如下:
Yk=∑iαiYi→k (5)
其中,Yk表示取样门k的IQ信号,它由不同偏转角度下的IQ信号加权而成,α表示加权系数,αi表示第i个取样门的偏转角度对选定取样门k的贡献权重,权重系数由信号幅值大小决定。Yi→k表示第k个取样门在第i个偏转角度下的信号强度。
通过上述过程重复发射宽波束,并接收回波信号进行波束合成处理,可得到不同取样门在时间方向上的扫描信息。
上述波束合成过程可概述为:对任一取样门,将对取样门执行聚焦处理时取样门反馈的回波信号称之为主回波信号,将对取样门之外的取样门执行聚焦处理时取样门反馈的回波信号成为辅回波信号。波束合成时,对取样门反馈的各辅回波信号分别进行相位补偿;然后对取样门的辅回波信号以及主回波信号进行加权求和,得到扫描信息。
步骤403:基于各取样门的扫描信息进行多普勒成像。
综上,本申请实施例中,通过对同一目标检测区域的不同取样门采用各自的偏转角度,能够提高对运动检测的灵敏度,进一步通过对同一采样门的回波信号进行波束合成,能够提高信号强度,提高对多采样门同时进行多普勒成像时的成像质量。
在另一个实施例中,用户可以基于需求,选择其中的一个采样门的超声图像进行查看,也可以选择多个采样门的超声图像同时进行查看和对比。在对比时,为了能够准确的比对出不同采样门之间的差距,可以对不同采样门的扫描信息的信号强度进行归一化,然后基于归一化的结果生成超声图像。
如图6所示,为本申请实施例提供的多普勒成像方法的整体流程图,可包括以下步骤:
步骤601:响应于用户操作,获得目标检测区域内各个取样门的设置信息。
该设置信息例如取样门的位置信息、尺寸信息、偏转角度、深度信息等。
步骤602:根据各个取样门的设置信息,确定发射振元的数量,并确定各个取样门的发送系数。
步骤603:基于预设的取样门序列,逐一对每个取样门执行聚焦处理,聚焦处理包括:以取样门为宽波束的聚焦位置,根据取样门的发射系数向目标检测区域发射宽波束;并接收每个取样门反馈的回波信号。
步骤604:当对目标检测区域的所有取样门完成一轮聚焦处理之后,分别对同一取样门反馈的回波信号进行波束合成,得到扫描信息,并返回执行步骤603。
步骤605:基于各取样门的扫描信息进行多普勒成像。
步骤606:接收到结束多普勒成像的指示时,结束对各采样门的多普勒成像。
在一些可能的实施方式中,本申请实施例提供的方法的各个方面还可以实现为一种程序产品的形式,其包括程序代码,当所述程序代码在计算机设备上运行时,所述程序代码用于使所述计算机设备执行本说明书中描述的根据本申请各种示例性实施方式的数据处理的方法中的步骤。
所述程序产品可以采用一个或多个可读介质的任意组合。可读介质可以是可读信号介质或者可读存储介质。可读存储介质例如可以是——但不限于——电、磁、光、电磁、红外线、或半导体的系统、装置或器件,或者任意以上的组合。可读存储介质的更具体的例子(非穷举的列表)包括:具有一个或多个导线的电连接、便携式盘、硬盘、随机存取存储器(RAM)、只读存储器(ROM)、可擦式可编程只读存储器(EPROM或闪存)、光纤、便携式紧凑盘只读存储器(CD-ROM)、光存储器件、磁存储器件、或者上述的任意合适的组合。
根据本申请的实施方式的用于执行数据处理的程序产品,其可以采用便携式紧凑盘只读存储器(CD-ROM)并包括程序代码,并可以在服务器设备上运行。然而,本申请的程序产品不限于此,在本文件中,可读存储介质可以是任何包含或存储程序的有形介质,该程序可以被信息传输、装置或者器件使用或者与其结合使用。
可读信号介质可以包括在基带中或者作为载波一部分传播的数据信号,其中承载了可读程序代码。这种传播的数据信号可以采用多种形式,包括——但不限于——电磁信号、光信号或上述的任意合适的组合。可读信号介质还可以是可读存储介质以外的任何可读介质,该可读介质可以发送、传播或者传输用于由周期网络动作系统、装置或者器件使用或者与其结合使用的程序。
可读介质上包含的程序代码可以用任何适当的介质传输,包括——但不限于——无线、有线、光缆、RF等,或者上述的任意合适的组合。
可以以一种或多种程序设计语言的任意组合来编写用于执行本申请操作的程序代码,所述程序设计语言包括面向对象的程序设计语言一诸如Java、C++等,还包括常规的过程式程序设计语言一诸如″C″语言或类似的程序设计语言。程序代码可以完全地在用户计算设备上执行、部分地在用户设备上执行、作为一个独立的软件包执行、部分在用户计算设备上部分在远程计算设备上执行、或者完全在远程计算设备或服务器上执行。在涉及远程计算设备的情形中,远程计算设备可以通过任意种类的网络——包括局域网(LAN)或广域网(WAN)一连接到用户计算设备,或者,可以连接到外部计算设备。
本申请实施例执行数据处理的方法还提供一种计算设备可读存储介质,即断电后内容不丢失。该存储介质中存储软件程序,包括程序代码,当所述程序代码在计算设备上运行时,该软件程序在被一个或多个处理器读取并执行时可实现本申请实施例上面任何一种数据处理的方案。
以上参照示出根据本申请实施例的方法、装置(系统)和/或计算机程序产品的框图和/或流程图描述本申请。应理解,可以通过计算机程序指令来实现框图和/或流程图示图的一个块以及框图和/或流程图示图的块的组合。可以将这些计算机程序指令提供给通用计算机、专用计算机的处理器和/或其它可编程数据处理装置,以产生机器,使得经由计算机处理器和/或其它可编程数据处理装置执行的指令创建用于实现框图和/或流程图块中所指定的功能/动作的方法。
相应地,还可以用硬件和/或软件(包括固件、驻留软件、微码等)来实施本申请。更进一步地,本申请可以采取计算机可使用或计算机可读存储介质上的计算机程序产品的形式,其具有在介质中实现的计算机可使用或计算机可读程序代码,以由指令执行系统来使用或结合指令执行系统而使用。在本申请上下文中,计算机可使用或计算机可读介质可以是任意介质,其可以包含、存储、通信、传输、或传送程序,以由指令执行系统、装置或设备使用,或结合指令执行系统、装置或设备使用。
尽管结合具体特征及其实施例对本申请进行了描述,显而易见的,在不脱离本申请的精神和范围的情况下,可对其进行各种修改和组合。相应地,本说明书和附图仅仅是所附权利要求所界定的本申请的示例性说明,且视为已覆盖本申请范围内的任意和所有修改、变化、组合或等同物。
显然,本领域的技术人员可以对本申请进行各种改动和变型而不脱离本申请的范围。这样,倘若本申请的这些修改和变型属于本申请权利要求及其等同技术的范围之内,则本申请也意图包括这些改动和变型在内。
Claims (10)
1.一种多普勒成像方法,其特征在于,目标检测区域中包括预先设置的多个取样门,且各取样门的偏转角度支持任意角度,所述方法包括:
逐一对每个取样门执行聚焦处理,所述聚焦处理包括:以所述取样门为宽波束的聚焦位置,根据所述取样门的发射系数向所述目标检测区域发射宽波束;并接收每个取样门反馈的回波信号;
当对所述目标检测区域的所有取样门完成一轮聚焦处理之后,分别对同一取样门反馈的回波信号进行波束合成,得到扫描信息;
基于各取样门的扫描信息进行多普勒成像。
2.根据权利要求1所述的方法,其特征在于,各取样门的设置信息中包括位置信息、尺寸信息、深度信息以及偏转角度;
所述方法还包括:
根据各取样门的位置信息、尺寸信息、深度信息以及偏转角度确定用于发送宽波束的发射振元的数量,以使所述发射振元发送的宽波束能够覆盖所述目标检测区域。
3.根据权利要求1所述的方法,其特征在于,所述分别对同一取样门反馈的回波信号进行波束合成,得到扫描信息,包括:
对任一取样门,对所述取样门反馈的各辅回波信号分别进行相位补偿,所述辅回波信号为对所述取样门之外的取样门执行聚焦处理时所述取样门反馈的回波信号;
对所述取样门的辅回波信号以及主回波信号进行加权求和,得到所述扫描信息;
其中,所述主回波信号为对所述取样门执行聚焦处理时所述取样门反馈的回波信号,且各回波信号的权重与回波信号的幅值成正比。
4.根据权利要求1所述的方法,其特征在于,所述逐一对每个取样门执行聚焦处理,包括:
根据预设规则确定各取样门的排序;
根据所述排序顺序,依序逐一对每个取样门执行聚焦处理。
5.根据权利要求1所述的方法,其特征在于,所述发射系数中包括变迹系数、延迟系数以及宽波束类型;
其中,对任一取样门,所述宽波束类型为以下中的任一种:强聚焦波束、弱聚焦波束和平面波。
6.根据权利要求1-5中任一所述的方法,其特征在于,所述发射系数中包括宽波束发射频率,所述方法还包括:
根据以下方法确定所述宽波束发射频率:
获取每个取样门的发射周期;其中所述发射周期与发射振元至取样门的距离成正比;
以各取样门中最大发射周期的倒数为所述宽波束发射频率。
7.一种超声设备,其特征在于,目标检测区域中包括预先设置的多个取样门,且各取样门的偏转角度支持任意角度,所述超声设备包括:
探头,被配置为发射宽波束,并接收各取样门反馈的回波信号;
显示单元,被配置为显示超声图像;
处理器,分别与所述探头以及所述显示单元相连接,被配置为:
逐一对每个取样门执行聚焦处理,所述聚焦处理包括:以所述取样门为宽波束的聚焦位置,根据所述取样门的发射系数控制所述探头向所述目标检测区域发射宽波束;并通过所述探头接收每个取样门反馈的回波信号;
当对所述目标检测区域的所有取样门完成一轮聚焦处理之后,分别对同一取样门反馈的回波信号进行波束合成,得到扫描信息;
基于各取样门的扫描信息进行多普勒成像。
8.根据权利要求7所述的超声设备,其特征在于,各取样门的设置信息中包括位置信息、尺寸信息、深度信息以及偏转角度;
所述处理器还被配置为:
根据各取样门的位置信息、尺寸信息、深度信息以及偏转角度确定所述探头中用于发送宽波束的发射振元的数量,以使所述探头发送的宽波束能够覆盖所述目标检测区域。
9.根据权利要求7所述的超声设备,其特征在于,所述处理器执行分别对同一取样门反馈的回波信号进行波束合成,得到扫描信息时,被配置为:
对任一取样门,对所述取样门反馈的各辅回波信号分别进行相位补偿,所述辅回波信号为对所述取样门之外的取样门执行聚焦处理时所述取样门反馈的回波信号;
对所述取样门的辅回波信号以及主回波信号进行加权求和,得到所述扫描信息;
其中,所述主回波信号为对所述取样门执行聚焦处理时所述取样门反馈的回波信号,且各回波信号的权重与回波信号的幅值成正比。
10.根据权利要求7所述的超声设备,其特征在于,所述处理器执行逐一对每个取样门执行聚焦处理时,被配置为:
根据预设规则确定各取样门的排序;
根据所述排序顺序,依序逐一对每个取样门执行聚焦处理。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010609264.2A CN113925528B (zh) | 2020-06-29 | 2020-06-29 | 多普勒成像方法和超声设备 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010609264.2A CN113925528B (zh) | 2020-06-29 | 2020-06-29 | 多普勒成像方法和超声设备 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113925528A true CN113925528A (zh) | 2022-01-14 |
CN113925528B CN113925528B (zh) | 2023-11-03 |
Family
ID=79272723
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010609264.2A Active CN113925528B (zh) | 2020-06-29 | 2020-06-29 | 多普勒成像方法和超声设备 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113925528B (zh) |
Citations (19)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH1099333A (ja) * | 1996-10-01 | 1998-04-21 | Toshiba Iyou Syst Eng Kk | 超音波カラードプラ診断装置および超音波カラードプライメージングの信号処理方法 |
JPH10165400A (ja) * | 1996-12-12 | 1998-06-23 | Furuno Electric Co Ltd | 超音波診断装置 |
JP2002336245A (ja) * | 2001-04-27 | 2002-11-26 | Ge Medical Systems Global Technology Co Llc | 超音波ビーム調整方法および装置並びに超音波診断装置 |
US20080306386A1 (en) * | 2007-06-05 | 2008-12-11 | Tatsuro Baba | Ultrasonic diagnostic apparatus and method of measuring velocity with ultrasonic waves |
JP2009247511A (ja) * | 2008-04-04 | 2009-10-29 | Hitachi Ltd | 超音波撮像装置 |
KR20110069190A (ko) * | 2009-12-17 | 2011-06-23 | 삼성메디슨 주식회사 | 신호 보정 처리를 수행하는 초음파 시스템 및 방법 |
CN102319087A (zh) * | 2011-07-13 | 2012-01-18 | 无锡祥生医学影像有限责任公司 | 触摸屏超声诊断仪及其脉冲波多普勒模式调节方法 |
CN102871645A (zh) * | 2011-07-11 | 2013-01-16 | 浙江大学 | 近红外成像超声血管治疗仪 |
US20130144166A1 (en) * | 2011-12-01 | 2013-06-06 | Donald F. Specht | Motion detection using ping-based and multiple aperture doppler ultrasound |
CN103356242A (zh) * | 2012-03-29 | 2013-10-23 | 深圳市蓝韵实业有限公司 | 脉冲重复频率调节方法及高脉冲多普勒成像装置 |
JP2014068980A (ja) * | 2012-10-01 | 2014-04-21 | Hitachi Aloka Medical Ltd | 超音波診断装置 |
CN103913513A (zh) * | 2014-03-26 | 2014-07-09 | 深圳大学 | 相控阵全局聚焦系统及其聚焦方法 |
CN104739449A (zh) * | 2015-04-15 | 2015-07-01 | 惠尔图像公司 | 超声波束合成方法 |
CN105167802A (zh) * | 2015-09-10 | 2015-12-23 | 深圳华声医疗技术有限公司 | 多普勒成像方法及装置 |
WO2016176855A1 (zh) * | 2015-05-07 | 2016-11-10 | 深圳迈瑞生物医疗电子股份有限公司 | 血流成像方法及系统 |
WO2017047328A1 (ja) * | 2015-09-17 | 2017-03-23 | 株式会社日立製作所 | 超音波診断装置、及び超音波撮像方法 |
CN108078590A (zh) * | 2018-01-03 | 2018-05-29 | 声泰特(成都)科技有限公司 | 基于超声频谱多普勒的血流动力学可视化方法与系统 |
CN109452954A (zh) * | 2017-09-06 | 2019-03-12 | 深圳迈瑞生物医疗电子股份有限公司 | 超声成像方法及装置 |
CN110613476A (zh) * | 2018-06-19 | 2019-12-27 | 青岛海信医疗设备股份有限公司 | 超声信号的处理方法和装置 |
-
2020
- 2020-06-29 CN CN202010609264.2A patent/CN113925528B/zh active Active
Patent Citations (19)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH1099333A (ja) * | 1996-10-01 | 1998-04-21 | Toshiba Iyou Syst Eng Kk | 超音波カラードプラ診断装置および超音波カラードプライメージングの信号処理方法 |
JPH10165400A (ja) * | 1996-12-12 | 1998-06-23 | Furuno Electric Co Ltd | 超音波診断装置 |
JP2002336245A (ja) * | 2001-04-27 | 2002-11-26 | Ge Medical Systems Global Technology Co Llc | 超音波ビーム調整方法および装置並びに超音波診断装置 |
US20080306386A1 (en) * | 2007-06-05 | 2008-12-11 | Tatsuro Baba | Ultrasonic diagnostic apparatus and method of measuring velocity with ultrasonic waves |
JP2009247511A (ja) * | 2008-04-04 | 2009-10-29 | Hitachi Ltd | 超音波撮像装置 |
KR20110069190A (ko) * | 2009-12-17 | 2011-06-23 | 삼성메디슨 주식회사 | 신호 보정 처리를 수행하는 초음파 시스템 및 방법 |
CN102871645A (zh) * | 2011-07-11 | 2013-01-16 | 浙江大学 | 近红外成像超声血管治疗仪 |
CN102319087A (zh) * | 2011-07-13 | 2012-01-18 | 无锡祥生医学影像有限责任公司 | 触摸屏超声诊断仪及其脉冲波多普勒模式调节方法 |
US20130144166A1 (en) * | 2011-12-01 | 2013-06-06 | Donald F. Specht | Motion detection using ping-based and multiple aperture doppler ultrasound |
CN103356242A (zh) * | 2012-03-29 | 2013-10-23 | 深圳市蓝韵实业有限公司 | 脉冲重复频率调节方法及高脉冲多普勒成像装置 |
JP2014068980A (ja) * | 2012-10-01 | 2014-04-21 | Hitachi Aloka Medical Ltd | 超音波診断装置 |
CN103913513A (zh) * | 2014-03-26 | 2014-07-09 | 深圳大学 | 相控阵全局聚焦系统及其聚焦方法 |
CN104739449A (zh) * | 2015-04-15 | 2015-07-01 | 惠尔图像公司 | 超声波束合成方法 |
WO2016176855A1 (zh) * | 2015-05-07 | 2016-11-10 | 深圳迈瑞生物医疗电子股份有限公司 | 血流成像方法及系统 |
CN105167802A (zh) * | 2015-09-10 | 2015-12-23 | 深圳华声医疗技术有限公司 | 多普勒成像方法及装置 |
WO2017047328A1 (ja) * | 2015-09-17 | 2017-03-23 | 株式会社日立製作所 | 超音波診断装置、及び超音波撮像方法 |
CN109452954A (zh) * | 2017-09-06 | 2019-03-12 | 深圳迈瑞生物医疗电子股份有限公司 | 超声成像方法及装置 |
CN108078590A (zh) * | 2018-01-03 | 2018-05-29 | 声泰特(成都)科技有限公司 | 基于超声频谱多普勒的血流动力学可视化方法与系统 |
CN110613476A (zh) * | 2018-06-19 | 2019-12-27 | 青岛海信医疗设备股份有限公司 | 超声信号的处理方法和装置 |
Non-Patent Citations (1)
Title |
---|
高晗 等: "基于线阵超声探头多波束并行发射的颈动脉血流多点频谱多普勒分析", 《中国医疗设备》, pages 1 - 5 * |
Also Published As
Publication number | Publication date |
---|---|
CN113925528B (zh) | 2023-11-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US11624816B2 (en) | Method and system for performing retrospective dynamic transmit focussing beamforming on ultrasound signals | |
US11717270B2 (en) | Method, apparatus and system for imaging in ultrasonic scanning | |
US10359515B2 (en) | Angle independent velocity spectrum determination | |
US11980503B2 (en) | Ultrasound signal processing device and ultrasound diagnostic device | |
KR102332464B1 (ko) | 분석 데이터로부터 전단파 추정 | |
US7666142B2 (en) | Ultrasound doppler diagnostic apparatus and image data generating method | |
US11160536B2 (en) | Ultrasound method and ultrasound system for real time automatic setting of parameters for doppler imaging modes | |
US20200222029A1 (en) | Information processing apparatus, information processing method, and storage medium | |
EP2766738B1 (en) | Three dimensional (3d) transverse oscillation vector velocity ultrasound imaging | |
KR20120121229A (ko) | 빔포머, 진단시스템, 의료영상시스템 및 진단영상을 표시하는 방법 | |
JP6819184B2 (ja) | 超音波測定装置および超音波測定装置の制御方法 | |
JP6746895B2 (ja) | 超音波診断装置、及び超音波信号処理方法 | |
US20200080969A1 (en) | Method for acquiring signals by ultrasound probing, corresponding computer program and ultrasound probing device | |
CN113925528B (zh) | 多普勒成像方法和超声设备 | |
EP3050514A1 (en) | Ultrasonic diagnostic apparatus and method for controlling the same | |
KR101356616B1 (ko) | 가중치 보간 방법 및 이를 이용한 초음파 진단장치 | |
JP6492230B2 (ja) | スペクトル解析装置、スペクトル解析方法及び超音波撮像装置 | |
CN111855824B (zh) | 超声装置及其控制方法 | |
US20200241134A1 (en) | Three-Dimensional (3D) and/or Four-Dimensional (4D) Ultrasound Imaging | |
WO2018195873A1 (zh) | 一种超声连续波多普勒成像方法及装置、存储介质 | |
US20240142620A1 (en) | Ultrasound imaging method and ultrasound imaging system | |
JP7066487B2 (ja) | 超音波診断装置、医用画像処理装置及び医用画像処理プログラム | |
US20150335314A1 (en) | Ultrasonic diagnostic device | |
US20210255321A1 (en) | Ultrasound imaging system with pixel extrapolation image enhancement | |
JP2004174236A (ja) | 超音波撮像方法及び超音波撮像装置 |
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 |