CN112641468B - 超声成像方法、超声图像处理方法、超声设备及存储介质 - Google Patents

超声成像方法、超声图像处理方法、超声设备及存储介质 Download PDF

Info

Publication number
CN112641468B
CN112641468B CN202011479604.0A CN202011479604A CN112641468B CN 112641468 B CN112641468 B CN 112641468B CN 202011479604 A CN202011479604 A CN 202011479604A CN 112641468 B CN112641468 B CN 112641468B
Authority
CN
China
Prior art keywords
sub
blood flow
receiving
dimensional model
receiving lines
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.)
Active
Application number
CN202011479604.0A
Other languages
English (en)
Other versions
CN112641468A (zh
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.)
Qingdao Hisense Medical Equipment Co Ltd
Original Assignee
Qingdao Hisense Medical Equipment Co Ltd
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 Qingdao Hisense Medical Equipment Co Ltd filed Critical Qingdao Hisense Medical Equipment Co Ltd
Priority to CN202011479604.0A priority Critical patent/CN112641468B/zh
Publication of CN112641468A publication Critical patent/CN112641468A/zh
Application granted granted Critical
Publication of CN112641468B publication Critical patent/CN112641468B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/06Measuring blood flow
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/48Diagnostic techniques
    • A61B8/488Diagnostic techniques involving Doppler signals
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/52Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/5215Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Molecular Biology (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Pathology (AREA)
  • Radiology & Medical Imaging (AREA)
  • Physics & Mathematics (AREA)
  • Biomedical Technology (AREA)
  • Veterinary Medicine (AREA)
  • Medical Informatics (AREA)
  • Biophysics (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • General Health & Medical Sciences (AREA)
  • Public Health (AREA)
  • Hematology (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Ultra Sonic Daignosis Equipment (AREA)

Abstract

本申请实施例提供了一种超声成像方法、超声图像处理方法、超声设备及存储介质。通过将预先选定的血流成像区域划分成多个子区域。按照划分顺序,向各子区域内分别发射超声波,以此得到每个子区域对应的接收线。在确定每个子区域对应的接收线后,按照每条接收线对应的子取样门位置对接收线进行划分,以此得到各子取样门对应的子接收线。通过上述步骤所得到的,各子区域的接收线和各子取样门对应的子接收线构建血流成像区域的四维模型。该四维模型能够描述血流成像区域内的血流在时间、血流速度、成像深度以及接收线这四个维度的变化情况。以此解决相关技术中存在只能对特定的血流区域进行评估和诊断,不能对血流区域进行全面且准确的评估的问题。

Description

超声成像方法、超声图像处理方法、超声设备及存储介质
技术领域
本申请涉及超声成像技术领域,特别涉及一种超声成像方法、超声图像处理方法、超声设备及存储介质。
背景技术
在医学超声领域,频谱多普勒技术被广泛的用于血流信息定量检测。相关技术中,针对指定区域内的血流信息进行检测时,预先在该指定区域内设有取样门,然后根据预先设定的时间和频率,以脉冲波的形式向取样门内发射聚焦超声波。在接收到取样门内反馈的超声回波数据之后,对该超声回波数据进行均值处理,以此得到随时间变化的血流信号。通过对血流信号进行频谱分析,来进一步得到血流区域内的血流随着时间变化的频谱图像。该频谱图像的横向表示时间,纵向表示血流速度。由于将取样门内的数据进行了平均,导致血流的空间信息丢失,因此相关技术中,也有提出将取样门分为多个子取样门,通过对各子取样门求均值的方式从血流信号中提取出血流的空间信息。
而上述的成像方法,都存在着只能对特定的血流区域进行评估和诊断,不能对血流区域进行全面且准确的评估的问题。
发明内容
本申请的目的是提供一种超声成像方法、超声图像处理方法、超声设备及存储介质,用于解决相关技术中存在着只能对特定的血流区域进行评估和诊断,不能对血流区域进行全面且准确的评估的问题。
第一方面,本申请实施例提供了一种超声成像方法,所述方法包括:
将预先选定的血流成像区域划分为多个子区域;
依序对每个所述子区域分别发射超声波得到每个所述子区域各自的接收线;
对每条接收线分别执行,按照所述接收线对应的子取样门位置,将所述接收线划分为多份,得到各子取样门分别对应的子接收线;
依据各子区域的接收线,以及各子取样门对应的子接收线,构建所述血流成像区域的四维模型,所述四维模型用于描述:时间、血流速度、成像深度以及接收线之间的关联关系。
在一些可能的实施例中,所述将预先选定的血流成像区域划分为多个子区域,包括:
基于单次发射超声波后得到的接收线条数和所述血流成像区域的总接收线条数,对所述血流成像区域进行划分操作,得到所述多个子区域;其中,所述多个子区域的个数与所述总接收线条数具有反比关系,且与单次发射超声波后得到的接收线条数具有反比关系。
在一些可能的实施例中,所述方法还包括:
针对不同子区域采用分时检测扫描的方式,在一个扫描周期内依序向各子区域分别发送超声波、重复n个所述扫描周期,获得各子区域在不同扫描周期内的各自的接收线;其中,n为正整数;
针对同一子区域,采用插值处理方式,得到相邻两个扫描周期之间的接收线。
在一些可能的实施例中,每条接收线的多个子取样门中位置相邻的子取样门之间具有重叠区域。
在一些可能的实施例中,针对每条接收线,所述接收线的子取样门的数量与所述接收线的取样门内的总点数成正比、与单个子取样门内设置的点数成反比,且与子取样门之间的重叠区域的大小成反比。
第二方面,本申请实施例提供了一种超声图像处理方法,所述方法包括:
接收针对血流成像区域的四维模型的查看指示,其中,所述四维模型用于描述:时间、血流速度、成像深度以及接收线之间的关联关系;
响应所述查看指示,展示所述四维模型中与所述查看指示对象的成像信息。
在一些可能的实施例中,所述查看指示用于从所述四维模型中获取与所述查看指示对应的数据构建所述四维模型中的时间、血流速度、成像深度以及接收线表示四个维度中的任一维度;
所述响应所述查看指示,展示所述四维模型中与所述查看指示对象的成像信息,包括:
从所述四维模型中,获取选定的维度对应的其他三个维度的信息,构建三维模型并展示。
在一些可能的实施例中,所述方法还包括:
构建所述四维模型包括:
将所述血流成像区域划分为多个子区域;
依序对每个所述子区域分别发射超声波得到每个所述子区域各自的接收线;
对每条接收线分别执行,按照所述接收线对应的子取样门位置,将所述接收线划分为多份,得到各子取样门分别对应的子接收线;
依据各子区域的接收线,以及各子取样门对应的子接收线,构建所述血流成像区域的四维模型。
在一些可能的实施例中,所述将预先选定的血流成像区域划分为多个子区域,包括:
基于单次发射超声波后得到的接收线条数和所述血流成像区域的总接收线条数,对所述血流成像区域进行划分操作,得到所述多个子区域;其中,所述多个子区域的个数与所述总接收线条数具有反比关系,且与单次发射超声波后得到的接收线条数具有反比关系。
第三方面,本申请实施例提供了一种超声设备,血流成像区域中包括预先设置的多个取样门,且各取样门的偏转角度支持任意角度,所述超声设备包括:
探头,被配置为发射宽波束,并接收各取样门反馈的回波信号;
显示单元,被配置为显示超声图像;
处理器,分别与所述探头以及所述显示单元相连接,被配置为:
将预先选定的血流成像区域划分为多个子区域;
依序对每个所述子区域分别发射超声波得到每个所述子区域各自的接收线;
对每条接收线分别执行,按照所述接收线对应的子取样门位置,将所述接收线划分为多份,得到各子取样门分别对应的子接收线;
依据各子区域的接收线,以及各子取样门对应的子接收线,构建所述血流成像区域的四维模型,所述四维模型用于描述:时间、血流速度、成像深度以及接收线之间的关联关系。
在一些可能的实施例中,所述处理器执行将预先选定的血流成像区域划分为多个子区域时,被配置为:
基于单次发射超声波后得到的接收线条数和所述血流成像区域的总接收线条数,对所述血流成像区域进行划分操作,得到所述多个子区域;其中,所述多个子区域的个数与所述总接收线条数具有反比关系,且与单次发射超声波后得到的接收线条数具有反比关系。
在一些可能的实施例中,所述处理器还被配置为:
针对不同子区域采用分时检测扫描的方式,在一个扫描周期内依序向各子区域分别发送超声波、重复n个所述扫描周期,获得各子区域在不同扫描周期内的各自的接收线;其中,n为正整数;
针对同一子区域,采用插值处理方式,得到相邻两个扫描周期之间的接收线。
在一些可能的实施例中,每条接收线的多个子取样门中位置相邻的子取样门之间具有重叠区域。
在一些可能的实施例中,针对每条接收线,所述接收线的子取样门的数量与所述接收线的取样门内的总点数成正比、与单个子取样门内设置的点数成反比,且与子取样门之间的重叠区域的大小成反比。
在一些可能的实施例中,所述处理器还被配置为:
接收针对血流成像区域的四维模型的查看指示,其中,所述四维模型用于描述:时间、血流速度、成像深度以及接收线之间的关联关系;
响应所述查看指示,展示所述四维模型中与所述查看指示对象的成像信息。
在一些可能的实施例中,所述查看指示用于从所述四维模型中获取与所述查看指示对应的数据构建所述四维模型中的时间、血流速度、成像深度以及接收线表示四个维度中的任一维度;
所述处理器响应所述查看指示,执行展示所述四维模型中与所述查看指示对象的成像信息时,被配置为:
从所述四维模型中,获取选定的维度对应的其他三个维度的信息,构建三维模型并展示。
在一些可能的实施例中,所述处理器还被配置为:
构建所述四维模型包括:
将所述血流成像区域划分为多个子区域;
依序对每个所述子区域分别发射超声波得到每个所述子区域各自的接收线;
对每条接收线分别执行,按照所述接收线对应的子取样门位置,将所述接收线划分为多份,得到各子取样门分别对应的子接收线;
依据各子区域的接收线,以及各子取样门对应的子接收线,构建所述血流成像区域的四维模型。
在一些可能的实施例中,所述处理器执行将预先选定的血流成像区域划分为多个子区域时,被配置为:
基于单次发射超声波后得到的接收线条数和所述血流成像区域的总接收线条数,对所述血流成像区域进行划分操作,得到所述多个子区域;其中,所述多个子区域的个数与所述总接收线条数具有反比关系,且与单次发射超声波后得到的接收线条数具有反比关系。
第四方面,本申请另一实施例还提供了一种计算机存储介质,所述计算机存储介质存储有计算机程序,所述计算机程序用于使计算机执行本申请实施例提供的任一超声成像方法和超声图像处理方法。
本申请实施例,将预先选定的血流成像区域划分成多个子区域。按照划分顺序,向各子区域内分别发射超声波,以此得到每个子区域对应的接收线。在确定每个子区域对应的接收线后,按照每条接收线对应的子取样门位置对接收线进行划分,以此得到各子取样门对应的子接收线。通过上述步骤所得到的,各子区域的接收线和各子取样门对应的子接收线构建血流成像区域的四维模型。该四维模型能够描述血流成像区域内的血流在时间、血流速度、成像深度以及接收线这四个维度的变化情况。以此解决相关技术中存在只能对特定的血流区域进行评估和诊断,不能对血流区域进行全面且准确的评估的问题。
本申请的其它特征和优点将在随后的说明书中阐述,并且,部分地从说明书中变得显而易见,或者通过实施本申请而了解。本申请的目的和其他优点可通过在所写的说明书、权利要求书、以及附图中所特别指出的结构来实现和获得。
附图说明
为了更清楚地说明本申请实施例的技术方案,下面将对本申请实施例中所需要使用的附图作简单地介绍,显而易见地,下面所介绍的附图仅仅是本申请的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为根据本申请一个实施例的超声设备硬件配置框图;
图2为根据本申请一个实施例的应用原理示意图;
图3a为根据本申请一个实施例的超声成像方法的流程示意图;
图3b为根据本申请一个实施例的血流成像区域示意图;
图3c为根据本申请一个实施例获取各子区域内接收线的示意图;
图3d为根据本申请一个实施例对子区域接收线进行插值处理的示意图;
图3e为根据本申请一个实施例子区域门频谱图像的示意图;
图3f为根据本申请一个实施例接收线划分子取样门的示意图;
图4为根据本申请一个实施例的超声图像处理方法的流程示意图。
具体实施方式
下面将结合附图对本申请实施例中的技术方案进行清楚、详尽地描述。在本申请实施例的描述中,除非另有说明,″/″表示或的意思,例如,A/B可以表示A或B;文本中的″和/或″仅仅是一种描述关联对象的关联关系,表示可以存在三种关系,例如,A和/或B,可以表示:单独存在A,同时存在A和B,单独存在B这三种情况,另外,在本申请实施例的描述中,″多个″是指两个或多于两个。
在本申请实施例的描述中,除非另有说明,术语″多个″是指两个或两个以上,其它量词与之类似应当理解,此处所描述的优选实施例仅用于说明和解释本申请,并不用于限定本申请,并且在不冲突的情况下,本申请的实施例及实施例中的特征可以相互组合。
为进一步说明本申请实施例提供的技术方案,下面结合附图以及具体实施方式对此进行详细的说明。虽然本申请实施例提供了如下述实施例或附图所示的方法操作步骤,但基于常规或者无需创造性的劳动在方法中可以包括更多或者更少的操作步骤。在逻辑上不存在必要因果关系的步骤中,这些步骤的执行顺序不限于本申请实施例提供的执行顺序。方法在实际的处理过程中或者控制设备执行时,可以按照实施例或者附图所示的方法顺序执行或者并行执行。
本申请的其它特征和优点将在随后的说明书中阐述,并且,部分地从说明书中变得显而易见,或者通过实施本申请而了解。本申请的目的和其他优点可通过在所写的说明书、权利要求书、以及附图中所特别指出的结构来实现和获得。
鉴于相关技术中超声成像都存在着只能对特定的血流区域进行评估和诊断,不能对血流区域进行全面且准确的评估的问题。本申请实施例提供了解决方案。本申请实施例中,能够从多维度描述使血流信息并进行展示。基于超声成像系统单次发射超声波,能够返回多条波束(即,多条接收线)的能力,本申请的发明构思为:基于超声成像系统单次发射超声波后所得接收线的条数和血流成像区域接收线的总条数将血流成像区域划分为若干子区域。依序向每个子区域发射超声波,并接收携带有血流信号的接收线(从该血流信号中能够分析出血流速度)。通过对每个子区域对应的接收线进行频谱分析能够获得血流成像区域内的血流在时间、接收线以及血流速度三个维度的血流状态变化。进一步的,对每条接收线的取样门划分出多个子取样门,并按照每条接收线对应的子取样门的位置将每条接收线进行划分,以此得到各子取样门对应的子接收线。针对每条子接收线进行频谱分析能够获得血流成像区域内的血流在时间、成像深度以及血流速度三个维度的血流变化状态,并进一步结合不同接收线,获得一个接收线维度。故此,本申请根据各子区域的接收线,以及各子取样门对应的子接收线构建出血流成像区域的四维模型,该四维模型能够描述血流成像区域内的血流在时间、血流速度、成像深度以及接收线之间的关联关系。通过上述方法能够从四个维度描述血流变化状态进行多维度展示,使得对血流区域的评估和诊断更加准确和全面。
图1示出了本申请一个实施例提供的超声设备100的结构示意图。下面以超声设备100为例对实施例进行具体说明。应该理解的是,图1所示超声设备100仅是一个范例,并且超声设备100可以具有比图1中所示的更多的或者更少的部件,可以组合两个或多个的部件,或者可以具有不同的部件配置。图中所示出的各种部件可以在包括一个或多个信号处理和/或专用集成电路在内的硬件、软件、或硬件和软件的组合中实现。
图1中示例性示出了根据示例性实施例中超声设备100的硬件配置框图。
如图1所示,超声设备100例如可以包括:处理器110、存储器120、显示单元130和探头140;其中,
探头140,被配置为发射宽波束,并接收各取样门反馈的回波信号;
显示单元130,被配置为显示超声图像;
存储器120被配置为存储用于超声图像所需的数据,可包括软件程序,应用界面数据等;
处理器110,分别与所述探头140以及所述显示单元130相连接,被配置为:
将预先选定的血流成像区域划分为多个子区域;
依序对每个所述子区域分别发射超声波得到每个所述子区域各自的接收线;
对每条接收线分别执行,按照所述接收线对应的子取样门位置,将所述接收线划分为多份,得到各子取样门分别对应的子接收线;
依据各子区域的接收线,以及各子取样门对应的子接收线,构建所述血流成像区域的四维模型,所述四维模型用于描述:时间、血流速度、成像深度以及接收线之间的关联关系。
在一些可能的实施例中,所述处理器执行将预先选定的血流成像区域划分为多个子区域时,被配置为:
基于单次发射超声波后得到的接收线条数和所述血流成像区域的总接收线条数,对所述血流成像区域进行划分操作,得到所述多个子区域;其中,所述多个子区域的个数与所述总接收线条数具有反比关系,且与单次发射超声波后得到的接收线条数具有反比关系。
在一些可能的实施例中,所述处理器还被配置为:
针对不同子区域采用分时检测扫描的方式,在一个扫描周期内依序向各子区域分别发送超声波、重复n个所述扫描周期,获得各子区域在不同扫描周期内的各自的接收线;其中,n为正整数;
针对同一子区域,采用插值处理方式,得到相邻两个扫描周期之间的接收线。
在一些可能的实施例中,每条接收线的多个子取样门中位置相邻的子取样门之间具有重叠区域。
在一些可能的实施例中,针对每条接收线,所述接收线的子取样门的数量与所述接收线的取样门内的总点数成正比、与单个子取样门内设置的点数成反比,且与子取样门之间的重叠区域的大小成反比。
在一些可能的实施例中,所述处理器还被配置为:
接收针对血流成像区域的四维模型的查看指示,其中,所述四维模型用于描述:时间、血流速度、成像深度以及接收线之间的关联关系;
响应所述查看指示,展示所述四维模型中与所述查看指示对象的成像信息。
在一些可能的实施例中,所述查看指示用于从所述四维模型中获取与所述查看指示对应的数据构建所述四维模型中的时间、血流速度、成像深度以及接收线表示四个维度中的任一维度;
所述处理器响应所述查看指示,执行展示所述四维模型中与所述查看指示对象的成像信息时,被配置为:
从所述四维模型中,获取选定的维度对应的其他三个维度的信息,构建三维模型并展示。
在一些可能的实施例中,所述处理器还被配置为:
构建所述四维模型包括:
将所述血流成像区域划分为多个子区域;
依序对每个所述子区域分别发射超声波得到每个所述子区域各自的接收线;
对每条接收线分别执行,按照所述接收线对应的子取样门位置,将所述接收线划分为多份,得到各子取样门分别对应的子接收线;
依据各子区域的接收线,以及各子取样门对应的子接收线,构建所述血流成像区域的四维模型。
在一些可能的实施例中,所述处理器执行将预先选定的血流成像区域划分为多个子区域时,被配置为:
基于单次发射超声波后得到的接收线条数和所述血流成像区域的总接收线条数,对所述血流成像区域进行划分操作,得到所述多个子区域;其中,所述多个子区域的个数与所述总接收线条数具有反比关系,且与单次发射超声波后得到的接收线条数具有反比关系。
在一些可能的实施例中,处理器110被配置为接收针对血流成像区域的四维模型的查看指示,其中,所述四维模型用于描述:时间、血流速度、成像深度以及接收线之间的关联关系;
响应所述查看指示,展示所述四维模型中与所述查看指示对象的成像信息。
在一些可能的实施例中,所述查看指示用于从所述四维模型中获取与所述查看指示对应的数据构建所述四维模型中的时间、血流速度、成像深度以及接收线表示四个维度中的任一维度;
所述处理器110响应所述查看指示,执行展示所述四维模型中与所述查看指示对象的成像信息时,被配置为:
从所述四维模型中,获取选定的维度对应的其他三个维度的信息,构建三维模型并展示。
在一些可能的实施例中,所述装置还被配置为:
构建所述四维模型包括:
将所述血流成像区域划分为多个子区域;
依序对每个所述子区域分别发射超声波得到每个所述子区域各自的接收线;
对每条接收线分别执行,按照所述接收线对应的子取样门位置,将所述接收线划分为多份,得到各子取样门分别对应的子接收线;
依据各子区域的接收线,以及各子取样门对应的子接收线,构建所述血流成像区域的四维模型。
图2为根据本申请一个实施例的应用原理的示意图。其中,该部分可由图1所示超声设备的部分模块或功能组件实现,下面将仅针对主要的部件进行说明,而其它部件,如存储器、控制器、控制电路等,此处将不进行赘述。
如图2所示应用环境中可以包括经由输入输出单元提供的待用户操作的用户界面310、用于显示所述用户界面的显示单元320以及处理器330。
显示单元320可以包括显示面板321、背光组件322。其中,显示面板321被配置为对超声图像进行显示,背光组件322位于显示面板321背面,背光组件322可以包括多个背光分区(图中未示出),各背光分区可以发光,以点亮显示面板321。
处理器330可以被配置为控制背光组件322中各背光分区的背光源亮度,以及控制探头发射宽波束和接收回波信号。
其中处理器330可以包括聚焦处理单元331、波束合成单元332、频谱生成单元333。其中聚焦处理单元331可以被配置为逐一对每个取样门执行聚焦处理,聚焦处理包括:以取样门为宽波束的聚焦位置,根据取样门的发射系数向目标检测区域发射宽波束;并接收每个取样门反馈的回波信号。波束合成单元332被配置为当对目标检测区域的所有取样门完成一轮聚焦处理之后,分别对同一取样门反馈的回波信号进行波束合成,得到扫描信息。频谱生成单元333被配置为基于各取样门的扫描信息进行多普勒成像。
在采用多普勒成像技术对血流信息进行定量检测时,为对血流区域能够进行全面准确的评估,本申请实施例提供的一种超声成像方法的流程示意图,具体如图3a所示,包括:步骤301:将预先选定的血流成像区域划分为多个子区域。
采用多普勒成像技术针对血流进行检测时,具体可如图3b所示,B成像区域为一处人体组织结构;B成像区域中,处于血流成像区域以外的区域为组织细胞等非血液物质;血流成像区域中为人体血液。在针对血流成像区域进行检测时,由于血流成像区域较宽,为便于获取血流成像区域中的血流状态,可以基于对血流成像区域单次发射超声波后得到的接收线条数和该区域的总接收线条数,对血流成像区域进行划分操作。
如,将总接收线条数看做血流成像区域的长度,则单次发送超声波所得到的接收线条数可看做单个子区域的长度。故此,从血流成像区域中划分出的子区域的个数应与单次发射超声波后得到的接收线条数具有反比关系,其与总接收线条数间具有正比关系。考虑到存在总接收线条数与单次发射超声波后得到的接收线条数间不是整数倍的情况(即,无法将血流成像区域均分成若干子区域),此时可将总接收线与单次接收线的比值进行取整后,根据取整结果对血流成像区域进行划分。例如,预设的血流成像区域中共有B条接收线,超声成像系统向血流成像区域单次发射超声波后,对接收到的超声回波信号进行波束合成,得到A条接收线。血流成像区域所划分成若干子区域的个数如公式(1)所示:
Figure BDA0002837110980000131
其中,N为血流成像区域所划分子区域的个数;B为总接收线条数;A为单次发射超声波后得到的接收线条数;floor()为向下取整函数;ceil()为向上取整函数。
在对血流成像区域进行划分时,其划分顺序可根据实际需求自行设置,如从右至左划分、从左至右划分等,本申请对划分顺序不作限定。通过上述公式将血流成像区域划分为N个子区域后,执行步骤302:依序对每个所述子区域分别发射超声波得到每个所述子区域各自的接收线。
在针对血流成像区域内的血流状态进行检测时,需要向血流成像区域内发射超声波,根据该区域所反馈的超声回波进行波束处理,得到该区域对应的接收线。将该接收线进行频谱分析即可得到该区域的血流状态。为提高检测结果的准确性,需要将血流成像区域划分为若干子区域后,获取不同时间下,每个子区域对应的接收线。实施时,可预设扫描周期,针对每个子区域采用分时检测扫描的方式来获取不同时间下,各子区域对应的接收线。具体可如图3c所示,图3c中依据对血流成像区域的划分顺序将子区域排序为子区域1、子区域2和子区域3。在每个扫描周期,按照排序顺序首先向子区域1发射超声波(图中标号1的箭头),接收子区域1反馈的超声回波(超声回波经波束处理即为该区域对应的接收线,图中标号2的箭头)后,向子区域2和子区域3依序重复发射、接收操作。以此得到每个扫描周期内,子区域1、子区域2和子区域3对应的接收线。重复执行操作,直至完成预设的n个执行周期。
由于采用分时检测扫描所获取的检测结果是按照子区域的排序顺序获得的,因而无法在同一时刻采集到每个子区域所对应的接收线。采集结果存在时间间隔会造成对接收线进行频谱分析得到的频谱存在间隙,导致频谱数据稀疏。为解决该问题,可针对同一子区域,将该子区域在相邻扫描周期下所反馈的接收线进行插值处理,得到该区域在相邻扫描周期间的接收线,针对每个子区域实际反馈的接收线和差值处理所得到的接收线进行频谱分析即可得到该区域内,连续时间的血流状态变化。具体可如图3d所示,图3d中V1下属一列中,每个实心圆表示在不同周期,子区域1所反馈的接收线。各实心圆间的空心圆表示经插值处理所得的接收线。V2下属一列中,每个实心圆表示在不同周期,子区域2所反馈的接收线,各实心圆间的空心圆表示经插值处理所得的接收线。V3下属一列中,每个实心圆表示在不同周期,子区域3所反馈的接收线,各实心圆间的空心圆表示经插值处理所得的接收线。以此类推至Vn,可得到N个子区域各自不同扫描周期的接收线。
步骤303:对每条接收线分别执行,按照所述接收线对应的子取样门位置,将所述接收线划分为多份,得到各子取样门分别对应的子接收线。
接收线是由超声回波经过波束处理得到的。超声波发射装置向血流成像区域内发射超声波,接收到预设取样门区域所反馈的超声回波后,对取样门内的超声回波进行采样(采样点的数量与取样门大小以及采样率成正比),得到采样点之后,对采样点进行波束合成得到多条接收线。接收线是由超声回波经过波束处理得到的。超声波发射装置向血流成像区域内发射超声波,接收到预设取样门区域所反馈的超声回波后,对取样门内的超声回波进行采样(采样点的数量与取样门大小以及采样率成正比),得到采样点之后,对采样点进行波束合成得到多条接收线。为了能够多维度的展示血流成像区域内的血流变化状态,考虑到血流成像区域内预设的取样门能够用于表示血流的成像深度,确定接收线在取样门内的总点数(该点数即为上述采样点),并根据实际需求预设子取样门内的点数,基于取样门内的总点数和子取样门内的点数具有反比关系,即可得到该接收线对应的子取样门个数。根据预设的子取样门内点数和子取样门个数可以确定该接收线对应的子取样门位置。将该接收线按照子取样门的位置进行划分,以此得到各子取样门对应的子接收线。针对每条子接收线进行频谱分析,即可得到每个子取样门对应的频谱图像,将每个取样门对应的频谱图像拼接即为该接收线对应的频谱图像。即,从成像深度的维度观测该接收线对应的各子取样门内的血流状态变化,具体如图3e所示。
考虑到子取样门位置划分不妥,会导致相邻子取样门内血流状态变化相差较大的情况,该情况会导致频谱图像中的血流变化曲线产生较大变化,不利于对数据的分析。为避免该情况,可在针对子取样门进行划分时,将位置相邻的子取样门间设置重叠区域。该重叠区域可根据实际情况预设数据重叠率,数据重叠率越大则说明相邻子取样门间的重叠区域越大。在对各子取样门对应的子接收线进行频谱分析时,由于相邻子取样门间存在重叠区域,会使得频谱中的血流变化曲线变化循序渐进,不会产生突兀变化的情况。为了便于理解上述划分流程,具体可如图3f所示,其中子取样门与子取样门3相邻,子取样门2为子取样门1和子取样门3之间的重叠区域。
在针对接收线划分子取样门时,若预先确定取样门内点数是预设子取样门内的点数与数据重叠率乘积的倍数,则可根据下述公式(2)进行划分:
Figure BDA0002837110980000161
其中,n为子取样门的个数;M为取样门内点数;N为预设子取样门内的点数;Ration为数据重叠率,数据重叠率的取值范围在0到1之间。
在针对接收线划分子取样门时,若预先确定取样门内点数不是预设子取样门内的点数与数据重叠率乘积的倍数,则可根据下述公式(3)进行划分:
Figure BDA0002837110980000162
其中,n为子取样门的个数;M为取样门内点数;N为预设子取样门内的点数;Ration为数据重叠率,数据重叠率的取值范围在0到1之间;floor()为向下取整函数。
用于表示血流成像区域内血流速度随时间变化状态的频谱图像是由接收线进行频谱分析所得,对血流成像区域所划分的各子区域对应的接收线进行频谱分析可以得到血流成像区域内,在时间、接收线以及血流速度三个维度上的血流状态变化。而对各子取样门对应的子接收线进行频谱分析可以得到血流成像区域内,在时间、成像深度以及血流速度三个维度上的血流状态变化。故此,可基于时间、血流速度、成像深度以及接收线四个维度构建四维模型,该四维模型能够描述血流成像区域内的血流在上述四个维度之间的关联关系。
步骤304:依据各子区域的接收线,以及各子取样门对应的子接收线,构建所述血流成像区域的四维模型,所述四维模型用于描述:时间、血流速度、成像深度以及接收线之间的关联关系。
针对该四维模型进行查看时,该四维模型可展示出某一维度下,其他三个维度的三维数据关系。针对向四维模型下达的查看指示中可包含上述四个维度中的任一维度信息。四维模型在接收到该查看指示后,将查看指示中的维度信息固定后,构建其余三个维度的三维模型,并展示。以此可以得到血流成像区域在某一维度下,其他三个维度的血流变化状态。
在一些可能的实施例中,四维模型从查看指示中获取到选定四维模型中的接收线维度,四维模型响应于查看指示,选取接收线,根据该接收线在四维模型中对应的时间、血流速度以及成像深度可以接收线位置,血流成像区域内的血流速度随时间和成像深度的变化状态。
在一些可能的实施例中,四维模型从查看指示中获取到选定四维模型中的成像深度,四维模型响应于查看指示,确定成像深度,根据该成像深度在四维模型中对应的时间、血流速度以及接收线可以得到固定成像深度下,血流成像区域内的血流速度随时间和接收线的变化状态;
在一些可能的实施例中,四维模型从查看指示中获取到选定四维模型中的时间维度,四维模型响应于查看指示,确定时间维度,根据该时间维度在四维模型中对应的血流速度、成像深度以及接收线可以得到固定时间下,血流成像区域内的血流速度随成像深度和接收线的变化状态。
图4示出了本申请一个实施例提供的超声图像处理方法,应用于超声设备的流程示意图,包括:
步骤401:接收针对血流成像区域的四维模型的查看指示,其中,所述四维模型用于描述:时间、血流速度、成像深度以及接收线之间的关联关系。
步骤402:响应所述查看指示,展示所述四维模型中与所述查看指示对象的成像信息。
在一些可能的实施例中,所述查看指示用于从所述四维模型中获取与所述查看指示对应的数据构建所述四维模型中的时间、血流速度、成像深度以及接收线表示四个维度中的任一维度;
所述响应所述查看指示,展示所述四维模型中与所述查看指示对象的成像信息,包括:
从所述四维模型中,获取选定的维度对应的其他三个维度的信息,构建三维模型并展示。
在一些可能的实施例中,所述方法还包括:
构建所述四维模型包括:
将所述血流成像区域划分为多个子区域;
依序对每个所述子区域分别发射超声波得到每个所述子区域各自的接收线;
对每条接收线分别执行,按照所述接收线对应的子取样门位置,将所述接收线划分为多份,得到各子取样门分别对应的子接收线;
依据各子区域的接收线,以及各子取样门对应的子接收线,构建所述血流成像区域的四维模型。
在一些可能的实施例中,所述将预先选定的血流成像区域划分为多个子区域,包括:
基于单次发射超声波后得到的接收线条数和所述血流成像区域的总接收线条数,对所述血流成像区域进行划分操作,得到所述多个子区域;其中,所述多个子区域的个数与所述总接收线条数具有反比关系,且与单次发射超声波后得到的接收线条数具有反比关系。
在一些可能的实施方式中,本申请提供的一种视频展示方法的各个方面还可以实现为一种程序产品的形式,其包括程序代码,当程序产品在计算机设备上运行时,程序代码用于使计算机设备执行本说明书上述描述的根据本申请各种示例性实施方式的一种视频展示方法中的步骤。
程序产品可以采用一个或多个可读介质的任意组合。可读介质可以是可读信号介质或者可读存储介质。可读存储介质例如可以是——但不限于——电、磁、光、电磁、红外线、或半导体的系统、装置或器件,或者任意以上的组合。可读存储介质的更具体的例子(非穷举的列表)包括:具有一个或多个导线的电连接、便携式盘、硬盘、随机存取存储器(RAM)、只读存储器(ROM)、可擦式可编程只读存储器(EPROM或闪存)、光纤、便携式紧凑盘只读存储器(CD-ROM)、光存储器件、磁存储器件、或者上述的任意合适的组合。
本申请的实施方式的用于视频展示的程序产品可以采用便携式紧凑盘只读存储器(CD-ROM)并包括程序代码,并可以在电子设备上运行。然而,本申请的程序产品不限于此,在本文件中,可读存储介质可以是任何包含或存储程序的有形介质,该程序可以被指令执行系统、装置或者器件使用或者与其结合使用。
可读信号介质可以包括在基带中或者作为载波一部分传播的数据信号,其中承载了可读程序代码。这种传播的数据信号可以采用多种形式,包括——但不限于——电磁信号、光信号或上述的任意合适的组合。可读信号介质还可以是可读存储介质以外的任何可读介质,该可读介质可以发送、传播或者传输用于由指令执行系统、装置或者器件使用或者与其结合使用的程序。
可读介质上包含的程序代码可以用任何适当的介质传输,包括——但不限于——无线、有线、光缆、RF等等,或者上述的任意合适的组合。
可以以一种或多种程序设计语言的任意组合来编写用于执行本申请操作的程序代码,程序设计语言包括面向对象的程序设计语言一诸如Java、C++等,还包括常规的过程式程序设计语言一诸如″C″语言或类似的程序设计语言。程序代码可以完全地在用户电子设备上执行、部分地在用户设备上执行、作为一个独立的软件包执行、部分在用户电子设备上部分在远程电子设备上执行、或者完全在远程电子设备或服务端上执行。在涉及远程电子设备的情形中,远程电子设备可以通过任意种类的网络——包括局域网(LAN)或广域网(WAN)一连接到用户电子设备,或者,可以连接到外部电子设备(例如利用因特网服务提供商来通过因特网连接)。
尽管已描述了本申请的优选实施例,但本领域内的技术人员一旦得知了基本创造性概念,则可对这些实施例做出另外的变更和修改。所以,所附权利要求意欲解释为包括优选实施例以及落入本申请范围的所有变更和修改。
显然,本领域的技术人员可以对本申请进行各种改动和变型而不脱离本申请的精神和范围。这样,倘若本申请的这些修改和变型属于本申请权利要求及其等同技术的范围之内,则本申请也意图包含这些改动和变型在内。

Claims (8)

1.一种超声成像方法,其特征在于,所述方法包括:
将预先选定的血流成像区域划分为多个子区域;
依序对每个所述子区域分别发射超声波得到每个所述子区域各自的接收线;
对每条接收线分别执行,按照所述接收线对应的子取样门位置,将所述接收线划分为多份,得到各子取样门分别对应的子接收线;
依据各子区域的接收线,以及各子取样门对应的子接收线,构建所述血流成像区域的四维模型,所述四维模型用于描述:时间、血流速度、成像深度以及接收线之间的关联关系;所述将预先选定的血流成像区域划分为多个子区域,包括:
基于单次发射超声波后得到的接收线条数和所述血流成像区域的总接收线条数,对所述血流成像区域进行划分操作,得到所述多个子区域;其中,所述多个子区域的个数与所述总接收线条数具有正比关系,且与单次发射超声波后得到的接收线条数具有反比关系。
2.根据权利要求1所述的方法,其特征在于,所述方法还包括:
针对不同子区域采用分时检测扫描的方式,在一个扫描周期内依序向各子区域分别发送超声波,重复n个所述扫描周期,获得各子区域在不同扫描周期内的各自的接收线;其中,n为正整数;
针对同一子区域,采用插值处理方式,得到相邻两个扫描周期之间的接收线。
3.根据权利要求1所述的方法,其特征在于,每条接收线的多个子取样门中位置相邻的子取样门之间具有重叠区域。
4.根据权利要求3所述的方法,其特征在于,针对每条接收线,所述接收线的子取样门的数量与所述接收线的取样门内的总点数成正比、与单个子取样门内设置的点数成反比,且与子取样门之间的重叠区域的大小成反比。
5.一种超声图像处理方法,其特征在于,所述方法包括:
接收针对血流成像区域的四维模型的查看指示,其中,所述四维模型用于描述:时间、血流速度、成像深度以及接收线之间的关联关系;
响应所述查看指示,展示所述四维模型中查看指示对象的成像信息;其中,所述四维模型是通过下述方式构建的:
基于单次发射超声波后得到的接收线条数和所述血流成像区域的总接收线条数,对所述血流成像区域进行划分操作,得到多个子区域;其中,所述多个子区域的个数与所述总接收线条数具有正比关系,且与单次发射超声波后得到的接收线条数具有反比关系;
依序对每个所述子区域分别发射超声波得到每个所述子区域各自的接收线;
对每条接收线分别执行,按照所述接收线对应的子取样门位置,将所述接收线划分为多份,得到各子取样门分别对应的子接收线;
依据各子区域的接收线,以及各子取样门对应的子接收线,构建所述血流成像区域的四维模型。
6.根据权利要求5所述的方法,其特征在于,所述查看指示用于从所述四维模型中获取与所述查看指示对应的数据构建所述四维模型中的时间、血流速度、成像深度以及接收线四个维度中的任一维度;
所述响应所述查看指示,展示所述四维模型中查看指示对象的成像信息,包括:
从所述四维模型中,获取选定的维度对应的其他三个维度的信息,构建三维模型并展示。
7.一种超声设备,其特征在于,血流成像区域中包括预先设置的多个取样门,且各取样门的偏转角度支持任意角度,所述超声设备包括:
探头,被配置为发射宽波束,并接收各取样门反馈的回波信号;
显示单元,被配置为显示超声图像;
处理器,分别与所述探头以及所述显示单元相连接,被配置为执行如权利要求1-6中任一所述的方法。
8.一种计算机存储介质,其特征在于,所述计算机存储介质存储有计算机程序,所述计算机程序用于使计算机执行如权利要求1-6中任一项所述的方法。
CN202011479604.0A 2020-12-15 2020-12-15 超声成像方法、超声图像处理方法、超声设备及存储介质 Active CN112641468B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011479604.0A CN112641468B (zh) 2020-12-15 2020-12-15 超声成像方法、超声图像处理方法、超声设备及存储介质

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011479604.0A CN112641468B (zh) 2020-12-15 2020-12-15 超声成像方法、超声图像处理方法、超声设备及存储介质

Publications (2)

Publication Number Publication Date
CN112641468A CN112641468A (zh) 2021-04-13
CN112641468B true CN112641468B (zh) 2022-11-04

Family

ID=75354081

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011479604.0A Active CN112641468B (zh) 2020-12-15 2020-12-15 超声成像方法、超声图像处理方法、超声设备及存储介质

Country Status (1)

Country Link
CN (1) CN112641468B (zh)

Citations (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5797846A (en) * 1996-12-30 1998-08-25 General Electric Company Method to control frame rate in ultrasound imaging
CN101011267A (zh) * 2006-01-30 2007-08-08 株式会社东芝 超声波诊断装置、超声波诊断方法及其图像处理程序
CN102613990A (zh) * 2012-02-03 2012-08-01 声泰特(成都)科技有限公司 三维超声频谱多普勒的血流速度及其空间分布显示方法
CN103181789A (zh) * 2011-12-28 2013-07-03 三星麦迪森株式会社 基于决定数据提供矢量多普勒图像的超声波系统和方法
CN103519845A (zh) * 2013-10-12 2014-01-22 清华大学 图像的高帧频数据获取及扫描域的时间配准方法和系统
CN108078590A (zh) * 2018-01-03 2018-05-29 声泰特(成都)科技有限公司 基于超声频谱多普勒的血流动力学可视化方法与系统
CN108433742A (zh) * 2018-03-29 2018-08-24 东南大学 一种便携式模数全深度频移搜索型经颅多普勒检测装置及方法
CN109895763A (zh) * 2018-05-17 2019-06-18 华为技术有限公司 基于超声波雷达的泊车车位检测方法和终端
CN109952135A (zh) * 2016-09-30 2019-06-28 索尼互动娱乐股份有限公司 具有差分渲染和声音定位的无线头戴式显示器
CN110327077A (zh) * 2019-07-09 2019-10-15 深圳开立生物医疗科技股份有限公司 一种血流显示方法、装置及超声设备和存储介质
CN110840484A (zh) * 2019-11-27 2020-02-28 深圳开立生物医疗科技股份有限公司 自适应匹配最优声速的超声成像方法、装置及超声设备

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6398735B1 (en) * 2000-03-07 2002-06-04 Koninklijke Philips Electronics N.V. Detecting a relative level of an ultrasound imaging contrast agent
JP5268374B2 (ja) * 2008-01-25 2013-08-21 株式会社東芝 超音波診断装置、及びその制御方法
JP5231828B2 (ja) * 2008-02-04 2013-07-10 株式会社東芝 超音波診断装置、超音波画像処理装置および超音波画像処理プログラム
US11717270B2 (en) * 2017-06-06 2023-08-08 Shenzhen Mindray Bio-Medical Electronics Co., Ltd. Method, apparatus and system for imaging in ultrasonic scanning
CN110613477B (zh) * 2018-12-29 2024-05-24 深圳迈瑞生物医疗电子股份有限公司 超声成像方法以及超声设备

Patent Citations (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5797846A (en) * 1996-12-30 1998-08-25 General Electric Company Method to control frame rate in ultrasound imaging
CN101011267A (zh) * 2006-01-30 2007-08-08 株式会社东芝 超声波诊断装置、超声波诊断方法及其图像处理程序
CN103181789A (zh) * 2011-12-28 2013-07-03 三星麦迪森株式会社 基于决定数据提供矢量多普勒图像的超声波系统和方法
CN102613990A (zh) * 2012-02-03 2012-08-01 声泰特(成都)科技有限公司 三维超声频谱多普勒的血流速度及其空间分布显示方法
CN103519845A (zh) * 2013-10-12 2014-01-22 清华大学 图像的高帧频数据获取及扫描域的时间配准方法和系统
CN109952135A (zh) * 2016-09-30 2019-06-28 索尼互动娱乐股份有限公司 具有差分渲染和声音定位的无线头戴式显示器
CN108078590A (zh) * 2018-01-03 2018-05-29 声泰特(成都)科技有限公司 基于超声频谱多普勒的血流动力学可视化方法与系统
CN108433742A (zh) * 2018-03-29 2018-08-24 东南大学 一种便携式模数全深度频移搜索型经颅多普勒检测装置及方法
CN109895763A (zh) * 2018-05-17 2019-06-18 华为技术有限公司 基于超声波雷达的泊车车位检测方法和终端
CN110327077A (zh) * 2019-07-09 2019-10-15 深圳开立生物医疗科技股份有限公司 一种血流显示方法、装置及超声设备和存储介质
CN110840484A (zh) * 2019-11-27 2020-02-28 深圳开立生物医疗科技股份有限公司 自适应匹配最优声速的超声成像方法、装置及超声设备

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
Ultrasound‐induced modifications of cytoskeletal components in osteoblast‐like SAOS‐2 cells;Joerg Hauser et al;《JOURNAL OF ORTHOPAEDIC RESEARCH MARCH》;20080827;全文 *
直方图匹配算法在超声弹性成像上的应用研究;何颖妮,邵党国,刘东权;《计算机应用研究》;20130430;第30卷(第4期);全文 *

Also Published As

Publication number Publication date
CN112641468A (zh) 2021-04-13

Similar Documents

Publication Publication Date Title
CN104510499A (zh) 医学超声成像中的剪切波检测
CN113040878B (zh) 超声穿刺针的位置信息处理方法、超声设备及存储介质
CN110353731B (zh) 超声弹性成像方法及系统
CN103984842A (zh) 用于提供诊断信息的诊断帮助设备和方法及其诊断系统
CN111698947B (zh) 多参数组织硬度量化
KR101656146B1 (ko) 초음파 진단 장치 및 그 제어 프로그램
US20200060653A1 (en) Ultrasonic diagnostic device and operating method thereof
US9773496B2 (en) Apparatus and system for adaptively scheduling ultrasound system actions
CN112641468B (zh) 超声成像方法、超声图像处理方法、超声设备及存储介质
US9999404B2 (en) Ultrasonic diagnostic apparatus
WO2018146296A1 (en) Ultrasound evaluation of anatomical features
US20160345937A1 (en) System and method for imaging using ultrasound
US20190043387A1 (en) Systems and methods for providing an interactive demonstration of an ultrasound user interface
JP6731321B2 (ja) 超音波診断装置
EP2915491A1 (en) Method, apparatus, and system for outputting information of blood flow included in region of interest based on selection information
CN110604595B (zh) 一种脂肪肝定量分析方法及脂肪肝定量分析系统
EP2853918B1 (en) Shear wave detection in medical ultrasound imaging
Doblhoff et al. Recognising small image quality differences for ultrasound probes and the potential of misdiagnosis due to undetected side lobes
KR20190121935A (ko) 초음파 변환자 접촉평가 방법 및 그 장치
US20240118245A1 (en) Ultrasonic measurement representation
US12125401B2 (en) Systems and methods for providing an interactive demonstration of an ultrasound user interface
Bradway et al. Real-time GPU implementation of transverse oscillation vector velocity flow imaging
CN106999152B (zh) 超声诊断设备和程序
JP2023101197A (ja) 超音波データ解析装置、超音波検査装置、超音波データ解析方法、超音波検査方法、プログラム、及び記録媒体
US20210338204A1 (en) Ultrasound system and methods for smart shear wave elastography

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