CN104114101B - 自动多普勒脉动周期选择 - Google Patents

自动多普勒脉动周期选择 Download PDF

Info

Publication number
CN104114101B
CN104114101B CN201280069732.5A CN201280069732A CN104114101B CN 104114101 B CN104114101 B CN 104114101B CN 201280069732 A CN201280069732 A CN 201280069732A CN 104114101 B CN104114101 B CN 104114101B
Authority
CN
China
Prior art keywords
cycle
equipment
parameter
stream
fragment
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
CN201280069732.5A
Other languages
English (en)
Other versions
CN104114101A (zh
Inventor
L·古普塔
P·瓦集内帕里
R·S·西索迪亚
G·拉马钱德兰
C·菲尔雄
J·彼得鲁齐洛
A·阿南德
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.)
Koninklijke Philips NV
Original Assignee
Koninklijke Philips Electronics NV
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 Koninklijke Philips Electronics NV filed Critical Koninklijke Philips Electronics NV
Publication of CN104114101A publication Critical patent/CN104114101A/zh
Application granted granted Critical
Publication of CN104114101B publication Critical patent/CN104114101B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/02Measuring pulse or heart rate
    • 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/08Detecting organic movements or changes, e.g. tumours, cysts, swellings
    • A61B8/0891Detecting organic movements or changes, e.g. tumours, cysts, swellings for diagnosis of blood vessels
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/13Tomography
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/46Ultrasonic, sonic or infrasonic diagnostic devices with special arrangements for interfacing with the operator or the patient
    • A61B8/461Displaying means of special interest
    • A61B8/463Displaying means of special interest characterised by displaying multiple images or images and diagnostic data on one display
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/46Ultrasonic, sonic or infrasonic diagnostic devices with special arrangements for interfacing with the operator or the patient
    • A61B8/461Displaying means of special interest
    • A61B8/466Displaying means of special interest adapted to display 3D data
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/46Ultrasonic, sonic or infrasonic diagnostic devices with special arrangements for interfacing with the operator or the patient
    • A61B8/467Ultrasonic, sonic or infrasonic diagnostic devices with special arrangements for interfacing with the operator or the patient characterised by special input means
    • A61B8/468Ultrasonic, sonic or infrasonic diagnostic devices with special arrangements for interfacing with the operator or the patient characterised by special input means allowing annotation or message recording
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/46Ultrasonic, sonic or infrasonic diagnostic devices with special arrangements for interfacing with the operator or the patient
    • A61B8/467Ultrasonic, sonic or infrasonic diagnostic devices with special arrangements for interfacing with the operator or the patient characterised by special input means
    • A61B8/469Ultrasonic, sonic or infrasonic diagnostic devices with special arrangements for interfacing with the operator or the patient characterised by special input means for selection of a region of interest
    • 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
    • A61B8/5223Devices 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
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S15/00Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
    • G01S15/88Sonar systems specially adapted for specific applications
    • G01S15/89Sonar systems specially adapted for specific applications for mapping or imaging
    • G01S15/8906Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques
    • G01S15/8979Combined Doppler and pulse-echo imaging systems
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16HHEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H50/00ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
    • G16H50/30ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for calculating health indices; for individual health risk assessment
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/08Detecting organic movements or changes, e.g. tumours, cysts, swellings
    • A61B8/0866Detecting organic movements or changes, e.g. tumours, cysts, swellings involving foetal diagnosis; pre-natal or peri-natal diagnosis of the baby
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/44Constructional features of the ultrasonic, sonic or infrasonic diagnostic device
    • A61B8/4444Constructional features of the ultrasonic, sonic or infrasonic diagnostic device related to the probe
    • A61B8/4455Features of the external shape of the probe, e.g. ergonomic aspects
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/44Constructional features of the ultrasonic, sonic or infrasonic diagnostic device
    • A61B8/4444Constructional features of the ultrasonic, sonic or infrasonic diagnostic device related to the probe
    • A61B8/4472Wireless probes
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/46Ultrasonic, sonic or infrasonic diagnostic devices with special arrangements for interfacing with the operator or the patient
    • A61B8/467Ultrasonic, sonic or infrasonic diagnostic devices with special arrangements for interfacing with the operator or the patient characterised by special input means
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S7/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/52Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
    • G01S7/52017Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 particularly adapted to short-range imaging
    • G01S7/52079Constructional features
    • G01S7/52084Constructional features related to particular user interfaces

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Public Health (AREA)
  • Medical Informatics (AREA)
  • Pathology (AREA)
  • General Health & Medical Sciences (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Biophysics (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Radiology & Medical Imaging (AREA)
  • Veterinary Medicine (AREA)
  • Molecular Biology (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • Remote Sensing (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Acoustics & Sound (AREA)
  • Hematology (AREA)
  • Vascular Medicine (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • General Physics & Mathematics (AREA)
  • Cardiology (AREA)
  • Physiology (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Computer Graphics (AREA)
  • General Engineering & Computer Science (AREA)
  • Primary Health Care (AREA)
  • Epidemiology (AREA)
  • Databases & Information Systems (AREA)
  • Data Mining & Analysis (AREA)
  • Ultra Sonic Daignosis Equipment (AREA)
  • Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)

Abstract

一种设备(308)被配置用于检查搏动流,用于基于所检查的流来导出谱特征并用于基于所导出的特征来确定选择哪一或哪些脉动周期作为所述流的表示。所选择的周期能够是连续的并且达到预定的周期数量,例如五个。最初可以基于波形异常来滤除受到选择的周期(200),其中,基于诸如波形卡规测量结果和其他类型的特征的参数来判断处于足够数量的连续的组中的存留周期。如果测得参数与在未经最初过滤的谱图周期上的每一各自参数的参数中值的变化不大,那么检测到良好周期(202)。所述技术可以根据用户选择考虑适于具体医疗应用的额外参数。用途包括通过名称正确地识别出动脉。

Description

自动多普勒脉动周期选择
技术领域
本发明涉及搏动流,更具体而言,涉及选择表示所述流的脉动周期。
背景技术
商业双工超声系统被广泛地用于对血管定位,并获得来自血管的流特征。例如,在产科中,存在针对子宫动脉、脐动脉、大脑中动脉、心脏内应用、颈动脉等的应用。双工超声系统将对解剖结构的实时、二维、脉动回波成像的模态与多普勒超声系统的模态组合,从多普勒超声系统能获得来自血管的多普勒频移或速度信息。
在血管应用中使用超声来执行多普勒速度测量需要准确的流参数的计算以生成一致、可重现并且可靠的诊断。
计算流参数的准确度取决于超声医师或医生所选出的周期,其中,人工选择良好的周期。
多普勒检查通常需要很高程度的技能才能获得临床上有用的测量结果。例如,探头相对于血管的正确取向对于确保射束-流角度小于60度是很重要的。当在速度确定中使用大于60度的角时,测量结果中的误差被放大。关于临床超声扫描器的标准工艺流程允许声谱仪操作员使用标准的B模式和彩色流显示器来确定探头相对于血管的取向。之后获得谱多普勒测量结果,因而确保了测得的速度是正确的。
在血管应用中使用超声来执行多普勒速度测量需要有技能型人员。
发明内容
在诸如印度的新兴市场国家,专家短缺限制了超声的可用性和使用性。因此,采集和评估用于临床诊断的多普勒信号的自动化方法(不需要用户解读超声扫描图像)将对作为首要处置提供者的诸如OB/GYN或心脏病专家的非放射科医师有用。
另外,低成本系统对于在新兴市场环境中提供有吸引力的解决方案也是很重要的。当前市售的用于产前检查和分娩的设备是超声扫描器和胎儿监测器/胎心宫缩图(CTG)机。然而,这些设备都相对昂贵。
需要一种低成本的易于使用的解决方案来提供多普勒速度测量,从而为高危妊娠进行筛查和监测。
另外,甚至就在除了产科以外的医疗应用中的双工超声系统以及类似系统而言,良好脉动周期的人工的即可视的选择也需要专业技能并且是一项繁琐而耗时的任务。具体而言,良好的周期是表示血管中的实际血液动力学概况的周期。
而且,关于对什么构成良好的周期的判断在观察者当中有所不同。
尤其是在诸如印度的新兴市场中更是需要自动选择。
发明名称为“Automated Doppler Velocimetry Using a Low-CostTransducer”和“Excitation Schemes for Low-Cost Transducer Arrays”的共同转让的专利申请公开了一种手持的、独立的、基于多普勒的超声探头,其检查面被不太精细地划分成了单独的换能器元件,即,为了使用相对少的单独元件,通过引用将两申请的全部公开内容并入本文。如本文所提到的,探头自动运转,而不必对解剖结构的视觉显示加以解读。
本申请涉及新颖的自动脉动周期选择,其具体应用于刚刚在上文中提到的探头以及诸如双工超声系统的用户交互成像系统。
根据本发明的一方面,一种设备被配置用于检查搏动流,用于基于所检查的流来导出谱特征并且用于基于所导出的特征来确定选择哪一或哪些脉动周期作为所述流的表示。
根据一个子方面,所选择的周期是连续的。
根据另一个子方面,所述选择选出预定数量的周期。该数量能够是五。
根据一个其他子方面,所述设备包括显示器,并且所述设备还被配置用于响应于所述确定来显示所选择的一个或多个周期。
根据一个不同的子方面,所述设备还被配置用于对选择的周期进行操作,以计算临床参数。通常针对选择的五个周期中的每个计算临床参数,并将从这五个周期得到的平均值作为临床参数值。
根据另一子方面,所述设备包括显示器,并且所述设备还被配置用于响应于所述计算来显示所述临床参数。
根据相关子方面,所述流是血管的流。
根据一个补充子方面,所述检查包括接收超声。所述设备还被配置用于从所接收的超声生成受到选择的周期。
在又一子方面中,手持的、独立的诊断装置包含所述设备。
在相关子方面中,所述设备以换能器元件为特征并被配置为不以集合的方式使用这些元件中的任何元件对用于为了执行所述导出而所做的检查的任何射束进行聚焦和导引。
根据一个额外的子方面,所述设备包括用于指定在所述确定中使用的血流参数和/或执行所述确定所针对的医疗应用的用户接口。
根据另一个额外的子方面,所述设备包括用户接口,所述用于接口包括显示器。所述设备还被配置用于:通过选出多个片段而进行所述选择,所述多个片段的每个由所述周期中的多个周期构成;并且用于在所述显示器上显示所述片段,以用于经由所述接口的用户选择。
根据另一个补充子方面,所述设备还被配置用于通过突出显示所显示的片段而区分所显示的片段。
在又一不同子方面中,谱图包括受到所述选择的所述周期。
在另一子方面中,受到所述选择的所述周期分别具有多个指定参数。如果一周期的指定参数中的任何参数偏离能从所述谱图导出的各自的平均值超过了预定量,那么所述确定从所述选择中排除这一周期。
在又一子方面中,所述平均值为中值。
作为一个具体的子方面,具有多个指定参数的周期还被配置用于:a)通过指定参数,使用所述参数与平均值的偏差的绝对值作为指数以形成项;并且b)在所述多个参数上对所述项求和,以得到周期质量量度。
在另一个子方面中,所述指数的底数是自然对数的底数的函数。
备选地或额外地,在另一子方面中,所述确定包括对所述量度求和,以得到周期系列片段质量量度。
在又一不同的子方面中,特征的范例是对周期的多普勒谱波形卡规测量结果。所述确定需要计算所述测量结果的偏差。
在又一其他子方面中,所选择的周期表示随时间变化的频移。
在又一额外的子方面中,将所述设备被配置用于基于一组条件中的至少一个来从受到所述选择的所述周期中滤除一周期,所述组包括:a)在所述周期中存在超过一个峰以及存在超过一个谷这两种情况中的至少一种;b)所述周期的谱带中的频率响应未超过预定阈值;c)所述周期在预定持续时间内未达到其峰收缩值;以及d)所述周期具有小于预定长度的长度。
下文将借助下述附图进一步阐述新颖的、自动的脉动周期选择设备的细节,所述附图并不是按比例绘制的。
附图说明
图1是示出了例如超声探头、包含血管的感兴趣体积和血流波形以及各自的临床多普勒指标的示意图;
图2是对良好的脉动周期以及周期系列片段质量量度进行举例说明的概念图;
图3是对脉动周期的片段的一种可能的视觉显示的示意图,所述片段用作选择的候选;并且
图4是用于选择良好的脉动周期的过程的示范性流程图。
具体实施方式
在对文中提出的关于自动脉动周期选择的内容的描述之前大部分是对上文中提到的共同转让的专利申请中公开的基于多普勒的探头的概述。
图1通过说明性而非限制性范例描绘了超声探头100以及包含血管108、110、112的体积或“感兴趣体积”106。还描绘了血流或者“谱多普勒超声”,波形114(即,谱图包络)以及各自的临床多普勒指标116、118。
可将探头100实施为自动的、手持的、独立的、自包含的超声检查设备。其具有换能器外壳120和手柄122。
在换能器外壳120内,非相控二维换能器阵列124包括换能器元件126,元件的数量由扫描体积和解剖结构确定。尽管数据采集可以单独地由元件126发生,但是如以下更为详细地描述的,多个元件可同时地运转以缩短总的采集时段。
如在图1中通过范例可以看出,元件126的数量为32。因此,就10mm的元件尺寸而言,大约覆盖了6cm×6cm的体积。换能器元件126的超声接收面132与外壳120的前表面130平齐,所述接收面同样用于发射(即,发出)超声。
总共只有32个元件126与常规医疗成像中覆盖同样的6cm×6cm的体积所需的大得多的元件数量形成了鲜明的对比。
在这点上,用于医疗成像的电子聚焦与相控阵列换能器的情况一样需要1/2波长(即,1/2λ)或更小的元件间间隔。用于成像的多普勒超声通常能够处于2×106个周期/秒到4×106个周期/秒(2MHz到4MHz)之间的范围内。超声以大约1540米/秒的速度通过身体软组织传播。波长(即,λ)等于速度除以频率。这里,其为1540m/s除以大约2×106周期/秒=0.8毫米。因此,用于显示器的医疗超声成像将需要小于0.4mm的元件间间隔以及小于(0.4mm)2的元件表面面积,该面积小于0.2mm2。因此,就大约1/2λ的小元件尺寸而言,将需要数千元件126来构建2D阵列,该阵列与图1中看到的阵列类似覆盖6cm×6cm的体积。
图1中的元件的间隔(尺寸)为10mm,如上文所讨论的,其通常超过用于对血管108、110、112所存在的感兴趣体积106进行检查的超声的12λ。
更一般而言,根据文中提出的内容的元件126间隔分开超过1/2λ,尽管元件间间隔128可以是λ、2λ或更大,如上文中所讨论的。对应地,面132的面积至少为0.6平方毫米(mm2),并且其可以更大,例如10mm2、25mm2或者100mm2,如图1所示。
有利地,自动超声设备100不依赖于医疗图像的显示就能达成诊断;但是其反而以包括更少换能器元件以及因此更少信道的阵列为特征。因此,制造成本低,同时借助于自动运转维持了可靠性。在必须以较快的速度执行医疗检查时,甚至可以改善可靠性。自动运转还倾向于减少检查时间,由此减轻工作量,从而使检查更加方便。
在多普勒数据采集期间,元件126以顺序方式被激活(fire),或者以一个或多个组的方式被激活,其中,注意来自一个元件的声学信号不会显著影响同时受到激励的其他信号。对于每个元件126而言,接收时段滞后于发射时段。将多普勒接收栅极对应地放置在接收时段中,从而使得能够从感兴趣体积106内的对应深度采样。
多个用户接口以及包括上面板136、左面板138和右面板140的输入-输出面板处于外壳120的后表面134上,从而使其面对用户。将开启-关闭开关142和音频扬声器面144设置在上面板136中。左面板138框入了功能导航/启动按钮146、显示器148、多普勒功率检测指示器150、胎儿心搏采集指示器152、孕妇心搏采集指示器154、正常血流指示器156以及异常血流指示器158。右面板140包括三个初始化参数输入反馈窗口160、162、164。
阵列124的元件126全部独立地被操作以进行成像。
这与(例如)以集合的方式使用多个单独的换能器元件进行成像或导引射束的相控阵列形成了对比。在相控阵列中,通过使元件的输入和/或输出相对于其他元件适当延迟而执行导引和聚焦。
根据文中提出的内容,同时激活一个组的换能器元件。该组元件继续同时地并且由元件独立地进行成像直到该组的数据采集时段终止。
一种用于按照组进行成像的设备被配置为不以集合的方式使用元件126中的任何元件对成像中使用的射束进行聚焦和导引。作为演示,图1中的换能器元件166、168、170、172每个具有其各自的信号174、176、178、180。信号174、180在发射时滞后于发射信号176、178,由此得到对所产生的超声射束的聚焦和/或导引。如在图1中通过“X”182指示的,探头100并不是针对这样的协议实施的。类似地,在接收时,也不向元件166、168、170、172差异性地施加延迟。
血流波形114是随时间变化的多普勒频移的图,并且其由此指示了随时间变化的血流速度。
诸如搏动指数(PI)116和阻力指数(RI)118的临床多普勒指标是血液搏动性的与角度无关的多普勒度量。标注血流波形114的符号S、D和C分别表示收缩期峰值频移、舒张末期频移以及一个心搏周期的长度。另一个常用的临床多普勒指标是收缩期/舒张期比S/D。符号S、D和C是多普勒谱波形的卡规测量结果。这些以及诸如临床多普勒指标的其他临床或血流参数是搏动流的谱特征的范例,能够基于其判断周期的质量。另一个谱特征是谱带内的频率响应,如果其足够低则指示周期质量不够。类似地,如果一周期没有在给定的时段内取得其峰收缩值或者其完整周期长度,那么认为其质量不够。
探头100能够利用上文识别的多普勒指标来识别血管以及评估血流的正常性。
对由探头100在感兴趣体积106中找到的血管108-112进行分类和分析所涉及的信号处理以及有关探头及其使用的更多细节已在上文提到的共同转让的专利申请中公开。
图2例证了多普勒信号的谱图中的良好脉动周期200以及用于从谱图选择一系列良好脉动周期的周期系列片段质量量度202。所述多普勒信号是已经从载波频率(即,要从血流产生回波的脉动的载波频率)中提取出来的。使用正交解调器进行所述提取。对来自所述解调器的输出进行带通滤波,以去除可能源自于诸如血管壁和软组织的缓慢移动结构的低频信号以及具有超出特定水平的频率的信号。所得到的多普勒信号指示多普勒频移以及血流速度。
对应于单次心搏的典型的良好脉动周期200在自从该周期开始的可接受的持续时间内具有处于两个谷206、208之间的显著峰204。可接受的持续时间在孕妇动脉和胎儿动脉中是不同的。孕妇动脉的持续时间能够在0.6秒到1.5秒之间变化,而对于胎儿动脉而言,其能够是从0.3秒到1.0秒。还要检查没有局部峰204或谷208的情况。因此,每个脉动周期200优选具有单个峰204和单个谷208。
由于在周期200期间频移的升高和降低分别表示血流速度的升高和降低,所以良好的周期将呈现平稳的升高和平稳的降低。周期200的下降部分(即,从峰到谷)经历另一测试。检查连接峰204和紧随该峰的谷208的线210是否与波形114交叉至少一次。图2示出了交叉点211。没有这样的交叉点211指示周期200不能充分地表示搏动血流。
而且,谱信息中的一些可能并不是可连续解读的,因为目标血管108-112可能由于患者或超声医师的运动或者由于血管108-112自身的运动而未处于样本体积内。如果超声探头100没有与血管108-112成适当的角度,或者由于来自相邻组织或者内置电子器件的噪声的原因,脉动周期200中的一些可能不具有用于解读的良好质量。
结果,谱图的一些部分可能很弱,并且在屏幕上看起来像是消退的或者缺失的,这指示处理输入多普勒信号的快速傅里叶变换(FFT)的相关频率样本的幅度低。在谱图轮廓的形状中观察到的其他类型的异常包括弱峰值强度、不完整的周期200、缺少峰204和尖高的峰。
基于这些考虑,最初滤除周期200中的一些。存留的周期200的谱特征被提取,并被用于对五个连续的存留周期的片段进行评分。依赖于五个周期是可接受的临床实践。本发明人凭经验发现要审慎地选择一系列连续的周期200。在图2中被示为Em的周期系列片段质量量度202用于评分。其为五个周期质量量度212的和201,其中,一个质量量度针对五个连续周期200中的每个。每个周期质量量度212包括三个项214,尽管针对除了产科之外的其他医疗应用也能够增加额外的项。每个项214具有底数216和指数218。底数216是通常被称为“e”的项目(即,自然对数(ln)的底数)的函数。在图2的范例中,所述函数是恒等函数。底数216能够是任何超过一的值,例如,“e”或任何更大的数。指数218是偏差222的绝对值220。偏差222是关对周期200的多普勒谱波形卡规测量结果224与该测量结果在未被作为异常排除的谱图的周期上的中值226之间的偏差。符号Si、Ci和Di表示周期i的收缩期、周期长度和舒张期。指数218是收缩期、周期长度和舒张期与对应的中值226的偏差222的各自的绝对值220。所述中值被用作表示值的度量,因为所述中值相对于极端观测值而言是鲁棒的,并且还充当某种类型的平均。项214的指数性质确保了即使单个周期200偏离了中值226的代价也会招致高惩罚。Em的方程是可扩展的,以包括针对其他医疗应用的其他项214。例如,颈动脉的特征是谱展宽,其为最大频率包络中的峰值与最小频率包络中的峰值之间的距离。最大频率包络是相对于血流中的最大频率应用于谱图的包络。具体而言,由于血液的粘滞度的原因,从横截面的角度上来讲血液以不同的速度流经动脉,其在邻近动脉壁的外围附近是滞后的。对于心脏应用而言,例如,对于那些涉及二尖瓣、三尖瓣和主动脉瓣的应用而言,所述特征包括等容量(isovolume)弛豫时间、等容量加速时间、收缩期早期的加速时间和舒张期早期的减速之间。在肾脏应用中,肾动脉的谱展宽能够充当特征。
图3描绘了一种实施方式,其包括对脉动周期200的片段304、306的视觉显示302,所述片段在根据Em评分之后充当供用户选择的候选。能够(但不需要)将这一实施方式中的超声探头308设计为独立的设备。探头308无线地或者有线地可连接至或被连接至处理器310,处理器310转而被连接至用户接口312。用户接口312具有显示屏314和用户可启动控制316。屏幕314的触敏特征能够被包含于控制316当中,如其他导航和选择设备(例如鼠标、按钮、滑键、旋钮和跟踪球)所能够的。
显示302从左到右跨越显示屏314连续移动。当前处于屏幕上的任何片段304、306都是可选择的,例如,通过触摸触摸屏进行选择。在滚动显示中,两个片段304、306之间的点线320表示时间上的已经排除了周期200或高评分(并且因此不太期望的)片段的位置。通过亮度、颜色或者如这里所示的下划线324、328对所显示的片段304、306突出显示。连同突出显示一起,可以在每个片段304、306跨越屏幕314平移的过程中,使每个片段附有每个周期200的卡规测量结果224和/或片段评分。当选择片段304、306时,所述片段中的周期200的多普勒参数116、118被计算,并出现在屏幕上。
作为使片段滚动的替代方案,所述系统能够自动选择所述片段或具有最低评分的(一个或多个)片段304、306,计算各自的参数116、118并显示所选择的片段和参数。
图4通过说明性而非限制性范例提供了一种用于选择良好的脉动周期的过程400。可以是临床医师、助产士、全科医生、产科医生/妇科医生或胎儿放射科医师的用户或操作员选择诸如产科的医疗应用,或者从一组参数中选择一个或多个参数(步骤S402)。在短时段期间从所接收的超声中捕获多普勒信号(步骤S404)。计算谱图(步骤S406)。
如果在考虑中的当前周期200的特定谱带上缺乏充分的频率响应,如从功率阈值来判断的(步骤S408),那么将当前周期从选择中排除(步骤S410)。
在任一情况下,如果有更多的来自在步骤S404中捕获的数据的周期200(步骤S412),那么处理下一周期(步骤S414),并且处理返回至步骤S408,并以该下一周期充当当前周期。
当没有更多需要考虑的周期200时(步骤S412),没有被以集合方式排除的周期构成了在已经排除的周期的位置处具有缺口的谱图。因此,具有缺口的谱图通常具有其内的周期没有被排除的一个或多个时间部分。通过任何已知的且适当的方法针对单独地包括五个或更多连续周期并因此具有供给五个连续的良好周期200的片段的潜力的(一个或多个)谱图时间部分来计算包络(步骤S415)。在颁发给Zhang的美国专利No.7611467中描述了一种用于从谱图提取包络的方法,通过引用将该专利的全部公开内容并入本文。
周期过滤继续。例如,对当前周期200做出卡规测量结果224,所述测量结果既针对其时间长度又针对直到其峰值的持续时间。
如果任一测量结果224超过了各自的预定正常范围(S416),那么排除当前周期200(步骤S418)。否则,如果范围两者都没有被超过,那么针对在当前周期中是否存在不只一个波峰做出询问(步骤S420),并针对在当前周期中是否存在不只一个波谷做出询问(步骤S422)。对于任一情况,都将当前周期排除(步骤S418)。
在任何情况下,如果有更多的周期可用于考虑(步骤S424),那么处理返回至步骤S416,并以下一周期充当当前周期(步骤S426)。
当没有更多的周期要进行上文描述的过滤流程时(步骤S424),计算所有存留周期上的中值226。这一操作是针对每个参数进行的,这样的直接选择的或者隐式地选择的卡规测量224是步骤S402(步骤S428)。所述计算是计算Em(周期系列片段质量量度202)的初始步骤,如图2所示。
对于从上文描述的过滤中剩余的周期200中的五个连续周期的当前片段304、306而言,进行有关所述片段的每个周期的偏差222的询问。具体而言,针对参数与其的中值226的偏差222的绝对值220是否超过了预定偏差阈值做出询问。这一询问是针对所选择的参数中的每个做出的(步骤S430)。可以将阈值设置为检测与各自的中值226的25%或更高的偏差。如果超过了偏差阈值(步骤S430),从而指示当前周期不可用,那么处理逐周期地移过偏离周期200(步骤S432)。
如果尽管存在移动,当前仍有在步骤S424之前的过滤中存留下来的五个或更多的周期200能够形成片段(S434),那么返回至步骤S430,并以下一片段304、306充当当前片段(S435)。
如果另一方面可获得的周期200少于五个(S434),那么返回至步骤S404。
如果对于当前片段304、306而言没有超过偏差阈值(步骤S430),那么对所述片段评分,从而得到周期系列片段质量量度202(步骤S436)。
如果期望更多的片段304、306(步骤S438),那么处理通过支路返回至步骤S434。
否则,如果不期望获得另外的经评分的良好片段304、306(步骤438),并且用户将从这些片段中做出选择(步骤S440),那么如上文描述的所述片段被显示并被突出显示以供选择(步骤S442)。片段304、306在屏幕上任选地附有卡规测量结果224和片段评分。之后,针对由用户选择的片段304、306计算诸如多普勒指标的参数(步骤S444),并将所述片段和指标都显示到屏幕314上(步骤S446)。另一方面,如果不期望获得另外的经评分的良好片段304、306(步骤S438),并且选择是自动的(步骤S440),那么选择具有最低评分的片段(步骤S448)。类似地,之后针对所选择的片段304、306计算参数(步骤S450),并显示所述片段和参数(步骤S452)。所述计算可以涉及对所选择的片段304、306的五个周期200的结果求平均,将平均值包含在显示屏314上。
借助卡规测量结果对来自所选择的周期200当中的周期执行操作,以计算临床参数,这在本新颖技术中全部是自动完成的。更一般而言,所述计算通常涉及关于每一周期200的这样的测量以及对结果求平均。
在探头100、308内,充当文中提出的设备的控制电路(未示出)能够采取一个或多个集成电路(IC)的形式。或者,所述一个或多个IC能够被配置为安装到诸如超声双工扫描器的现有装置中。
根据文中提出的技术,用于对所述设备(即,(一个或多个)IC)、探头或双工系统进行操作的信号能够是内部形成的,能够是通过改变施加至所述设备的有线输入端或施加至天线的电流而形成的,所述天线用于所述信号的无线发射并用于由所述设备的接收天线进行的接收。
一种设备被配置用于检查搏动流,用于基于所检查的流来导出谱特征并用于基于所导出的特征来确定选择哪一或哪些脉动周期作为所述流的表示。所选择的周期能够是连续的并且达到预定的周期数量,例如五个。最初可以基于波形异常来滤除受到选择的周期,其中,基于诸如波形卡规测量结果和其他类型的特征的参数来判断处于足够数量的连续的组中的存留周期。如果测得参数与在未经最初过滤的谱图周期上的每一各自参数的参数中值的变化不大,那么检测到良好周期。所述技术可以根据用户选择考虑适于具体医疗应用的额外参数。用途包括通过名称正确地识别出动脉。
尽管能够将本发明的方法有利地应用于针对人或动物对象的医疗诊断的提供中,但是本发明的范围不限于此。更宽泛地来说,文中公开的技术涉及有效地找到体内的、体外的或者离体的身体组织中的血管并使其受到改进的流体流分析。文中提出的内容涉及选择谱多普勒波形的良好周期,所选择的周期表示血流,其用于基于所选择的周期的分析结果提出临床诊断。所述技术对于通过名称准确地识别出动脉尤其有用,如相同发明人的所共同拥有的发明名称为“Automatic Blood Vessel Identification by Name”的专利中所述。除了胎儿健康状况评估之外,文中提出的技术的应用还包括颈动脉和肾动脉筛查、用于检测外周动脉疾病(PAD)的踝臂指数(ABI)测量、经颅和心脏检查、外伤或其他出血中的血液检查。
尽管已经在附图和前面的描述中详细说明和描述了本发明,但这样的说明和描述被认为是说明性或示范性的而非限制性的;本发明不限于公开的实施例。
例如,与周期峰的形状相关的波形异常能够被检测并被用于周期的滤除,例如,步骤S418。
通过研究附图、说明书和权利要求书,本领域的技术人员在实践所主张的本发明时能够理解和实现所公开的实施例的其他变型。在权利要求中,词语“包括”不排除其他元件或步骤,并且,量词“一”或“一个”并不排除多个。权利要求中的任何附图标记不得被解释为对范围的限制。
计算机程序能够被暂时地、临时地存储,或者可以在诸如光存储介质或固态介质的适当的计算机可读介质上被存储较长的一段时间。这样的介质只是从并非瞬态传播信号的意义上来讲是非暂态的,但是其包括其他形式的计算机可读介质,例如寄存器存储器、处理器高速缓存以及RAM。
单个处理器或其他单元可以实现权利要求书中记载的若干项目的功能。在互不相同的从属权利要求中记载了特定措施并不指示不能有利地使用这些措施的组合。

