CN103415259B - 计算至少两种组织类型中超声的速度 - Google Patents
计算至少两种组织类型中超声的速度 Download PDFInfo
- Publication number
- CN103415259B CN103415259B CN201280011510.8A CN201280011510A CN103415259B CN 103415259 B CN103415259 B CN 103415259B CN 201280011510 A CN201280011510 A CN 201280011510A CN 103415259 B CN103415259 B CN 103415259B
- Authority
- CN
- China
- Prior art keywords
- ultrasonic
- medical image
- subject
- ultrasound
- kinds
- 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.)
- Expired - Fee Related
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/08—Detecting organic movements or changes, e.g. tumours, cysts, swellings
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/103—Detecting, measuring or recording devices for testing the shape, pattern, colour, size or movement of the body or parts thereof, for diagnostic purposes
- A61B5/107—Measuring physical dimensions, e.g. size of the entire body or parts thereof
- A61B5/1075—Measuring physical dimensions, e.g. size of the entire body or parts thereof for measuring dimensions by non-invasive methods, e.g. for determining thickness of tissue layer
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/08—Detecting organic movements or changes, e.g. tumours, cysts, swellings
- A61B8/0825—Detecting organic movements or changes, e.g. tumours, cysts, swellings for diagnosis of the breast, e.g. mammography
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/08—Detecting organic movements or changes, e.g. tumours, cysts, swellings
- A61B8/0858—Detecting organic movements or changes, e.g. tumours, cysts, swellings involving measuring tissue layers, e.g. skin, interfaces
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/40—Positioning of patients, e.g. means for holding or immobilising parts of the patient's body
- A61B8/406—Positioning of patients, e.g. means for holding or immobilising parts of the patient's body using means for diagnosing suspended breasts
-
- 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
- A61B8/4483—Constructional features of the ultrasonic, sonic or infrasonic diagnostic device characterised by features of the ultrasound transducer
- A61B8/4494—Constructional features of the ultrasonic, sonic or infrasonic diagnostic device characterised by features of the ultrasound transducer characterised by the arrangement of the transducer elements
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61N—ELECTROTHERAPY; MAGNETOTHERAPY; RADIATION THERAPY; ULTRASOUND THERAPY
- A61N7/00—Ultrasound therapy
- A61N7/02—Localised ultrasound hyperthermia
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01H—MEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
- G01H5/00—Measuring propagation velocity of ultrasonic, sonic or infrasonic waves, e.g. of pressure waves
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B2562/00—Details of sensors; Constructional details of sensor housings or probes; Accessories for sensors
- A61B2562/02—Details of sensors specially adapted for in-vivo measurements
- A61B2562/0204—Acoustic sensors
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B2562/00—Details of sensors; Constructional details of sensor housings or probes; Accessories for sensors
- A61B2562/04—Arrangements of multiple sensors of the same type
- A61B2562/046—Arrangements of multiple sensors of the same type in a matrix array
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/05—Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves
- A61B5/055—Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves involving electronic [EMR] or nuclear [NMR] magnetic resonance, e.g. magnetic resonance imaging
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
- A61B6/50—Clinical applications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
- A61B6/52—Devices using data or image processing specially adapted for radiation diagnosis
- A61B6/5211—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data
- A61B6/5229—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data combining image data of a patient, e.g. combining a functional image with an anatomical image
- A61B6/5247—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data combining image data of a patient, e.g. combining a functional image with an anatomical image combining images from an ionising-radiation diagnostic technique and a non-ionising radiation diagnostic technique, e.g. X-ray and ultrasound
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/58—Testing, adjusting or calibrating the diagnostic device
- A61B8/585—Automatic set-up of the device
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61N—ELECTROTHERAPY; MAGNETOTHERAPY; RADIATION THERAPY; ULTRASOUND THERAPY
- A61N7/00—Ultrasound therapy
- A61N2007/0078—Ultrasound therapy with multiple treatment transducers
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61N—ELECTROTHERAPY; MAGNETOTHERAPY; RADIATION THERAPY; ULTRASOUND THERAPY
- A61N7/00—Ultrasound therapy
- A61N2007/0086—Beam steering
- A61N2007/0095—Beam steering by modifying an excitation signal
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Animal Behavior & Ethology (AREA)
- Veterinary Medicine (AREA)
- Public Health (AREA)
- Physics & Mathematics (AREA)
- General Health & Medical Sciences (AREA)
- Engineering & Computer Science (AREA)
- Biomedical Technology (AREA)
- Pathology (AREA)
- Biophysics (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Heart & Thoracic Surgery (AREA)
- Radiology & Medical Imaging (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Medical Informatics (AREA)
- Gynecology & Obstetrics (AREA)
- Dentistry (AREA)
- Oral & Maxillofacial Surgery (AREA)
- General Physics & Mathematics (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
一种医学设备(400),包括用于采集描述沿至少两条路径(606、1114)的超声的速度(304、306)的超声数据(476)的超声发射器(444、602)和接收器(446、604)系统(600)。所述医学设备还包括用于采集医学图像数据(468)的医学成像系统(402)和包含使所述处理器采集(100、200)医学图像数据的指令(490、492、494、496、498、500、502、504)的存储器(464)。所述指令还使所述处理器采集(102、202)超声数据。所述指令还使所述处理器将所述医学图像数据分割(104、204)成至少两种组织类型(434、436、610、612、708、710)。所述指令还使所述处理器确定(106、206)对应于所述受试体内的至少两条路径的至少两个距离。所述指令还使所述处理器计算(108、208)至少两种组织类型中的超声的速度。
Description
技术领域
本发明涉及医学成像,尤其涉及使用超声和医学成像来计算组织中超声的速度。
背景技术
可以使用来自聚焦超声换能器的超声有选择地处置身体内部的区域。超声波是作为高能机械振动传输的。这些振动在受到阻尼时将诱发组织发热,而且还能够导致空化作用。所述组织发热和空化作用两者都能够用于在临床环境中对组织进行破坏。然而,利用超声对组织进行加热比空化作用更易于控制。能够使用超声处置对组织进行消融,从而有选择地杀死癌细胞的区域。这种技术已经被应用于子宫肌瘤的处置,并且其降低了子宫切除术的必要性。能够以较低的功率或者以脉冲模式使用超声有选择地向某一区域递送基因材料或药物。
为了执行超声治疗,能够使用聚焦超声换能器将超声聚焦到特定的处置体积上。通常将所述换能器安装在能够传输超声的诸如脱气水的介质之内。然后,使用致动器调整超声换能器的位置,并由此调整接受处置的组织区域。
美国专利申请US2004/0122323示出了一种用于对穿过组织的超声能量进行聚焦的系统和用于对组织区域进行成像的成像器。采用相位和幅度校正因子使所述超声能量聚焦。相应的校正因子将补偿各种组织类型之内的各组织类型的声音速度的变化。通过对图像进行分析并从图像识别组织特征来确定这些校正因子。
发明内容
本发明在独立权利要求中提供了一种医学设备、一种方法和一种计算机程序产品。在从属权利要求中给出了实施例。
在通过高强度聚焦超声(HIFU)对乳腺癌进行处置的过程中,超声行进通过若干种介质,诸如,分别以超声波速1451±36m/s和1538m/s±22m/s为特征的皮下脂肪和腺组织。在射束路径之内,这两种组织类型的1cm的厚度的变化将诱发210°的相位像差(aberration),其足以使射束散焦。HIFU处置经常与磁共振成像(MRI)结合,从而对肿瘤进行定位、规划处置以及监测高热流程。
也可以在组织波速和衰减的先验知识的基础上使用由MRI或另一成像模态提供的解剖学图像来处理相位和幅度像差补偿。然而,乳房组织波速(用于计算重聚焦的基本参数)取决于患者和年龄。本发明的实施例可以提供一种通过在医学成像系统内增加超声发射器和接收器而在处置之前或者在处置期间表征组织类型的波速和衰减的设备和/或方法。使用超声发射器和接收器测量超声传播时间和沿若干方向的衰减,而解剖学图像则提供每种组织类型沿这些方向的表征。因而通过结合两种数据获得了波速和衰减值。那些在高热之前或者在高热期间快速获得的组织表征将改善射束重聚焦,从而获得更高的处置效率和准确性。
通过MR引导的HIFU对乳腺癌进行处置最近显示出了进入临床实践的极大的可能性。然而,由于还存在着相当大的对处于超声射束的远场内的胸廓、心脏和肺造成不希望的组织损伤的风险,阻碍了广泛的临床接受。本发明的实施例可以提供一种使用具有横向主射束方向的换能器设计来规避这些风险的方式。另外,射束近场中的较低的能量密度允许对更大的肿瘤进行处置,而不会超过安全范围从而导致皮下组织层的不希望的温度升高。
大孔径换能器设计的焦点的质量对非均匀传播介质的影响更为敏感。由于从声学的角度来讲乳房组织是由腺体组织、脂肪组织和纤维组织构成的复杂结构,这些组织表现出的声传播速度的变化高达6%,使用这样的设计通常得到很差的射束。
本发明的实施例可以提供重新建立射束质量的校正函数。穿越颅骨的HIFU消融系统的设计就遇到了类似的问题,在所述系统中,使用高分辨率CT成像作为校正函数的计算的基础,所述校正函数假定对于所有患者而言类似的组织声学特性。然而,这些方法不能够直接转用到乳房组织上,因为这些方法被优化为表征和校正由通过骨组织的射束传播导致的像差。另外,多种模态的使用(CT和MRI)诱发了与系统可用性、工作流程持续时间和图像配准相关的问题。因此,允许使用以HIFU对乳腺癌进行最佳非侵入性处置的大孔径换能器设计的关键是表征临床相关时间范围内的乳房组织的声学性质的有效的、非侵入性的方法。
MR成像表现出了优越的软组织对比度,并且因此能够以高空间分辨率绘制出异质组织结构。可以采用这一信息与不同乳房组织类型的活体内波速测量的预定义表格的组合来处理重聚焦方法。然而,这一方法存在限制,因为皮下脂肪、腺组织的波速取决于患者和年龄。
有人已经提出了飞行时间超声断层摄影平台,来对乳房组织的波速量化。然而,所得到的测量需要长扫描时间(对于物品而言为每切片9min/s),并且其对比度和空间分辨率局限于超声的现有技术水平。本发明的实施例可以通过使表征超声组织特性的机构与一种或若干种成像模态(除了超声波速测量之外)相结合而使得针对HIFU处置的组织异质性补偿更具患者特异性。
能够使用与患者的解剖学图像结合的超声发射器和接收器来表征超声波速以及超声衰减。可以使用一种成像模态或者结合若干种成像模态来获得解剖学图像,所述成像模态诸如是:
·磁共振成像(MRI)
·辐射摄影
-荧光透视检查
-投影辐射摄影(X射线)
·核医学
-闪烁照相
-SPECT
-正电子发射断层摄影(PET)
·光声成像
·乳房热敏成像法
·断层摄影
-线性断层摄影
-复式断层摄影
-厚层断层摄影法
-矫正角膜学(OPT或OPG)
-计算机断层摄影(CT)
-计算机轴向断层摄影(CAT)
·超声检查法
所述超声发射器和接收器能够测量传播延迟和跨若干组织层的衰减。另外,另一成像模态能够沿超声的传播方向提供下述特征的量化:
·每一组织层的厚度
·组织密度分布
·温度分布
超声发射器和接收器测量与另一成像模态的结合能够在处理最佳重聚焦所必需的患者特异性的基础上改善每种组织类型的超声波速和衰减的评估。
就本发明的实施例而言,可以在患者筛查阶段期间、在规划阶段期间、或者在HIFU处置过程中在没有显著的流程延迟的情况下即时地(onthefly)处理组织表征。超声回波的采集需要数毫秒,其比大部分其他医学成像模态的采集时间少得多。另外,如果将每个系统集成到同一平台中并且如果将采集系统设计为避免相互之间存在的干扰,那么能够在使用其他成像模态的同时执行超声回波的采集。此外,也可以按照任何顺序单独采集数据,并在采集到所有必需的数据之后将其结合到一起,以用于处理步骤。
由于例如可以重新使用最初为了对处置进行规划而采集的MRI图像来处理选定目标区域的像差校正,因而可以不将相关的采集持续时间看作是额外的持续时间。
然而,在同一会话期间同时或者顺次测量解剖学特性(即组织厚度)以及采集超声回波是有益处的:它们将自动处于同一参考系中。如果在两个不同的会话期间单独测量组织厚度和回波采集的测量结果,那么必须在两次测量过程中完全确定患者位置,从而对用于评估飞行时间和解剖学信息的超声回波互配准。
乳腺癌的处置是本发明的实施例的很好的应用,因为能够使超声发射器和接收器具有相互面对的取向,从而容易地沿直线表征乳房组织。另外,利用大孔径换能器的乳腺癌处置对组织像差更为敏感,这使得这一应用愈发重要。
文中使用的“计算机可读存储介质”包含可以存储能够由计算装置的处理器执行的指令的任何有形存储介质。可以将所述计算机可读存储介质称为非暂时计算机可读存储介质。也可以将所述计算机可读存储介质称为有形计算机可读介质。在一些实施例中,计算机可读存储介质还可以能够存储数据,所述数据能够被所述计算装置的处理器访问。计算机可读存储介质的范例包括,但不限于:软盘、磁硬盘驱动器、固态硬盘、闪速储存器、USB优盘、随机存取储存器(RAM)、只读储存器(ROM)、光盘、磁光盘和处理器的寄存器堆。光盘的范例包括紧致磁盘(CD)和数字通用盘(DVD),例如CD-ROM、CD-RW、CD-R、DVD-ROM、DVD-RW或DVD-R盘。计算机可读存储介质一词还指计算机装置能够通过网络或通信链路进行访问的各种类型的记录介质。例如,可以通过调制调解器、互联网或者局域网检索数据。
“计算机存储器”或“存储器”是计算机可读存储介质的范例。计算机存储器是能够由处理器的直接访问的任何储存器。计算机存储器的范例包括,但不限于:RAM储存器、寄存器和寄存器堆。
“计算机储存器”或“储存器”是计算机可读存储介质的范例。计算机储存器可以是任何非易失计算机可读存储介质。计算机储存器的范例包括,但不限于:硬盘驱动器、USB优盘、软盘驱动器、智能卡、DVD、CD-ROM和固态硬盘驱动器。在一些实施例中,所述计算机储存器还可以是计算机存储器,反之亦然。
文中使用的“处理器”涵盖了能够执行程序或机器可执行指令的电子部件。在提及包括“处理器”的计算装置时,应当将其解释为可以包含不只一个处理器或处理核。例如,所述处理器可以是多核处理器。处理器还可以指处于单个计算机系统之内的或者分布于多个计算机系统中的处理器的集合。还应当将计算装置一词解释为可以是指每者均包括一个或多个处理器的计算装置的集合或网络。很多程序的指令都是由多个处理器执行的,这些处理器可以处于同一计算装置之内,甚至可以跨多个计算装置分布。
文中使用的“用户接口”是允许用户或操作人员与计算机或计算机系统交互的接口。也可以将“用户接口”称为“人接口装置”。用户接口可以向操作人员提供信息或数据,并且/或者从操作人员接收信息或数据。用户接口可以使计算机接收到来自操作人员的输入,并且可以将来自计算机的输出提供给用户。换言之,所述用户接口可以允许操作人员控制或操纵计算机,并且所述接口可以允许计算机指示操作人员的控制或操纵的效果。数据或信息在显示器或图形用户界面上的显示是向操作人员提供信息的范例。通过键盘、鼠标、跟踪球、触控板、指示杆、图形输入板、操纵杆、游戏键盘、网络摄影、耳机、变速杆、方向盘、踏板、有线手套、跳舞毯、遥控器和加速度计的数据接收均为能够从操作人员接收信息或数据的用户接口部件的范例。
文中使用的“硬件接口”包含能够使计算机系统的处理器与外部计算装置和/或设备交互和/或对其进行控制的接口。硬件接口可以允许处理器向外部计算装置和/或设备发送控制信号或指令。硬件接口还可以使处理器与外部计算装置和/或设备交换数据。硬件接口的范例包括,但不限于:通用串行总线、IEEE1394端口、并行端口、IEEE1284端口、串行端口、RS-232端口、IEEE-488端口、蓝牙连接、无线局域网连接、TCP/IP连接、以太网连接、控制电压接口、MIDI接口、模拟输入接口和数字输入接口。
“医学图像数据”涵盖了使用医学成像系统采集的一维、二维或三维空间数据。文中使用的“医学成像系统”涵盖了适于采集有关患者的身体结构的信息并构建一维、二维、三维医学图像数据集的设备。能够使用医学图像数据构建用于供医生进行分析的可视化表示。可以使用计算机或控制器执行这一可视化。
文中将磁共振(MR)数据定义为在磁共振成像扫描过程中由磁共振设备的天线记录的原子自旋发射的射频信号的测量结果。文中将磁共振成像图像定义为磁共振成像数据内包含的解剖学数据的重建一维、二维或三维可视化。能够使用处理器或计算机执行这一可视化。磁共振图像是一种示范性的医学图像数据。
文中使用的“超声数据”涵盖了描述穿过对象传输的超声射束的数据。所述超声数据可以包括如下项的任意组合,所述项不限于:发射器的发射时间、所发射的超声的频率、接收器的接收时间、发射与接收之间的延迟、发射幅度、反射幅度以及接收器测得的超声的幅度。例如,如果超声受到骨骼或者其他硬组织结构弹回或震开,那么反射的幅度和延迟可能尤其重要。
文中使用的“超声的速度”包含描述穿过固体或液体介质传播的超声的速度的量。超声的速度一词与下述词语是同义词:超声波速、声音速度、音速、声速、超声速度和超声速。
文中使用的“超声窗口”包含能够发射超声或能量的窗口。通常将薄膜或膜片用作超声窗口。例如,超声窗口可以由BoPET(双轴取向聚对苯二甲酸乙二醇酯)薄膜构成。
在一个方面中,本发明提供了一种医学设备,其包括用于采集超声数据的超声发射器和接收器。所述超声数据描述沿受试体中的至少两条路径的超声的速度。所述超声发射器发射超声,并且所述接收器系统接收超声。例如,能够通过超声发射器发射超声脉冲,并且然后由超声接收器对超声脉冲进行接收。可以由超声接收器确定接收到的超声脉冲的幅度和/或延迟。
所述医学设备还包括用于采集描述受试体中的至少两条路径的医学图像数据的医学成像系统。所述医学设备还包括用于控制所述医学设备的处理器。所述医学设备还包括包含机器可执行指令的存储器。
机器可执行指令的执行使处理器使用所述医学成像系统采集医学图像数据。例如,所述处理器可以向所述医学成像系统发送控制信号,其对所述医学成像系统加以控制以采集医学图像数据。机器可执行指令的执行还使所述处理器使用所述超声发射器和接收器系统采集超声数据。医学图像数据和超声数据的采集顺序是不相关的。例如,可以首先采集医学图像数据,可以首先采集超声数据,或者可以同时采集医学图像数据和超声数据。
机器可执行指令的执行还使所述处理器将医学图像分割成至少两种组织类型。医学图像数据可以描述受试体的解剖特征。在医学图像数据中识别或者分割至少两种不同的组织类型。可以使用标准的分割方法对医学图像数据进行分割。例如,可以识别医学图像数据的解剖学数据中的界标。在另一范例中,将可变形模型拟合至所述医学图像数据。所述机器可执行治疗的执行还使所述处理器确定针对受试体中的至少两条路径的每者在受试体中的至少两个距离。所述至少两个距离对应于超声发射器和接收器系统生成的超声在至少两种组织类型中行进的不同距离。所述机器可执行指令的执行还使所述处理器根据超声数据以及至少两条路径中的每者的至少两个距离来计算超声在至少两种组织类型中的每者中的速度。
这一实施例尤其有利,因为可以使用其准确地确定受试体的多种组织类型中的声音速度或超声速度。例如,如果受试体的一部分包含不只一种组织类型,那么可以预见到在所述不同类型的组织中声音的速度是不同的。通过准确地确定在不同的组织区域中的超声的速度,可以使超声更加准确地聚焦。与宽孔径换能器的情况一样,当各个相控阵列换能器元件的射束路径存在很大不同时,这种做法尤其有价值。沿多条路径发射超声,并且之后使用医学图像数据确定超声在不同组织类型中的每者内行进了多长时间和多大距离。对这些距离的了解与对超声用了多久穿越每条路径的了解相结合允许计算超声的速度。例如,超声数据可以包含超声发射器发射的超声脉冲与接收器接收的超声脉冲之间的延迟。
所述实施例也可以是有利的,因为组织内的超声的速度可能具有诊断价值。
在另一实施例中,所述超声发射器和接收器系统具有第一坐标系。所述医学成像系统具有第二坐标系。所述第一坐标系与所述第二坐标系配准。这一做法是有利的,因为能够使由医学成像系统和超声发射器和接收器系统采集的数据相关。这允许将数据合并到一个公共坐标系中,这是准确确定声音速度所必需的。
在另一实施例中,所述医学设备还包括用于对受试体内的目标区进行声波处理的高强度聚焦超声系统。所述高强度聚焦超声系统包括用于生成超声能量的超声换能器。所述超声换能器包括多个超声换能器元件。所述多个超声换能器元件中的每者的相位是能控制的。例如,可以通过交流电源驱动每一超声换能器元件。提供给多个超声换能器元件中的每者的电源的相位和/或幅度可以是能控制的。
所述指令的执行还使所述处理器根据经分割的医学图像数据和至少两种组织类型中的每者中的超声的速度来计算一组换能器元件相位。例如,如果希望对特定位置进行声波处理,那么可以在经分割的医学图像数据中识别出这一位置。对高强度聚焦超声换能器的位置与相对于两种组织类型的超声的速度之间的关系的了解允许为一组换能器元件中的每者选择适当的相位,从而使所述高强度聚焦超声适当地聚焦。如果不这样做,目标区中的相位可能受到不同组织类型中具有各种不同的超声速度的影响。这一实施例尤其有利,因为其允许将聚焦超声更加准确地聚焦到目标区中。
在另一实施例中,所述医学设备包括箱(tank)。所述箱具有用于接收超声传导液体和受试体的部分的开口。例如,所述箱可以是利用超声传导液体填充的器皿或容器。文中使用的超声传导液体包含适于传输超声能量的任何液体。例如,其可以是凝胶、油或水。利用超声传导液体部分填充所述箱,并且之后将受试体的部分插到箱内。所述部分包含目标区。所述超声发射器和接收器系统适于主要在穿过所述箱的第一平面内发送和接收超声。主要是指将所述超声发射器和接收器系统的发射器和接收器布置为使所述超声发射器生成的超声大部分通过所述第一平面行进。将大部分超声生成为使其通过所述第一平面行进。所述第一平面处于所述开口的下方。所述多个超声换能器元件处于第二平面的下方,其中,所述第二平面处于所述第一平面的下方。
在另一实施例中,使所述超声换能器聚焦到所述第二平面之上。这一实施例是有利的,因为所述超声换能器向上聚焦,并且可以聚焦到箱的顶部或者箱外。
在另一实施例中,使所述超声换能器聚焦到所述第二平面的下方。在这一实施例中,在所述声波处理区域之上采集超声数据。
在另一实施例中,所述至少两条路径穿过受试体。
在另一实施例中,所述指令的执行还使所述处理器使用所述超声发射器和接收器系统采集基线超声数据。沿通过所述超声传导液体的至少两条基线路径来采集所述基线超声数据。所述基线超声数据包括超声幅度测量结果。所述超声数据包括超声幅度测量结果。所述指令的执行还使所述处理器计算至少两种组织类型中的每者的衰减量。根据所述超声数据和所述基线超声数据的超声幅度测量结果并且根据至少两条路径中的每者的至少两个距离来计算衰减值。能够使用这一信息计算每一组织类型中的衰减。
在另一实施例中,所述医学设备还包括用于对受试体和/或超声发射器和接收器系统进行定位的定位系统。所述指令的执行还使所述处理器在经分割的医学图像数据中识别骨组织。所述指令的执行还使所述处理器对受试体和/或超声发射器和接收器系统进行定位,使得所述至少两条路径从被所述骨组织反射开并被接收器系统接收。例如,可以将所述处理器连接至机械致动系统,所述机械致动系统将操纵所述超声发射器和接收器系统的位置。或者,所述处理器可以能够控制患者定位系统,所述患者定位系统能够控制受试体的位置。在这一实施例中,超声进入受试体并受到骨骼反射。例如,在这样的实施例中,可以相对于与受试体接触的超声窗口来定位所述发射器和接收器。这一实施例可以是非常有利的,因为能够在超声不能通过的受试体的附近区域来确定通过组织的超声速度和/或组织内的超声衰减。例如,使超声射束穿过胸廓将是很难的。
在另一实施例中,所述指令的执行还使所述处理器根据至少两条路径中的每者的至少两个距离来计算至少两种组织类型中的每者中的超声衰减。
在另一实施例中,所述指令的执行还使所述处理器根据一组换能器元件相位生成超声控制信号。文中使用的超声控制信号是能够控制所述高强度聚焦超声系统的命令和/或信号。所述超声控制信号使所述高强度聚焦超声系统对目标区进行声波处理。所述指令的执行还使所述处理器向所述高强度聚焦超声系统发送超声控制信号。
在另一实施例中,所述医学成像系统是用于从成像区采集磁共振数据的磁共振成像系统。所述目标区处在所述成像区之内。所述指令的执行还使所述处理器沿从超声换能器到目标区生成的路径来采集磁共振温度测定数据。所述指令的执行还使所述处理器根据经分割的医学图像数据和磁共振温度测定数据来计算一组温度补偿的换能器元件相位。
所述指令的执行还使所述处理器根据所述一组换能器元件相位来生成温度校正的超声控制信号。所述指令的执行还使所述处理器向所述高强度聚焦超声系统发送温度校正的超声控制信号。这一实施例尤其有利,因为可以对通过至少两种组织类型的超声的速度进行温度补偿。身体组织的温度保持相对恒定,然而在使用高强度聚焦超声系统对目标区加热时,可能会发展出区域或热点。另外,目标区之内的声波处理点也可能发热并改变超声的速度。这可能导致轻微的散焦。超声换能器的近场中的发热也可能导致散焦。
在另一实施例中,所述医学成像系统是磁共振成像系统。在这一实施例中,所述医学图像数据可以是一组磁共振图像数据。所述磁共振成像系统可以采集磁共振数据,然后将其重建成磁共振图像数据。
在另一实施例中,所述医学成像系统是计算机断层摄影系统。
在另一实施例中,所述医学成像系统是诊断超声系统。
在另一实施例中,所述至少两条路径是通过受试体的多条或大量路径。然后采用统计学方法准确地确定通过受试体内的至少两种组织类型的超声的幅度和/或速度。
在另一方面中,本发明提供了一种操作医学设备的方法。所述方法还可以提供一种计算机实施的方法。所述医学设备包括用于采集超声数据的超声发射器和接收器系统。所述超声数据描述了沿受试体中的至少两条路径的超声的速度。所述医学设备还包括用于采集描述受试体内的至少两条路径的医学图像数据的医学成像系统。所述方法包括使用医学成像系统采集医学图像数据的步骤。所述方法还包括使用超声发射器和接收器系统采集超声数据的步骤。所述方法还包括将医学图像数据分割成至少两种组织类型的步骤。所述方法还包括确定针对受试体中的至少两条路径中的每者而言超声在受试体中的至少两个距离的步骤。所述至少两个距离对应于超声发射器和接收器系统生成的超声在至少两种组织类型中行进的距离。所述方法还包括根据超声数据以及至少两条路径的每者的至少两个距离来计算至少两种组织类型的每者中的超声的速度的步骤。
所述医学设备还包括用于对受试体内的目标区进行声波处理的高强度聚焦超声系统。所述高强度聚焦超声系统包括用于生成超声能量的超声换能器。所述超声换能器包括多个超声换能器元件。所述多个超声换能器元件中的每者的相位和/或幅度是能控制的。所述指令的执行还使所述处理器根据经分割的医学图像数据和至少两种组织类型中的每者中的超声的速度来计算一组换能器元件相位和/或幅度。
在另一方面中,本发明提供了一种包括机器可执行指令的计算机程序产品,所述指令由控制医学设备的处理器执行。例如,可以将所述计算机程序产品存储到计算机可读存储介质上。所述医学设备包括用于采集超声数据的超声发射器和接收器系统。所述超声数据描述沿受试体内的或者通过受试体的至少两条路径的超声的速度。所述医学设备还包括用于采集描述受试体内的至少两条路径的医学图像数据的医学成像系统。
所述指令的执行使所述处理器使用医学成像系统采集医学图像数据。所述指令的执行还使所述处理器使用超声发射器和接收器系统采集超声数据。所述指令的执行还使所述处理器将医学图像数据分割成至少两种组织类型。所述指令的执行还使所述处理器确定针对受试体中的至少两条路径中的每者而言超声在受试体中的至少两个距离。所述至少两个距离对应于由超声发射器和接收器系统生成的超声在至少两种组织类型中行进的距离。所述指令的执行还使所述处理器根据超声数据和至少两条路径中的每者的至少两个距离来计算至少两种组织类型中的每者内的超声的速度。
附图说明
在下文中,将参考附图仅通过举例的方式描述本发明的优选实施例,在附图中:
图1示出了说明根据本发明的实施例的方法的流程图;
图2示出了根据本发明的实施例的方法的另一实施例;
图3示出了图示说明人类乳房中的不同组织类型的超声波速的差异的曲线图;
图4图示了根据本发明的实施例的医学设备;
图5示出了图4的实施例的高强度聚焦超声系统的更加详细的图示;
图6示出了根据本发明的实施例的超声发射器和接收器系统的范例;
图7示出了根据本发明的另一实施例的超声发射器和接收器系统;
图8示出了根据本发明的实施例的超声换能器元件的布置的等角视图;
图9示出了图8的实施例的超声换能器元件的布置的顶视图;
图10示出了图8的实施例的超声换能器元件的布置的侧视图;
图11图示了与根据本发明的实施例的高强度聚焦超声系统的超声换能器相关的超声发射器和接收器系统的相对位置;
图12图示了与根据本发明的另一实施例的高强度聚焦超声系统的超声换能器相关的超声发射器和接收器系统的相对位置;
图13图示了与根据本发明的另一实施例的高强度聚焦超声系统的超声换能器相关的超声发射器和接收器系统的相对位置;
图14图示了与根据本发明的另一实施例的高强度聚焦超声系统的超声换能器相关的超声发射器和接收器系统的相对位置;以及
图15图示了与根据本发明的另一实施例的高强度聚焦超声系统的超声换能器相关的超声发射器和接收器系统的相对位置。
附图标记列表
具体实施方式
在这些附图中具有类似附图标记的元件要么是等价的元件,要么执行相同的功能。如果元件功能相当,那么前面论述过的元件在后面的附图中将未必再对其进行论述。
图1示出了图示说明根据本发明的实施例的方法的流程图。在步骤100中,使用医学成像系统采集医学图像数据。接下来,在步骤102中,使用超声发射器和接收器采集超声数据。可以按照任何顺序或者彼此同时地执行步骤100和102。接下来,在步骤104中,将医学图像数据分割成至少两种组织类型。在步骤106中,确定在至少两种组织类型的每者中的行进距离。可以通过跟踪由超声发射器和接收器生成的超声行进的至少两条路径的距离来进行这一操作。最后,在步骤108中,计算在所述至少两种组织类型的每者中的超声的速度。根据所述超声数据以及针对至少两条路径在至少两种组织类型的每者中的行进距离来进行所述计算。所述超声数据可以包含所述发射器和接收器对超声的发送和接收之间的延迟的相关信息。
图2示出了根据本发明的实施例的方法的另一实施例的流程图。步骤200-208对应于图1的步骤100-108。在步骤200中,使用所述医学成像系统采集医学图像数据。在步骤202中,使用超声发射器和接收器采集超声数据。针对通过受试体的至少两条路径采集数据。在步骤204中,将医学图像数据分割成至少两种组织类型。通过受试体的路径的数量应当至少与组织类型的数量一样大。在步骤206中,确定在至少两种组织类型的每者中行进的距离。在步骤208,计算至少两种组织类型的每者中的超声的速度。在步骤210中,计算一组换能器元件相位。在步骤212中,使用一组换能器元件相位生成超声控制信号。在步骤214中,将超声控制信号发送至高强度聚焦超声系统。
图3示出了图示说明人乳房中的不同组织类型的超声波速的差异的曲线图。x轴上是温标。所述温标是相对于正常身体内核温度以摄氏度量度的。y轴302以米/秒示出了波速或超声速度。具有附图标记304的线是乳房脂肪中的超声的速度。具有附图标记306的线是腺组织中的超声的速度。在这一曲线图中能够看出,乳房脂肪304和腺组织306的速度是不同的。随着温度的升高,所述差异也变大。
为了说明使超声测量机构与解剖学图像以外的其他图像类型结合的可能性以及使若干种成像模态结合到一起的可能性,我们可以考虑波速与温度的相关性。
对于大多数生物组织而言,在截至50℃的温度范围内声速以大约2.2ms-1℃-1的速度随着温度升高。与之不同的是,对于间质组织中的脂质而言,速度随温度降低,对于人类乳房脂肪而言,以-3.1m/(s℃)的估算斜率降低。
曾有文献报道,在无高热的情况下,乳房腺体和乳房脂肪的速度分别是1538m/s和1451m/s。因而,图3基于文献给出了函数温度升高中的波速,对于腺体组织和脂肪组织而言,其分别假定2.2mms-1℃-1和-3.1ms-1℃-1的恒定斜率。脂肪和腺体组织之间的波速差异伴随着+16.4℃的温度升高以因子2从87m/s升高到174m/s。
由于被设计为补偿乳房组织中的声学像差的相位校正是波速的函数,因此其也是温度的函数。在16.4℃的典型高热过程中,乳房组织之间的波速差异以因子2升高,因而针对温度的调整进行温度分布的量化似乎与不存在高热的情况下对波速的初始量化具有同样的必要性。
根据本发明的方法的实施例可以包括测量乳房组织内的温度分布以进行波速量化调整。也可以利用基于质子共振频率偏移的一项或若干项温度测绘(尤其是对于腺组织而言)和/或基于T1、T2或T2*(star)驰豫变化的热图(尤其是对于脂肪组织而言)来测量空间温度分布。
为了实现对每种组织类型的波速的温度相关性的患者特异性调整,可以通过在声波处理和超声测量之间快速切换或者在不发生干扰的情况下同时执行两操作而与治疗流程同时执行超声测量。能够使用即时采集的热图和超声测量结果实时估算每种组织类型的超声组织特性。例如,可以在将每一组织层细分成具有类似的温度的较小区域的情况下重新使用先前的方程。因此,针对每次出现温度变化获得了每一小区域内的波速和衰减。然而,大量的小区域将需要大量的超声发射器和接收器。为了保持这一系统简单,一种求解这一方程的备选方式可以是对与超声测量结果最佳匹配的每一组织类型的波速的适当温度相关性进行迭代搜索。
图4示出了图示说明根据本发明的实施例的医学设备400的图解。在这一实施例中,医学设备400包括磁共振成像系统402和高强度聚焦超声系统404。
图5示出了高强度聚焦超声系统404的更加详细的图示。在这一实施例中使用磁共振成像系统402。磁共振成像系统402旨在起到例示作用。在其他实施例中也可以使用其他医学成像系统。例如,可以利用计算机断层摄影系统和诊断超声系统来代替磁共振成像系统402。将图4和图5放到一起说明。
磁共振成像系统402包括磁体406。磁体406是圆柱形超导类型磁体并且具有贯穿其的腔膛408。磁体406具有利用超导线圈的液氦冷却低温恒温器。也可能使用永磁体或电阻性磁体。也可能使用不同类型的磁体,例如,也可能使用分裂圆柱形磁体和所谓的开放磁体两者。分裂圆柱形磁体与标准圆柱形磁体类似,只是低温恒温器分裂成了两个部分,从而允许访问磁体的等平面,这样的磁体可以例如与带电粒子束治疗结合使用。开放磁体可以具有两个磁体部分,所述两个部分中的一个处于另一个之上,其间具有大到足够容纳受试体的空间:所述两个部分的布置与亥姆霍兹线圈的布置类似。开放磁体是普遍使用的,因为其对受试体造成限制更少。在所述圆柱形磁体的低温恒温器内有超导线圈的集合。在圆柱形磁体的腔膛408内有成像区418,在成像区中,磁场强,并且具有足以执行磁共振成像的均匀性。
在所述磁体的腔膛内,还有磁场梯度线圈410,其用于采集磁共振数据,并对磁体406内的成像区418之内的磁自旋进行空间编码。将磁场梯度线圈410连接至磁场梯度线圈电源412。磁场梯度线圈410旨在起到代表性作用。典型的磁场梯度线圈410包含三个独立的线圈组,它们用于沿三个正交的空间方向进行空间编码。磁场梯度电源412向磁场梯度线圈供应电流。提供给磁梯度场线圈410的电流是作为时间的函数加以控制的,并且可以使其坡变或者形成脉冲。
射频线圈414与成像区418相邻,其用于操作成像区418之内的磁自旋的取向,并且用于接收来自同样处于成像区418之内的自旋的射频发射。射频天线可以包含多个线圈元件。也可以将射频线圈称为信道、天线和射频天线。射频线圈被连接至射频收发器416。可以使用单独的发射和接收线圈以及单独的发射器和接收器来替代射频线圈414和射频收发器416。应当理解,射频线圈414和射频收发器416只是代表性的。所述射频天线还旨在表示专用发射天线和专用接收天线。类似地,收发器416也可以表示单独的发射器和接收器。
受试体420被安置在磁体406的腔膛408之内的受试体支撑物422上。受试体420的部分处在成像区418之内。高强度聚焦超声系统404具有箱424。利用超声传导液体426填充箱424。由标记428指示液体或液面的顶部。受试体420的部分430通过受试体支撑物422和箱424中的开口423。部分430伸展并被超声传导液体426局部围绕。在这一范例中,两个超声换能器元件432被示为处在液体426之内。受试体420的部分430为乳房。存在两个组织区域:脂肪组织434区域和腺组织436区域。虚线438描绘了从超声换能器元件432通过液体426、通过脂肪组织434并且通过腺组织436到位于目标体积440之内的声波处理体积442的超声路径438。
在所述超声换能器元件之上是超声发射器444和超声接收器446。超声射束448从超声发射器444行进至超声接收器446。存在由虚线指示的第一平面450和第二平面452。超声发射器444和超声接收器446主要处在第一平面450之内。超声主要穿过所述第一平面行进。还存在第二平面452。第二平面452处于第一平面450的下方。超声换能器元件432处于第一平面452的下方,但是声波处理体积442处于第二平面452的上方。这一实施例是有利的,因为其允许直接测量紧密靠近声波处理体积的乳房组织之内的超声的速度。实际上,在一些实施例中,声波处理体积442主要位于第一平面450中。平面450当然是二维空间。超声射束448和声波处理体积442是三维体积。如果某一平面截穿所述体积的一部分,那么本文中使用的所述体积主要处于该平面内。
高强度聚焦超声系统404、磁场梯度线圈电源412和收发器416被示为连接至计算机系统454的硬件接口456。硬件接口456连接至处理器458。处理器458还连接至用户接口460、计算机储存器462和计算机存储器464。
计算机储存器462被示为包含由磁共振成像系统402采集的磁共振数据466。存在从磁共振数据466重建的磁共振图像数据468。磁共振图像数据468是医学图像数据的范例。磁共振图像数据468的图像分割470也存储在计算机储存器462之内。图像分割470包含描述至少两种组织类型434、436的位置的信息。在所述计算机储存器之内有热图472。所述热图是从磁共振数据466重建的。所述磁共振数据包括磁共振温度测定数据。所述计算机储存器还被示为包含距离数据474。所述距离数据描述超声射束448穿过第一和第二组织类型434、436行进的距离。第一组织类型是脂肪组织434,而第二组织类型是腺组织436。所述计算机储存器被示为包含超声数据476。所述超声数据是由超声接收器446采集的数据。
所述计算机储存器被示为包含从超声数据476提取的脉冲延迟数据478和幅度数据480。脉冲延迟数据478描述从发射器444行进至接收器446的超声射束448中的延迟。幅度延迟数据480描述从发射器444到接收器446的幅度变化。所述计算机储存器还被示为包含超声控制信号484。所述超声控制信号包含用于操作高强度聚焦超声系统404的指令。所述计算机储存器还被示为包含使用脉冲延迟数据478和距离数据474计算的超声速度486。距离数据474和幅度延迟数据480可以用于计算两种组织类型434、436中的超声衰减。
计算机存储器464被示为包含控制模块490。控制模块490包含用于控制医学设备400的操作和功用的计算机可执行代码。计算机存储器464还被示为包含脉冲序列492,其包含由磁共振成像系统402执行以采集磁共振数据466的操作。还有磁共振成像系统控制模块494,其使用脉冲序列492来控制磁共振成像系统402。计算机存储器464还被示为包含图像模块496。图像模块496包含用于将磁共振数据466重建成磁共振图像数据468和/或用于从磁共振图像数据468执行和生成图像分割470的计算机可执行代码。图像模块496还可以包含用于从磁共振数据466生成热图472的计算机可执行代码。
计算机存储器464还被示为包含高强度聚焦超声控制模块498。所述高强度聚焦超声控制模块498包含用于控制高强度聚焦超声系统404的操作和功能的计算机可执行代码。所述高强度聚焦超声控制模块498适于生成超声控制信号484。计算机存储器464还被示为包含图像分析模块500。所述图像分析模块能够使用图像分割470和磁共振图像数据468来计算距离数据474。
计算机存储器464还被示为包含速度计算模块502和衰减计算模块504。所述速度计算模块包含适于从距离数据474和脉冲延迟数据478计算超声速度486的计算机可执行代码。衰减计算模块504包含用于使用幅度数据480和距离数据474来计算超声衰减的计算机可执行代码。
图6示出了根据本发明的实施例的超声发射器和接收器系统600的范例。在该图中示出了四个发射器602。接收器604处于每一发射器602的对面。利用606指示的虚线跟踪从发射器602到接收器604的路径。路径606穿过超声传导液体608、第一组织类型610和第二组织类型612。为了确定第一和第二组织类型610、612中的速度,医学成像系统可以采集处于与路径606所在的相同的平面内的医学图像数据。能够针对每条路径计算在超声传导液体606、第一组织类型610和第二组织类型612中行进的距离。通过确定从发射器602到接收器604的传输之间的延迟能够计算出超声传导液体和每种组织类型610、612中的波速。
在任何MR-HIFU处置之前,患者经受各种类型的医学检查,以确立诊断。这一筛查阶段对于相对于其他处置类型预测HIFU治疗的成功而言是很关键的。能够使用超声组织特性的表征手段与筛查成像模态(除了超声波速测量以外)的结合使用来改进处置选项。
对于乳腺癌的HIFU处置而言,了解组织的波速和形状是很重要的,因为这样能够仅在有限的值范围内校正所产生的相位像差。因而,以通过先前描述的方法获得的组织模型为基础的声场模拟能够为通过HIFU对患者进行处置的可能性做出答复。
对于这样的筛查过程而言,能够将超声测量机构集成到独立于HIFU系统的合并了一种和若干种成像模态(除了超声波速测量以外)的系统内,以便利用与先前描述的方法相同的方法来确定超声波速和每种组织类型的衰减。
例如,能够将图6中描述的环结构插入到通常用于检测乳腺癌的MRI乳房造影术轨道台内。MRI乳房造影术轨道台相当于图4和图5所示的受试体支撑物442。在这种情况下,能够将超声发射器和接收器放置为在MR线圈接收器环之间传播声波。
一旦在患者筛查步骤中对每种组织类型的超声波速和衰减做出了量化并且该患者适于进行HIFU处置,就能够将所得到的组织的表征再次用于处置过程。由于每种组织的超声波速和衰减在筛查日和治疗日在很大的概率上是相同的,因而仅需对经分割组织层进行配准或者对组织层进行新的分割,以处理相位和强度校正,使之适用于使射束重新聚焦。当然,如有必要,在想要确保超声波速和衰减测量的可靠性的情况下,可以在治疗日重新采集所述测量结果。然而,如果那些测量很耗时,而且在筛查步骤中已经采集了,那么可以在治疗日避免对超声波速和衰减进行新的量化,以降低总的处置持续时间。
磁共振高强度聚焦超声(MR-HIFU)平台能够包括超声发射器/接收器系统,其负责对由MRI分割的组织体积的波速进行量化。在图6所示的范例中,这一系统包括放置在附接至刚性环上的4个接收器的前面的4个发射器构成,从而沿4个方向表征乳房组织。
超声脉冲的发射和接收提供了沿每一方向i的超声行进时间Ti和压力幅度Pi。类似地,也能够在先前的系统校准步骤中根据乳房与换能器之间的耦合介质仅使用水或者仅使用油来完成TRef i和PRef i的量化。沿一条线的测量持续时间是非常快的,因为经过15cm的距离的超声行进时间大约为0.1ms。因此,能够在不会给治疗流程持续时间带来显著的延迟的情况下出于求平均的目的重复测量。这样的求平均提高了延迟和压力幅度的测量的准确性。
MR图像的分割算法提供了对脂肪组织F、腺组织G和水W的描绘。跨每种介质表面从发射器延伸至接收器的线的交叉提供了长度LW i、LF i、LG i的量化。然而,对LW i的量化是非强制性的。
行进时间根据如下线性方程通过水、脂肪和腺体介质的波速CF、CG、CW和衰减AF、AG、AW相互联系。
两方程的相减消去了长度LW i,其使得乳房组织波速的量化更加准确。通过MRI对长度LW i的测量结果由于换能器附近的与磁化率场的急剧变化相关的图像失真的原因而变得不准确。TRef i和Ti之间的传播时间的差异对应于两个类似的波形的比较,可以通过互相关非常精确度估计所述波形。
其结果是通过前述方程的反演而获得了波速CF和CG:
根据如下线性方程通过水、脂肪和腺体介质的衰减αF、αG和αW,以及发射的压力幅度P0使类似的各个压力幅度相互联系:
两方程相减消去了长度LW i和压力幅度P0,其使得组织衰减的量化更加准确。
其结果是通过前述方程获得了波速αF和αG:
能够使这种方法一般化以表征2种以上的组织类型,直到组织类型的数量小于等于测量方向的数量为止。然而,在方向的数量大于组织类型的数量时,例如,这一范例具有2种组织类型和4个方向,此时将获得更高的组织特征测量准确性。
图7示出了根据本发明的实施例的超声发射器和接收器系统700。在这一范例中,有单个发射器702和多个接收器704、超声传导液体区域706、第一组织类型708、第二组织类型720和骨骼组织712。超声从发射器702通过超声传导液体、第一组织类型708、第二组织类型710沿路径714行进。然后,所述超声受到骨骼712的反射。虚线716示出了回到接收器704的超声的路径。能够使用在路径714和716的平面内采集的医学图像数据来计算在所述超声传导液体、第一组织类型708和第二组织类型710内行进的距离。可以将其与超声的发射702和不同接收器704对其的接收之间的延迟结合使用,以计算超声在所述组织类型中的波速。
可以考虑备选方法,其使用超声在骨骼上的反射来表征其他身体区域内的皮下脂肪和肌肉。在骨转移以及肝肾恶性肿瘤的处置过程中将频繁地遇到骨骼结构。
图8、图9和图10示出了超声换能器元件802的布置的等角视图800、顶视图900和侧视图1000。小圈中的每者是单个超声换能器元件。虚线804示出了通往声波处理体积808的超声路径。形状806通过三个不同的视图800、900和1000示出了相同的组织结构以供参考。
图11-15示出了与高强度聚焦超声系统的超声换能器相关的超声发射器和接收器系统的不同构造。
在图11中,区域1100表示受试体。其示出了两个超声换能器元件。利用1104标记的箭头和周围的锥形区域示出了治疗超声1104的路径,在该路径中,将所述超声聚焦到受试体1100的声波处理体积1106之内。围绕声波处理体积1106存在受热区域1108。在这一实施例中,使超声换能器元件1102向上聚焦。在位于超声换能器元件1102之上的平面中,存在用于沿穿过受试体的路径1114传输超声脉冲的发射器1110和接收器1112。在这一实施例中,在受试体1100与发射器110之间存在空间,在受试体1100与接收器1112之间也存在空间1116。
图12示出了非常类似的实施例,只是在这一实施例中不存在任何空间。发射器1110与受试体1100直接接触。接收器1112也与受试体1100直接接触。在前面的图11所示的实施例中,可以使用凝胶垫或者其他类似的超声传导介质来填充空间1116。
在图13所示的实施例中,超声换能器元件1102仍然具有与它们在图11和12中所具有的相同的关系。然而,在这一实施例中,没有单独的发射器1110。超声换能器元件1102起着图11和12所示的发射器1110的作用。有两个接收器1112,其现在处于沿治疗超声1104的路径的位置上。
在图14中,示出了具有不同的构造的超声换能器元件。在这一实施例中,超声换能器元件1102被示为处于受试体1100的侧面。治疗超声1104的路径再次聚焦到受试体1100之内的声波处理体积1106。受热区域1108围绕声波处理体积1106。治疗超声1104离开区域1108并传播至接收器1112。与图13所示的实施例一样没有发射器。超声换能器元件1102起着发射器1110的作用。
图15所示的实施例与图14所示的类似。在这一实施例中,有两个超声换能器元件1102。这些元件处于受试体1100的相邻侧面。通过两超声换能器元件1102将治疗超声1104的路径聚焦到声波处理体积1106。而且,在这一实施例中,由超声换能器元件1102替代发射器。接收器1112处于每一超声换能器元件1102上。通过受试体1114的路径等同于治疗超声1104的路径。
尽管在附图和上述说明中详细图示说明和描述了本发明,但是应当将这样的图示说明和描述看作是说明性或者示范性的,而非限制性的;本发明不限于所公开的实施例。
通过研究附图、说明书和权利要求书,本领域技术人员能够在对所要求保护的本发明的实践中理解并实施针对所公开的实施例的其他变型。在权利要求中,“包括”一词不排除其他元件或步骤,单数冠词“一”或“一个”不排除复数。单个处理器或其他单元可以完成权利要求中列举的若干项的功能。在互不相同的从属权利要求中记载某些措施不表示不能有利地采用这些措施的组合。可以将计算机程序存储/分布在适当的介质中,例如,所述介质可以是光存储介质或者与其他硬件一起提供的或者作为其他硬件的部分的固体介质,但是,也可以使所述计算机程序通过其他形式分布,例如,通过互联网或者其他有线或无线电信系统。不应将权利要求中的附图标记推断为有限制范围的作用。
Claims (13)
1.一种医学设备(400),包括:
-超声发射器(444、602)和接收器(446、604)系统(600),其用于采集超声数据(476),其中,所述超声数据描述沿受试体(430、1100)中的至少两条路径(606、1114)的超声的速度(304、306);
-医学成像系统(402),其用于采集描述所述受试体中的所述至少两条路径的医学图像数据(468);
-存储器(464),其包含机器可执行指令(490、492、494、496、498、500、502、504);以及
-处理器(458),其用于控制所述医学设备并且被配置为:
-使用所述医学成像系统采集(100、200)所述医学图像数据;
-使用所述超声发射器和接收器系统采集(102、202)所述超声数据;
-对所述医学图像数据进行分割(104、204),识别至少两种组织类型(434、436、610、612、708、710);
-确定(106、206)在所述受试体中的针对所述受试体中的所述至少两条路径中的每者的至少两个距离,其中,所述至少两个距离等于由所述超声发射器和接收器系统生成的超声在所述至少两种组织类型中行进的距离;
-根据所述超声数据并且根据针对所述至少两条路径中的每者的所述至少两个距离来计算(108、208)在所述至少两种组织类型中的每者中的超声的速度,其中,所述医学设备还包括用于对所述受试体之内的目标区(440)进行声波处理的高强度聚焦超声系统(404),其中,所述高强度聚焦超声系统包括用于生成超声能量的超声换能器,其中,所述超声换能器包括多个超声换能器元件(432、802、1102),其中,所述多个超声换能器元件中的每者的相位是能控制的;其中,所述处理器还被配置为根据经分割的医学图像数据(470)和所述至少两种组织类型中的每者中的所述超声的速度来计算(210)一组换能器元件相位(482)。
2.根据权利要求1所述的医学设备,其中,所述医学设备包括箱(424),其中,所述箱具有用于接收超声传导液体(426)和所述受试体的部分(430)的开口(423),其中,所述部分包含所述目标区,其中,所述超声发射器和接收器系统适于主要在穿过所述箱的第一平面(450)之内发射和接收超声,并且其中,所述第一平面处于所述开口的下方,其中,所述多个超声换能器元件处于第二平面(452)的下方,其中,所述第二平面处于所述第一平面的下方。
3.根据权利要求2所述的医学设备,其中,所述超声换能器聚焦在所述第二平面的上方。
4.根据权利要求2所述的医学设备,其中,所述超声换能器聚焦在所述第二平面的下方。
5.根据权利要求2、3或4中的任一项所述的医学设备,其中,所述至少两条路径穿过所述受试体。
6.根据权利要求5所述的医学设备,其中,所述处理器被配置为:
-使用所述超声发射器和接收器系统采集基线超声数据,其中,沿穿过所述超声传导液体的至少两条基线路径来采集所述基线超声数据,其中,所述基线超声数据包括超声幅度测量结果,其中,所述超声数据包括超声幅度测量结果;并且
-计算针对所述至少两种组织类型中的每者的衰减值,其中,所述衰减值是根据所述超声数据和所述基线超声数据的所述超声幅度测量结果并且根据针对所述至少两条路径中的每者的所述至少两个距离来计算的。
7.根据权利要求1到4中的任一项所述的医学设备,其中,所述医学设备还包括用于对所述受试体和/或所述超声发射器和接收器系统进行定位的定位系统,其中,所述处理器被配置为:
-在所述经分割的医学图像数据中识别骨组织(712);并且
-对所述受试体和/或所述超声发射器和接收器系统进行定位,使得所述至少两条路径从所述骨组织反射开并被所述接收器系统接收。
8.根据权利要求1到4中的任一项所述的医学设备,其中,所述处理器被配置为根据针对所述至少两条路径中的每者的所述至少两个距离来计算在所述至少两种组织类型中的每者中的超声衰减。
9.根据权利要求1到4中的任一项所述的医学设备,其中,所述处理器被配置为:
-根据所述一组换能器元件相位生成超声控制信号(484),其中,所述超声控制信号使所述高强度聚焦超声系统对所述目标区进行声波处理;并且
-将所述超声控制信号发送至所述高强度聚焦超声系统。
10.根据权利要求9所述的医学设备,其中,所述医学成像系统是用于从成像区采集磁共振数据的磁共振成像系统,其中,所述目标区处在所述成像区之内,其中,所述处理器被配置为:
-沿从所述超声换能器到所述目标区生成的路径来采集磁共振温度测定数据;
-根据所述经分割的医学图像数据和所述磁共振温度测定数据来计算一组温度补偿的换能器元件相位;并且
-根据所述一组换能器元件相位来生成温度校正的超声控制信号;并且
-将所述温度校正的超声控制信号发送至所述高强度聚焦超声系统。
11.根据权利要求1到4中的任一项所述的医学设备,其中,所述医学成像系统是如下中的任一项:磁共振成像系统(406)、计算机断层摄影系统和诊断超声系统。
12.一种操作医学设备的方法,其中,所述医学设备包括用于采集超声数据(476)的超声发射器(444、602)和接收器(446、604)系统(600),其中,所述超声数据描述沿受试体(430、1100)中的至少两条路径(606、1114)的超声的速度(304、306);其中,所述医学设备还包括用于采集描述所述受试体中的所述至少两条路径的医学图像数据(468)的医学成像系统(402),其中,所述方法包括如下步骤:
-使用所述医学成像系统采集(100、200)所述医学图像数据;
-使用所述超声发射器和接收器系统采集(102、202)所述超声数据;
-将所述医学图像数据分割(104、204)成至少两种组织类型(434、436、610、612、708、710);
-确定(106、206)超声在所述受试体中的针对所述受试体中的所述至少两条路径中的每者的至少两个距离,其中,所述至少两个距离对应于由所述超声发射器和接收器系统生成的超声在所述至少两种组织类型中行进的距离;并且
-根据所述超声数据并且根据针对所述至少两条路径中的每者的所述至少两个距离来计算(108、208)在所述至少两种组织类型中的每者中的超声的速度,其中,所述医学设备还包括用于对所述受试体之内的目标区(440)进行声波处理的高强度聚焦超声系统(404),其中,所述高强度聚焦超声系统包括用于生成超声能量的超声换能器,其中,所述超声换能器包括多个超声换能器元件(432、802、1102),其中,所述多个超声换能器元件中的每者的相位是能控制的;其中,所述方法还包括根据经分割的医学图像数据和所述至少两种组织类型中的每者中的所述超声的速度来计算(210)一组换能器元件相位。
13.一种操作医学设备的装置,其中,所述医学设备包括用于采集超声数据(476)的超声发射器(444、602)和接收器(446、604)系统(600),其中,所述超声数据描述沿受试体(430、1100)中的至少两条路径(606、1114)的超声的速度(304、306),其中,所述医学设备还包括用于采集描述所述受试体中的所述至少两条路径的医学图像数据(468)的医学成像系统(402),所述装置包括:
-用于使用所述医学成像系统采集(100、200)所述医学图像数据的模块;
-用于使用所述超声发射器和接收器系统采集(102、202)所述超声数据的模块;
-用于将所述医学图像数据分割(104、204)成至少两种组织类型(434、436、610、612、708、710)的模块;
-用于确定(106、206)超声在所述受试体中的针对所述受试体中的所述至少两条路径中的每者的至少两个距离的模块,其中,所述至少两个距离对应于由所述超声发射器和接收器系统生成的超声在所述至少两种组织类型中行进的距离;并且
-用于根据所述超声数据并且根据针对所述至少两条路径中的每者的所述至少两个距离来计算(108、208)在所述至少两种组织类型中的每者中的所述超声的速度的模块,并且
-用于控制包括用于生成超声能量的超声换能器的高强度聚焦超声系统的模块,所述超声换能器包括多个超声换能器元件(432、802、1102),其中,所述多个超声换能器元件中的每者的相位是受到控制的,并且根据经分割的医学图像数据和所述至少两种组织类型中的每者中的所述超声的速度来计算(210)一组换能器元件相位。
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
EP11156819.2 | 2011-03-03 | ||
EP11156819A EP2494925A1 (en) | 2011-03-03 | 2011-03-03 | Calculating the speed of ultrasound in at least two tissue types |
PCT/IB2012/050886 WO2012117328A1 (en) | 2011-03-03 | 2012-02-27 | Calculating the speed of ultrasound in at least two tissue types |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103415259A CN103415259A (zh) | 2013-11-27 |
CN103415259B true CN103415259B (zh) | 2016-01-20 |
Family
ID=44305084
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201280011510.8A Expired - Fee Related CN103415259B (zh) | 2011-03-03 | 2012-02-27 | 计算至少两种组织类型中超声的速度 |
Country Status (5)
Country | Link |
---|---|
US (1) | US9420989B2 (zh) |
EP (2) | EP2494925A1 (zh) |
JP (1) | JP6092788B2 (zh) |
CN (1) | CN103415259B (zh) |
WO (1) | WO2012117328A1 (zh) |
Families Citing this family (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2722012B1 (en) * | 2012-10-18 | 2015-06-03 | Storz Medical AG | Device for shock wave treatment of the human brain |
JP6040106B2 (ja) * | 2013-06-25 | 2016-12-07 | 株式会社日立製作所 | 超音波診断装置、方法、及び超音波プローブ |
US11191519B2 (en) * | 2014-08-05 | 2021-12-07 | HABICO, Inc. | Device, system, and method for hemispheric breast imaging |
WO2016115270A1 (en) * | 2015-01-14 | 2016-07-21 | Qi2 Elements, Llc | Automatic transducer operating parameter selection |
WO2017040866A1 (en) * | 2015-09-01 | 2017-03-09 | Delphinus Medical Technologies, Inc. | Tissue imaging and analysis using ultrasound waveform tomography |
US20190046158A1 (en) * | 2015-09-16 | 2019-02-14 | Koninklijke Philips N.V. | Ultrasound apparatus and method for medical examination of a subject |
CN109640830B (zh) * | 2016-07-14 | 2021-10-19 | 医视特有限公司 | 基于先例的超声聚焦 |
US10765892B1 (en) * | 2017-06-05 | 2020-09-08 | Insightec, Ltd. | Systems and methods for optimizing transcranial ultrasound focusing |
WO2019219549A1 (en) * | 2018-05-15 | 2019-11-21 | Koninklijke Philips N.V. | Synthetic transmit focusing ultrasound system with speed of sound aberration correction |
WO2020036876A1 (en) * | 2018-08-12 | 2020-02-20 | The Trustees Of Columbia University In The City Of New York | System, method, and computer-accessible medium for non-invasive temperature estimation |
US20220233890A1 (en) * | 2019-05-31 | 2022-07-28 | Sunnybrook Research Institute | Systems and methods for reducing thermal skull-induced aberrations during transcranial ultrasound therapeutic procedures |
CN114072067A (zh) * | 2019-09-04 | 2022-02-18 | 深圳迈瑞生物医疗电子股份有限公司 | 确定声速的方法及超声装置 |
CN110680398A (zh) * | 2019-10-11 | 2020-01-14 | 深圳开立生物医疗科技股份有限公司 | 一种超声回波成像方法、系统及超声设备 |
CN117065239B (zh) * | 2023-10-14 | 2023-12-12 | 北京小超科技有限公司 | 一种超声聚焦补偿的方法和系统 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5522869A (en) * | 1994-05-17 | 1996-06-04 | Burdette; Everette C. | Ultrasound device for use in a thermotherapy apparatus |
CN1297527A (zh) * | 1998-03-03 | 2001-05-30 | 阳光医学有限公司 | 骨中音速的确定 |
CN1732031A (zh) * | 2002-12-23 | 2006-02-08 | 因赛泰克-特克斯索尼克斯公司 | 超声波治疗中的组织偏差校正 |
CN1968656A (zh) * | 2004-06-11 | 2007-05-23 | 株式会社日立医药 | 超声波治疗装置 |
CN101273896A (zh) * | 2007-03-30 | 2008-10-01 | 通用电气公司 | 使用便携平板检测器的双能射线照相的图像采集和处理链 |
Family Cites Families (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP3860227B2 (ja) | 1993-03-10 | 2006-12-20 | 株式会社東芝 | Mriガイド下で用いる超音波治療装置 |
DE69840444D1 (de) | 1997-05-23 | 2009-02-26 | Prorhythm Inc | Wegwerfbarer fokussierender ultraschallapplikator hoher intensität |
US7112173B1 (en) * | 1998-03-03 | 2006-09-26 | Sunlight Medical Ltd. | Determination of acoustic velocity in bone |
US7379769B2 (en) * | 2003-09-30 | 2008-05-27 | Sunnybrook Health Sciences Center | Hybrid imaging method to monitor medical device delivery and patient support for use in the method |
JP2008509777A (ja) * | 2004-08-17 | 2008-04-03 | テクニオン リサーチ アンド ディベロップメント ファウンデーション リミテッド | 超音波を用いた画像誘導による組織損傷の処置 |
KR101328103B1 (ko) * | 2004-10-06 | 2013-11-13 | 가이디드 테라피 시스템스, 엘.엘.씨. | 비침습적 미용 증진 방법 및 시스템 |
EP1987774A1 (de) * | 2007-05-03 | 2008-11-05 | BrainLAB AG | Messung der Sonographie-Schallgeschwindigkeit mittels Markereinrichtung |
JP5191183B2 (ja) | 2007-07-31 | 2013-04-24 | 株式会社日立メディコ | 超音波診断装置 |
CN102170938B (zh) * | 2008-09-30 | 2015-01-14 | 皇家飞利浦电子股份有限公司 | 用于超声治疗处置的系统和方法 |
-
2011
- 2011-03-03 EP EP11156819A patent/EP2494925A1/en not_active Withdrawn
-
2012
- 2012-02-27 US US14/002,270 patent/US9420989B2/en active Active
- 2012-02-27 CN CN201280011510.8A patent/CN103415259B/zh not_active Expired - Fee Related
- 2012-02-27 EP EP12710790.2A patent/EP2680756B1/en not_active Not-in-force
- 2012-02-27 WO PCT/IB2012/050886 patent/WO2012117328A1/en active Application Filing
- 2012-02-27 JP JP2013555968A patent/JP6092788B2/ja not_active Expired - Fee Related
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5522869A (en) * | 1994-05-17 | 1996-06-04 | Burdette; Everette C. | Ultrasound device for use in a thermotherapy apparatus |
CN1297527A (zh) * | 1998-03-03 | 2001-05-30 | 阳光医学有限公司 | 骨中音速的确定 |
CN1732031A (zh) * | 2002-12-23 | 2006-02-08 | 因赛泰克-特克斯索尼克斯公司 | 超声波治疗中的组织偏差校正 |
CN1968656A (zh) * | 2004-06-11 | 2007-05-23 | 株式会社日立医药 | 超声波治疗装置 |
CN101273896A (zh) * | 2007-03-30 | 2008-10-01 | 通用电气公司 | 使用便携平板检测器的双能射线照相的图像采集和处理链 |
Also Published As
Publication number | Publication date |
---|---|
WO2012117328A1 (en) | 2012-09-07 |
EP2680756A1 (en) | 2014-01-08 |
US9420989B2 (en) | 2016-08-23 |
JP2014517704A (ja) | 2014-07-24 |
EP2680756B1 (en) | 2018-04-11 |
CN103415259A (zh) | 2013-11-27 |
US20130338485A1 (en) | 2013-12-19 |
JP6092788B2 (ja) | 2017-03-08 |
EP2494925A1 (en) | 2012-09-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103415259B (zh) | 计算至少两种组织类型中超声的速度 | |
Marsac et al. | Ex vivo optimisation of a heterogeneous speed of sound model of the human skull for non-invasive transcranial focused ultrasound at 1 MHz | |
Sastry et al. | Applications of ultrasound in the resection of brain tumors | |
Jensen et al. | Spatiotemporal monitoring of high-intensity focused ultrasound therapy with passive acoustic mapping | |
CN108135565A (zh) | 用于配准使用各种成像模态获得的图像并验证图像配准的系统和方法 | |
Miller et al. | Ultrashort echo‐time MRI versus CT for skull aberration correction in MR‐guided transcranial focused ultrasound: In vitro comparison on human calvaria | |
EP2519324B1 (en) | Therapeutic apparatus | |
CN102448547B (zh) | Mr成像引导的治疗 | |
Arthur et al. | Temperature dependence of ultrasonic backscattered energy in motion compensated images | |
Fontanarosa et al. | A CT based correction method for speed of sound aberration for ultrasound based image guided radiotherapy | |
Mougenot et al. | High intensity focused ultrasound with large aperture transducers: a MRI based focal point correction for tissue heterogeneity | |
Zhang et al. | Real‐time sonoelastography of hepatic thermal lesions in a swine model | |
WO2020058757A1 (en) | Ultrasound focusing utilizing a 3d-printed skull replica | |
Kaye et al. | Application of Zernike polynomials towards accelerated adaptive focusing of transcranial high intensity focused ultrasound | |
Holbrook et al. | In vivo MR acoustic radiation force imaging in the porcine liver | |
Canney et al. | A multi-element interstitial ultrasound applicator for the thermal therapy of brain tumors | |
Daoud et al. | Needle detection in curvilinear ultrasound images based on the reflection pattern of circular ultrasound waves | |
Angla et al. | Transcranial ultrasound simulations: A review | |
Teixeira et al. | Feasibility of non-invasive temperature estimation by the assessment of the average gray-level content of B-mode images | |
Liu et al. | Focal beam distortion and treatment planning for transrib focused ultrasound thermal therapy: A feasibility study using a two‐dimensional ultrasound phased array | |
Antoniou et al. | Tumor phantom model for MRI‐guided focused ultrasound ablation studies | |
Yang et al. | Comparison of three dimensional strain volume reconstructions using SOUPR and wobbler based acquisitions: A phantom study | |
Han | Development of a Harmonic Motion Imaging guided Focused Ultrasound system for breast tumor characterization and treatment monitoring | |
Audigier et al. | Physics-based simulation to enable ultrasound monitoring of hifu ablation: An mri validation | |
Hertzberg et al. | Rapid method for assessing relative tissue stiffness using MR acoustic radiation force imaging |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20160120 |