Claims (15)

1.一种用于选择脉动周期的设备(308),其被配置用于检查搏动流(S404),用于基于所检查的流来导出多个谱特征并且用于基于所导出的特征来确定选择谱图(114)中的哪一或哪些脉动周期作为所述流的表示,受到所述选择的所述周期分别具有多个指定参数,如果一周期的所述指定参数中的任何参数偏离(222)从所述谱图导出的各自的平均值超过了预定的量,那么所述确定将来自所述周期当中的这一周期从所述选择中排除。
2.如权利要求1所述的设备,所选择的周期(200)是连续的。
3.如权利要求1所述的设备,所述选择选出预定数量的周期。
4.如权利要求3所述的设备,所述数量为五。
5.如权利要求1所述的设备,包括显示器(314),并且所述设备还被配置用于响应于所述确定来显示已经选择的所述周期当中的周期。
6.如权利要求1所述的设备,还被配置用于对来自已经选择的所述周期当中的周期进行操作,以计算临床参数。
7.如权利要求6所述的设备,包括显示器,并且所述设备还被配置用于响应于所述计算来显示所述临床参数(116、118)。
8.如权利要求1所述的设备,包括用户接口,所述用户接口包括显示器,所述设备还被配置用于:通过选出多个片段(304、306)而进行所述选择,所述多个片段的每个由所述周期中的多个周期构成;并且用于在所述显示器上显示所述片段,以用于经由所述接口对来自所述片段中当中的片段的用户选择。
9.如权利要求8所述的设备,还被配置用于通过对所显示的片段(324、328)进行突出显示而区分所显示的片段(324、328)。
10.如权利要求1所述的设备,所述平均值为中值(226)。
11.一种用于选择脉动周期的设备(308),其被配置用于检查搏动流(S404),用于基于所检查的流来导出多个谱特征并且用于基于所导出的特征来确定选择哪一或哪些脉动周期作为所述流的表示,来自所述周期当中的一周期具有多个指定参数,所述设备还被配置用于:a)通过指定参数,使用所述参数与平均值的偏差的绝对值作为指数(218)来形成项;并且b)在所述多个指定参数上对所述项求和以得到周期质量量度。
12.如权利要求11所述的设备,所述指数的底数(216)为自然对数的底数的函数。
13.如权利要求11所述的设备,周期系列片段包括分别具有所述多个指定参数的多个周期,所述周期在所述多个周期当中,所述设备还被配置用于关于所述多个周期进行a)和b),所述确定包括对所述量度进行求和以得到周期系列片段质量量度(202)。
14.如权利要求1所述的设备,还被配置用于如果候选周期的谱带中的频率响应未超过预定阈值,则从受到所述选择的所述周期当中滤除所述候选周期。
15.一种用于选择脉动周期的设备,其被配置用于检查搏动流(S404);用于从搏动流导出多个谱特征(S408、S416、S422)并基于所导出的特征来确定选择所述周期中的哪些周期作为所述流的表示;以及用于如果一周期在预定的持续时间内未达到其峰收缩值,则滤除受到所述选择的所述周期中的这一周期。
CN201280069732.5A 2011-12-16 2012-12-06 自动多普勒脉动周期选择 Active CN104114101B (zh)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
US201161576630P 2011-12-16 2011-12-16
US61/576,630 2011-12-16
PCT/IB2012/057033 WO2013088314A1 (en) 2011-12-16 2012-12-06 Automated doppler pulse cycle selection

Publications (2)

Publication Number Publication Date
CN104114101A CN104114101A (zh) 2014-10-22
CN104114101B true CN104114101B (zh) 2016-08-24

Family

ID=47563558

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201280069732.5A Active CN104114101B (zh) 2011-12-16 2012-12-06 自动多普勒脉动周期选择

Country Status (8)

Country Link
US (1) US20140358000A1 (zh)
EP (1) EP2790588B1 (zh)
JP (1) JP6114757B2 (zh)
CN (1) CN104114101B (zh)
BR (1) BR112014014413B1 (zh)
IN (1) IN2014CN04916A (zh)
RU (1) RU2674241C2 (zh)
WO (1) WO2013088314A1 (zh)

Families Citing this family (17)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US11678808B2 (en) 2009-03-13 2023-06-20 Bluedop Medical, Ltd. Haemodynamic data estimation apparatus and method of use
GB0904435D0 (en) 2009-03-13 2009-04-29 King David H Haemodynamic data estimation
JP6265909B2 (ja) 2011-12-16 2018-01-24 コーニンクレッカ フィリップス エヌ ヴェKoninklijke Philips N.V. 名前による自動血管識別
WO2015011585A1 (en) * 2013-07-24 2015-01-29 Koninklijke Philips N.V. Non-imaging two dimensional array probe and system for classifying carotid stenosis
WO2015011599A1 (en) * 2013-07-24 2015-01-29 Koninklijke Philips N.V. Method for aligning spatially different subvolumes of ultrasonic data of a blood vessel
WO2015130705A1 (en) 2014-02-25 2015-09-03 Icu Medical, Inc. Patient monitoring system with gatekeeper signal
EP2989986B1 (en) * 2014-09-01 2019-12-18 Samsung Medison Co., Ltd. Ultrasound diagnosis apparatus and method of operating the same
WO2017070120A1 (en) 2015-10-19 2017-04-27 Icu Medical, Inc. Hemodynamic monitoring system with detachable display unit
US11660063B2 (en) * 2015-11-18 2023-05-30 Bluedop Medical, Ltd. System for determining peripheral artery disease and method of use
EP3189776A1 (en) * 2016-01-08 2017-07-12 Koninklijke Philips N.V. An apparatus and method for generating fetal heart rate data
EP3463094B1 (en) * 2016-06-07 2020-10-21 Viewcare Technologies 1 ApS Method and system for measuring a central pulse wave velocity in a pregnant woman
US11717255B2 (en) 2016-08-05 2023-08-08 Cimon Medical As Ultrasound blood-flow monitoring
EP3554379B1 (en) 2016-12-15 2022-09-07 Koninklijke Philips N.V. Prenatal ultrasound imaging
FR3060966B1 (fr) * 2016-12-23 2019-05-31 Azoth Systems Dispositif de mesure du flux sanguin
US11690590B2 (en) * 2017-11-13 2023-07-04 Edwards Lifesciences Corporation Non-invasive heart valve screening device and method
WO2019135096A1 (en) 2018-01-02 2019-07-11 Bluedop Medical, Ltd. System for determining peripheral artery disease and method of use
JP7281210B2 (ja) * 2018-02-07 2023-05-25 サイモン メディカル アーエス 超音波血流モニタリング

Family Cites Families (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4476874A (en) * 1982-06-01 1984-10-16 Sri International Ultrasonic imaging with volume flow measuring method and apparatus
US5628321A (en) * 1995-12-18 1997-05-13 Diasonics Ultrasound, Inc. Processing velocity information in an ultrasonic system
US6423006B1 (en) * 2000-01-21 2002-07-23 Siemens Medical Solutions Usa, Inc. Method and apparatus for automatic vessel tracking in ultrasound systems
WO2007085999A1 (en) * 2006-01-27 2007-08-02 Koninklijke Philips Electronics N.V. Automatic ultrasonic doppler measurements
CN101378699B (zh) * 2006-02-03 2012-06-20 皇家飞利浦电子股份有限公司 用于测量或检测非正弦周期性血流行为的超声方法和装置
US8323198B2 (en) * 2007-01-29 2012-12-04 Siemens Medical Solutions Usa, Inc. Spatial and temporal alignment for volume rendering in medical diagnostic ultrasound
WO2009013686A2 (en) * 2007-07-26 2009-01-29 Koninklijke Philips Electronics, N.V. Systems and methods for automated image selection in doppler ultrasound imaging systems
WO2009158399A1 (en) * 2008-06-26 2009-12-30 Verasonics, Inc. High frame rate quantitative doppler flow imaging using unfocused transmit beams
RU2407008C1 (ru) * 2009-06-04 2010-12-20 Федеральное государственное учреждение "Ивановский научно-исследовательский институт материнства и детства имени В.Н. Городкова Федерального агентства по высокотехнологичной медицинской помощи" Способ морфологической диагностики степени тяжести гестоза
US8394025B2 (en) * 2009-06-26 2013-03-12 Uab Vittamed Method and apparatus for determining the absolute value of intracranial pressure
CN103442649B (zh) * 2010-12-22 2015-11-25 皇家飞利浦电子股份有限公司 使用低成本换能器的自动多普勒速度测量法

Also Published As

Publication number Publication date
BR112014014413B1 (pt) 2022-09-06
BR112014014413A2 (pt) 2017-06-13
BR112014014413A8 (pt) 2017-06-13
CN104114101A (zh) 2014-10-22
JP6114757B2 (ja) 2017-04-12
WO2013088314A1 (en) 2013-06-20
US20140358000A1 (en) 2014-12-04
RU2674241C2 (ru) 2018-12-05
IN2014CN04916A (zh) 2015-09-18
EP2790588B1 (en) 2019-11-13
RU2014129043A (ru) 2016-02-10
JP2015500120A (ja) 2015-01-05
EP2790588A1 (en) 2014-10-22

Similar Documents

Publication Publication Date Title
CN104114101B (zh) 自动多普勒脉动周期选择
US11635514B2 (en) Imaging methods and apparatuses for performing shear wave elastography imaging
JP6994494B2 (ja) エラストグラフィ測定システム及びその方法
JP7132925B2 (ja) イメージングデータ及び組織情報を導出するためのニューラルネットワークを用いた超音波イメージングシステム
JP6265909B2 (ja) 名前による自動血管識別
CN103442649B (zh) 使用低成本换能器的自动多普勒速度测量法
JP6547612B2 (ja) 画像処理装置、画像処理方法、および、画像処理装置を備える超音波診断装置
US20140243673A1 (en) Excitation schemes for low-cost transducer arrays
EP3570754A1 (en) Method and apparatus to characterise non-invasively images containing venous blood vessels
JP2016531622A5 (zh)
CN109069121A (zh) 用于ctg超声换能器的定位支持和胎儿心率配准支持
EP3713497A1 (en) Ultrasonic pulmonary assessment
CN107690312A (zh) 超声成像装置
WO2021212578A1 (zh) 粘弹性测量方法、弹性测量方法和超声测量系统
JP5443781B2 (ja) 超音波診断装置
WO2007072720A1 (ja) 医用画像診断装置および生体組織の同定方法
CN111093522B (zh) 剪切波弹性成像方法及超声成像设备
JP2020199162A (ja) 超音波診断装置、および、超音波診断装置の制御方法
Mesin et al. METHOD AND APPARATUS TO CHARACTERISE NON-INVASIVELY IMAGES CONTAINING VENOUS BLOOD VESSELS
CN116919462A (zh) 超分辨率超声成像方法和超声成像系统

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant