CN112566694B - 包括用于使空化泡成像的装置的超声处理设备 - Google Patents

包括用于使空化泡成像的装置的超声处理设备 Download PDF

Info

Publication number
CN112566694B
CN112566694B CN201980035031.1A CN201980035031A CN112566694B CN 112566694 B CN112566694 B CN 112566694B CN 201980035031 A CN201980035031 A CN 201980035031A CN 112566694 B CN112566694 B CN 112566694B
Authority
CN
China
Prior art keywords
ultrasound
transducer
ultrasound transducer
imaging
therapeutic
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
CN201980035031.1A
Other languages
English (en)
Other versions
CN112566694A (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.)
Cardiawave
Centre National de la Recherche Scientifique CNRS
Institut National de la Sante et de la Recherche Medicale INSERM
Ecole Superieure de Physique et Chimie Industrielles de Ville Paris
Original Assignee
Cardiawave
Centre National de la Recherche Scientifique CNRS
Institut National de la Sante et de la Recherche Medicale INSERM
Ecole Superieure de Physique et Chimie Industrielles de Ville Paris
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 Cardiawave, Centre National de la Recherche Scientifique CNRS, Institut National de la Sante et de la Recherche Medicale INSERM, Ecole Superieure de Physique et Chimie Industrielles de Ville Paris filed Critical Cardiawave
Publication of CN112566694A publication Critical patent/CN112566694A/zh
Application granted granted Critical
Publication of CN112566694B publication Critical patent/CN112566694B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B17/00—Surgical instruments, devices or methods, e.g. tourniquets
    • A61B17/22—Implements for squeezing-off ulcers or the like on the inside of inner organs of the body; Implements for scraping-out cavities of body organs, e.g. bones; Calculus removers; Calculus smashing apparatus; Apparatus for removing obstructions in blood vessels, not otherwise provided for
    • A61B17/225—Implements for squeezing-off ulcers or the like on the inside of inner organs of the body; Implements for scraping-out cavities of body organs, e.g. bones; Calculus removers; Calculus smashing apparatus; Apparatus for removing obstructions in blood vessels, not otherwise provided for for extracorporeal shock wave lithotripsy [ESWL], e.g. by using ultrasonic waves
    • A61B17/2256—Implements for squeezing-off ulcers or the like on the inside of inner organs of the body; Implements for scraping-out cavities of body organs, e.g. bones; Calculus removers; Calculus smashing apparatus; Apparatus for removing obstructions in blood vessels, not otherwise provided for for extracorporeal shock wave lithotripsy [ESWL], e.g. by using ultrasonic waves with means for locating or checking the concrement, e.g. X-ray apparatus, imaging means
    • A61B17/2258—Implements for squeezing-off ulcers or the like on the inside of inner organs of the body; Implements for scraping-out cavities of body organs, e.g. bones; Calculus removers; Calculus smashing apparatus; Apparatus for removing obstructions in blood vessels, not otherwise provided for for extracorporeal shock wave lithotripsy [ESWL], e.g. by using ultrasonic waves with means for locating or checking the concrement, e.g. X-ray apparatus, imaging means integrated in a central portion of the shock wave apparatus
    • 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
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61N—ELECTROTHERAPY; MAGNETOTHERAPY; RADIATION THERAPY; ULTRASOUND THERAPY
    • A61N7/00—Ultrasound therapy
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B17/00—Surgical instruments, devices or methods, e.g. tourniquets
    • A61B17/22—Implements for squeezing-off ulcers or the like on the inside of inner organs of the body; Implements for scraping-out cavities of body organs, e.g. bones; Calculus removers; Calculus smashing apparatus; Apparatus for removing obstructions in blood vessels, not otherwise provided for
    • A61B17/22004—Implements for squeezing-off ulcers or the like on the inside of inner organs of the body; Implements for scraping-out cavities of body organs, e.g. bones; Calculus removers; Calculus smashing apparatus; Apparatus for removing obstructions in blood vessels, not otherwise provided for using mechanical vibrations, e.g. ultrasonic shock waves
    • A61B17/22012—Implements for squeezing-off ulcers or the like on the inside of inner organs of the body; Implements for scraping-out cavities of body organs, e.g. bones; Calculus removers; Calculus smashing apparatus; Apparatus for removing obstructions in blood vessels, not otherwise provided for using mechanical vibrations, e.g. ultrasonic shock waves in direct contact with, or very close to, the obstruction or concrement
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B17/00—Surgical instruments, devices or methods, e.g. tourniquets
    • A61B17/22—Implements for squeezing-off ulcers or the like on the inside of inner organs of the body; Implements for scraping-out cavities of body organs, e.g. bones; Calculus removers; Calculus smashing apparatus; Apparatus for removing obstructions in blood vessels, not otherwise provided for
    • A61B17/225—Implements for squeezing-off ulcers or the like on the inside of inner organs of the body; Implements for scraping-out cavities of body organs, e.g. bones; Calculus removers; Calculus smashing apparatus; Apparatus for removing obstructions in blood vessels, not otherwise provided for for extracorporeal shock wave lithotripsy [ESWL], e.g. by using ultrasonic waves
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/13—Tomography
    • A61B8/14—Echo-tomography
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/52—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/5207—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of raw data to produce diagnostic data, e.g. for generating an image
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/52—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/5269—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving detection or reduction of artifacts
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/54—Control of the diagnostic device
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B17/00—Surgical instruments, devices or methods, e.g. tourniquets
    • A61B2017/00017—Electrical control of surgical instruments
    • A61B2017/00022—Sensing or detecting at the treatment site
    • A61B2017/00106—Sensing or detecting at the treatment site ultrasonic
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B17/00—Surgical instruments, devices or methods, e.g. tourniquets
    • A61B2017/00017—Electrical control of surgical instruments
    • A61B2017/00137—Details of operation mode
    • A61B2017/00154—Details of operation mode pulsed
    • A61B2017/00172—Pulse trains, bursts, intermittent continuous operation
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61N—ELECTROTHERAPY; MAGNETOTHERAPY; RADIATION THERAPY; ULTRASOUND THERAPY
    • A61N7/00—Ultrasound therapy
    • A61N2007/0039—Ultrasound therapy using microbubbles
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61N—ELECTROTHERAPY; MAGNETOTHERAPY; RADIATION THERAPY; ULTRASOUND THERAPY
    • A61N7/00—Ultrasound therapy
    • A61N2007/0082—Scanning 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
    • 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/0091—Beam steering with moving parts, e.g. transducers, lenses, reflectors

Landscapes

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

Abstract

一种提供超声治疗的设备，包括：治疗超声换能器(ATA)，其适于生成聚焦的超声波；成像超声换能器(UID)，其与所述治疗超声换能器相关联；以及电子系统，其被配置为：控制治疗超声波换能器以便发射超声波能量的脉冲串、产生空化泡云(BC)；控制成像超声换能器，以采集待治疗区域的至少一个图像；获取由治疗超声换能器发射的超声波脉冲的多个回波信号，所述信号由所述成像超声换能器捕获；处理回波信号，以以重建空化泡云的图像IBC，该处理包括时空滤波；以及在待处理区域的图像上叠加显示所述空化泡云的图像。

Description

包括用于使空化泡成像的装置的超声处理设备
技术领域
本发明涉及一种用于通过超声提供对人体或动物体的治疗的设备。该设备特别适用于通过高强度聚焦超声(HIFU)、组织摧毁术、碎石术、血栓摧毁术等对移动器官(即，通常为腹部和胸部器官)进行治疗。
背景技术
组织摧毁术是一种通过聚焦超声脉冲来机械地摧毁组织的技术，该聚焦超声脉冲生成空化气泡/空化泡(为了简单起见，也称为“空化泡云”)。使用组织摧毁术治疗各种病变是可以代替高危外科手术的有前景的方法。空化允许高度受控的聚焦区域中的不期望的组织被非侵入性地破坏并且没有热效应。然而，人体是非常异质/不均匀的介质，并且在没有直接视觉反馈的情况下在体内进行治疗的事实使得精确地监测空化区域是至关重要的。特别地，在心脏应用中，胸腔可能导致超声的治疗路径的显著畸变。
对适用于在概念上类似于组织摧毁术的技术，诸如碎石术(肾结石的碎裂)和血栓摧毁术(血栓的碎裂)，或者甚至适用于通过高强度聚焦超声的治疗领域，都存在类似的考虑。
反向散射超声回波的频率的谐波的分析允许被动地检测空化，但不允许精确定位气泡云。
可以使用传统的B型超声成像技术，其允许实时地观察空化效应和被治疗的解剖结构。然而，这种方法并不完全令人满意。具体地，由于气泡在不均匀的生物组织内生成，所述可能难以在来自组织的所有回波中检测到它们的回波。因此，气泡的识别仍然是非常主观的、非定量的，并且对空化泡云的轮廓的定义成问题。
文章[1]教导了通过采用与时空滤波相关联的平行(或“超快速”)超声成像技术可以对空化泡进行更好辨别。然而却证明在被治疗的组织持续移动的情况下(诸如在心脏、肝脏或肾脏的情况下)，该方法并不是很适用。具体地，用于成像的超声波声穿透所有感兴趣区域(或者甚至比感兴趣区域大的区域)，并且它们的回波信号对于在该区域中发生的任何时空变化非常敏感，而不仅仅是由于空化引起的那些变化。这在组织是静止或准静止时尚可接受，但是在存在运动组织的情况下，时空滤波变成识别来自空化泡的回波信号的低得多的有效方式。然而，准确地说，在后一种情况下，能够精确地观察空化泡云及其相对于组织的位置是特别重要的。
被动探测方法也已经用于估计空化泡的位置。例如：
-[2]描述了一种用于声学映射由热消融治疗换能器生成的气泡的被动方法。该方法允许对长持续时间(几百振荡)的波列的连续波发射或发射所生成的气泡进行绘图。由于波列的持续时间，不可能精确地识别气泡云的位置。[2]提出了一种用于计算具有等于所使用的波长的分辨率的反向散射能量地图的方法，但这不允许精确地识别云的边界。
-[3]描述了[2]的替代方法，用于解决监测连续排放的相同问题。
-[4]描述了方法[2]在没有特定修改的情况下对组织摧毁的应用。
-[5]描述了一种对超声脉冲施加时间反转的方法。如在参考文献[2]、[3]和[4]中，该方法使用长脉冲并且不允许准确识别空化泡云的边界。
发明内容
本发明的目的是克服现有技术的上述缺点，并且更具体地提供一种成像技术，其允许对空化泡云进行精确定位和划分，包括在移动组织的治疗期间。
根据本发明，该目的通过将组织的主动超声成像、同步被动空化成像和时空滤波(例如，奇异值分解)相关联来实现。同步被动成像基于与治疗脉冲的发射序列同步地重建治疗超声波束与介质的相互作用的回波。其对于组织移动是稳健的，并且因此即使在成像的组织正在移动时也可以获得良好的性能。
本发明的技术允许以这样的空间分辨率定位空化泡云：该空间分辨率随着用于成像的超声脉冲(即超声治疗脉冲)的持续时间的减小而增大。因此，优选地，这些脉冲将具有落在0.1μs与50μs之间的持续时间(这在任何情况下都是聚焦治疗脉冲所期望的)。在参考文献[2]至[5]中描述的成像方法将不会从使用短脉冲获得任何益处。
在文献[6]中已经提出将时空滤波应用于被动空化成像，但是在完全不同的背景下，通过超声打开血脑屏障。该技术使用注入的微泡和强度不足以生成空化泡的超声。此外，通过头骨获得回波，这使得使用主动超声技术对治疗的组织进行成像几乎是不可能的。
因此，本发明的一个主题是用于提供超声治疗的设备，该设备包括：
-适于生成聚焦超声波的治疗超声换能器；
-与治疗超声换能器相关联的成像超声换能器；以及
-电子系统，其被配置为：
-控制治疗超声换能器以发射超声波能量的脉冲串，该脉冲串的持续时间适于在所述焦斑位于人体或动物体的待治疗区域内时导致在换能器的焦斑中生成空化泡云；
-控制成像超声换能器，以在由治疗超声换能器发射的两个超声波脉冲之间发射指向待治疗区域的超声波，以采集所述超声波的回波并对其进行处理以重建待治疗区域的至少一个图像；
-获取由治疗超声换能器发射的N>1个超声波脉冲的多个回波信号，所述信号由成像超声换能器捕获；
-使用波束成形算法来处理所获取的信号，以便形成相应的回波图像，并且使用时空滤波来提取回波图像的表示来自空化泡的超声波脉冲的反向散射的分量，将它们与代表来自待治疗区域的组织的反向散射的分量分离，以重建空化泡云的图像(这些操作可以按此顺序或相反的顺序进行)；以及
-在所述待处理区域的图像上叠加显示所述空化泡云的图像。
根据这种设备的特定实施例：
-所述电子系统还可以被配置为，在空化泡云的图像的重建期间，引入对治疗超声换能器和焦斑之间以及焦斑和成像超声换能器之间的超声波传播时间差的补偿。
-所述电子系统还可以被配置为控制该治疗超声换能器以便发射具有落在0.1μs与50μs之间并且优选地落在0.5μs与20μs之间的持续时间的超声波脉冲串。
-N可以落在2与10000之间并且优选在2与1000之间。
-所述电子系统还可以被配置为借助于并行波束成形算法来重建空化泡云的图像。
-所述电子系统和治疗超声波换能器可以被配置为以落在1和1000Hz之间的重复率发射具有落在0.1μs和50μs之间的持续时间、落在100kHz和5MHz之间的中心频率的超声波脉冲。
-所述电子系统和成像超声换能器可以被配置为发射具有落在2和15MHz之间的中心频率的超声波脉冲。
-所述电子系统还可以被配置为使用奇异值分解来执行所述时空滤波。
-治疗超声换能器和成像超声换能器可以同轴布置。
-该设备还可以包括用于在治疗区域中移动由治疗超声换能器发射的超声波脉冲的焦斑的设备。
-所述电子系统还可以被配置为根据空化泡云的图像来调节由治疗超声换能器发射的超声波脉冲的功率水平。
附图说明
通过阅读参考附图给出的描述，本发明的其他特征、细节和优点将变得显而易见，所述附图通过实例给出并且分别示出：
-图1是根据本发明的一个实施例的设备；
-图2是这种设备的功能示意图；
-图3是根据本发明的一个实施例的被动成像序列；
-图4是奇异值分解算法的图示；
-图5是使用根据本发明的一个实施例的设备获得的空化地图像；以及
-图6是说明本发明的技术效果的曲线图。
具体实施方式
图1的设备包括成像超声换能器UID，该成像超声换能器UID被布置在治疗换能器ATA的中心，该治疗换能器ATA由一组优选地数量落在5与20之间的环形的基本换能器组成，所述基本换能器被以对齐且同心的方式布置在聚焦超声波的球形碗SB上。成像换能器UID位于治疗换能器的最内环的中心。
换能器UID包括布置成一维(通常为线性)或二维、周期性或非周期性阵列的多个超声检测元件。例如，它可以是包括64个元件的阵列的二维回波描记探针。
成像换能器UID紧固到支架F，允许其连接到机械臂(未示出)，该机械臂通过机械连杆ML连接到治疗换能器，允许围绕轴线Ox的相对旋转运动(标记PJ为枢轴关节)。因此，可以在借助于马达(未示出)修改治疗换能器的取向(轴线Oz’)的同时保持成像换能器的取向(轴线Oz)恒定，从而允许治疗换能器的焦斑垂直于轴线Oz移动。此外，通过精细地控制各种基本换能器的控制信号之间的偏差，可以修改治疗换能器的焦距，并因此修改焦斑沿轴线Oz’的位置。这允许通过混合机械和电子设备来二维地扫描待治疗区域。
变型是可能的；例如，机械连杆还可以允许治疗换能器围绕第二轴线Oy旋转，该第二轴线同时垂直于Ox和Oz。
在欧洲专利申请EP 3236467 A1中更详细地描述了图1的组件。但是其仅作为示例给出，因为许多变型允许实施本发明。例如，也可以经由纯机械装置，或相反地纯电子装置(借助于使用基本换能器的矩阵阵列)来实施由治疗换能器的焦斑对待治疗的区域的扫描。此外，成像换能器与治疗换能器同轴、甚至与后者地机械连接均不是必要的：只要它们的相对位置是可控的，并且超声治疗脉冲的焦斑位于成像换能器的观察区域内部就足够了。
此外，成像换能器可以包括不同数量的元件，或者是双平面或矩阵阵列类型。
在任何情况下，必须提供电子系统SEL以控制成像和治疗换能器，以便：
-生成适于引发空化泡云形成的聚焦超声波脉冲；
-移动这些超声波的焦斑以便扫描待治疗的区域；
-获得并显示待处理区域和空化泡云的回波描记图像，并将它们显示在屏幕E上。
该电子系统SEL包括执行存储在存储器中的程序的一个或多个处理器、以及在这个或这些处理器的控制下操作的模拟和/或数字电子电路。图2示出了功能示意图。在该示意图中，标记PTD表示数据处理器；通常，它可以是包括微处理器或微控制器的板、计算机、或更复杂的可编程数字电子电路的组。该数据处理器通过接口设备IC(键盘、计算机鼠标等)从用户或操作者接收命令和/或参数；控制函数发生器GF的操作、用于控制成像换能器的电路CPTI的操作、波束形成电路GFF的操作、以及用于控制马达的电路CPM的操作；并处理从成像换能器的采集电路输出的信号以重建在屏幕E上显示的图像。
函数发生器GF生成被传递到电子系统的其他功能单元的多个电子信号。第一信号S1(例如频率为1MHz的正弦信号)被传送到波束形成电路GFF，该波束形成电路GFF将第一信号S1分解成多个具有不同相移的相同频率的单独信号；这些各种信号在被施加到成像换能器的基本换能器之前由功率放大器AP放大。如在上面引用的欧洲专利申请EP 3236467 A1中详细解释的，由波束形成电路GFF引入的相移允许治疗换能器的焦距被改变，并且因此超声焦斑TF沿着轴线Oz’的位置被改变。函数发生器还生成频率低得多(例如100Hz)的方波信号S2，其激活/去激活功率放大器；以此方式，该治疗换能器以100Hz的重复率发射超声脉冲(“治疗脉冲”)，这些脉冲例如具有8μs的持续时间。
函数发生器GF还生成高频(例如2MHz)的第三信号S3，该第三信号S3与信号S2一起被递送到用于电路CPTI以用于控制成像换能器。该电路控制成像换能器以在由治疗换能器发射的更强脉冲之间的间隔内发射低强度超声脉冲。这些脉冲的传播限定了观察区域RO，在二维成像换能器的情况下，该观察区域RO具有高度与Oz轴线一致的梯形形状。
成像换能器的获取电路CATI获取由成像换能器UID检测到的回波信号，将它们转换成数字格式，并将它们发送到处理器PTD，该处理器PTD继续处理它们。在跟随超声治疗脉冲的发射并且其持续时间(例如250μs)取决于焦斑TF的最大深度的第一时间窗口中，成像换能器检测治疗脉冲的回波，并且因此以被动模式操作；如下面将详细解释的，这些回波信号允许处理器PTD重建空化泡云的图像。在从第一窗口的尾端延伸到后续治疗脉冲的发射的第二时间窗口中，成像换能器检测其自身发射的脉冲的回波；这允许处理器PTD使用传统的主动超声成像技术来重建待治疗区域的解剖结构的图像。
在用于主动超声成像的脉冲的情况下，同一换能器用作超声源和检测器。这不是被动成像的情况，在被动成像中，超声脉冲由治疗换能器发射并且其回波由成像换能器检测。因此，为了能够重建空化泡云的图像——其图像以被动模式获取，需要知道“向外”传播时间(从治疗换能器到焦斑)与“返回”传播时间(从焦斑到成像换能器)之间的差异。如果足够精确地知道设备的机械配置，则这可以通过计算来获得，但是作为一般规则，优选地执行校准。为此，可以以多种不同的方式进行。
-首先，可以对应于焦斑放置水听器并测量由两个换能器发射的脉冲的到达时间。
-其次，可以对应于焦斑放置反射器，并在接收模式下使用成像换能器以测量由成像换能器自身和由治疗换能器发射的脉冲的回波的到达时间。
混合方法(使用水听器测量“向外”传播时间，使用反射器测量“返回”传播时间)也是可能的。
无论使用何种方法，都需要计算或测量与对整个待治疗区域进行采样的多个点相对应的传播时间差。
图3的时序图示出：以100Hz的重复率(10ms的周期)发射持续时间为8μs的治疗脉冲IUT，其生成空化泡的云BC；回波采集第一时间窗FT1，其从每个治疗脉冲开始起延伸250μs；以及主动成像第二时间窗FT2，其开始于每个第一窗的结束处并且延伸至下一治疗脉冲的发射。
马达控制电路CPM致动电动马达，允许治疗换能器枢转；它与波束形成电路GFF相互作用以移动超声脉冲的焦斑TF，从而扫描待治疗区域。在通过纯电子装置实现焦斑移动的实施例中可以省略马达控制电路CPM(相反，在其他实施例中，可以省略波束形成电路)。
上述各种功能单元(GF、CPTI、CATI、CPM、GFF、AP)不必对应于物理上分离的部件。例如，单个集成电路或电路板可以执行多个这些单元的全部或一些功能。相反，可以由多个集成电路和/或电路板来执行单个块的功能。
通过非限制性示例给出了脉冲特性。更一般地，成像换能器可以发射脉冲串，该脉冲串的持续时间典型地落在0.1μs和50μs之间(并且优选地落在0.5μs和20μs之间)，中心频率落在100kHz和5MHz之间，并且重复率落在1Hz和1000Hz之间，这些脉冲适于在焦点处生成落在50MPa和100MPa之间的峰值正压和落在-2.5MPa和-30MPa之间的峰值负压。成像换能器通常以高于治疗换能器的频率操作，典型地为落在2和15MHz之间的频率；例如，它可以发射2MHz的超声波，并且在接收模式中具有8MHz的采样频率。
如上所述，根据本发明的设备在主动模式和被动模式下都采集回波描记图像，主动模式是为了观察待治疗区域的组织，而被动模式是为了观察空化泡云。然后，图像被融合以便显示在屏幕E上，这允许视觉检查空化泡云相对于组织的位置。优选地，这些处理操作是实时执行的。在主动模式中获得的信号还允许估计空化所吸收的能量，并且因此在过程中评估治疗的有效性(和/或危险)。
处理器PTD因此可以自动控制治疗脉冲的强度，和/或在危险的情况下停止处理。例如，处理器可以逐渐增加超声脉冲的强度，直到探测到具有期望特性(形状、尺寸、能量等)的空化泡云，或者相反地，降低强度，同时确保空化的持续。该处理器还可以在已经施加了所希望的剂量时停止该处理，或者如果在所希望的位置中没有检测到空化泡云并且因此有损害有待治疗组织的风险时停止该处理。
主动成像信号可以以完全正常的方式被处理，因此不需要更详细地描述。相反，为了能够有效地提取空化泡的图像，处理器PTD必须对被动获取的信号进行特定的处理。该处理包括：
-获取N个通常是连续的治疗脉冲的回波信号，N通常落在2和10000之间，并且优选地落在2和1000之间；
-对这些回波信号应用波束形成重建算法，以重建N个图像；以及
-对如此重建的图像应用时空滤波，以便从中提取由空化泡导致的影响。
N的值必须严格地高于1，以便允许回波信号的时域滤波。可接受的上限取决于超声脉冲的重复率；在高重复率、每秒1000个脉冲并且N＝10000的情况下，每10秒获得一个滤波图像，这是对于大多数外科应用的最低可接受重复率。
优选地，使用本身已知的并行或超快类型的波束形成算法(不要与相对于治疗脉冲执行的波束形成混淆)来实现图像重建，但是该波束形成算法被修改为在其中包括对向外和返回路径上的传播时间差的补偿(这在上面讨论)。该算法基本上包括对由成像探头的各个元件检测到的信号的相干求和，这些信号在时间上偏移，以补偿超声波在向外和返回路径上的传播延迟(电子聚焦原理)。应当理解，考虑超声源和检测器之间的空间偏移对于有效聚焦是必要的。在频谱域中，通过信号相位的移位来实现对偏移的补偿。
空间-时间滤波可以使用任何盲源分离算法，其允许将源自空化泡的回波与那些尤其是从组织发出的回波区分开。例如，其可以是奇异值分解、(稀疏或独立的)主分量分析或非负矩阵因式分解等。下面，将考虑奇异向量分解的情况。
奇异值分解(SVD)是一种矩阵代数分解技术。它适用于图像的局部统计，并将最大能量集中到少量特征向量中。它包括将矩阵Xmxn因式分解为以下形式：U是mxm大小的正交矩阵、V是nxn正交矩阵、S是对角元素σi是X的奇异值并且其他元素是零的mxn矩阵。如果n<m：
Figure BDA0002795687750000081
其中V*是V的共轭转置矩阵，并且S的奇异值的数量等于X的秩。
当应用SVD方法进行图像处理时，最小奇异值与噪声相关联，并且图像的大部分能量被压缩成高值的奇异值。在本发明的情况下，SVD时空滤波用于通过明智地选择奇异向量将不期望的低频分量从与空化泡云相关联的高频分量中分离。这在图4中示出。
考虑与相应的治疗脉冲相关联并使用并行波束成形算法重建的一组N个“被动”图像F1、F2、…FN，每个图像具有(nx，ny)大小并被存储在缓冲存储器MT中。这些数据可被重新布置成称为Casora矩阵的(nx·ny，N)大小的二维时空矩阵X。该矩阵(图4中的标记SVD)的奇异值分解在于找到分别形成矩阵U和V的列的时间奇异向量和空间奇异向量，以及形成矩阵S的相应奇异值。有利的是，这些向量通过降低能量来排序。
实际的时空滤波在于仅使用描述气泡云的奇异向量来重建图像。假定这些向量与包含在p和q>p之间的索引的邻接奇异值相关联。过滤后的图像XBC因此由下式给出：
Figure BDA0002795687750000082
(图4中的参考标号RI)。
可以将滤波后的Casorati矩阵的项重新排列成三维矩阵XBC(nx、ny、N)，并且可以计算其功率积分或空化地图CM，其允许定位图像的最高能量区域：
CM(nx，nz)＝∫|XBC(nx，nz，N)|2 (等式3)
(图4中的参考标号IP)。
重要的是注意到N个图像给出单个空化地图；换句话说，获取速率被除以N。
将空化地图用作空化泡云的最终图像，该图像旨在叠加在通过主动回波描记术(图4中的标记IBC)获得的组织的图像上。可以使用伪色来表示地图CM的值
地图CM的空间分辨率受超声波脉冲在传播方向上的长度限制，该长度由其持续时间乘以其传播速度的乘积给出。因此，将优选地使用这样的脉冲：其尽可能短并且最短的持续时间等于超声波的一个周期。脉冲的持续时间意指用于控制治疗换能器的电子脉冲的持续时间(即，时间提示)，治疗换能器的有限带宽不可避免地导致实际发射的超声脉冲的延长。实际上，脉冲的持续时间不仅根据成像分辨率来选择，而且、甚至更重要的是根据治疗方案的需要来选择。
参数N、p和q的最优值的选择取决于所讨论的具体应用。发明人已经观察到，源于组织的回波信号(其是期望消除的)主要集中在第一奇异向量中，后面的N-1个奇异向量主要包含空化泡的贡献(期望被分离出来)。因此，p可以被设置为等于2(p＝2)并且q被设置为等于N(q＝N)。关于参数N，发明人观察到，如果该参数的值保持低于或等于约10，则对比度噪声比(下面定义)随N快速增加，然后趋于稳定。下面，考虑了三种情况：N＝6；N＝10和N＝14。
定义图像(更准确地说，空化地图)质量的对比度噪声比由下式给出：
Figure BDA0002795687750000091
其中，<>i表示区域i中的空间平均值，σi表示这两个区域的标准偏差(i＝1：气泡云；i＝2：背景)。将手动识别区域1和2。
作为变型，可以在波束形成之前实施时空滤波步骤。在这种情况下，滤波的空间分量涉及成像探头的元件。
为了测试本发明，用1.75L体积的水和8％聚乙烯醇(PVA)制备模型。使用1％纤维素(Sigmacell，20μm，USA)添加超声散射材料。使用实验室加热设备将自来水加热至90℃，并将所需体积倒入具有磁性转向枢轴销的烧杯中。然后将PVA溶解在水中。将混合物冷却至40℃、加入纤维素，将溶液倒入由塑料制成的方形容器中，并在冰箱中放置8小时。然后将模型解冻，放回冷冻机中8小时。在该过程结束时，将模型脱模并放置在水箱中用于空化实验。
根据本发明的设备的换能器组件附接到双轴台(PI，Micos，德国)并放置在含有模型的水箱内。为了模拟待治疗的身体区域的生理运动，以不同的速度引起换能器的轴向和横向运动，对于轴向运动，最高速度为10mm/s、运动幅度为10mm，对于横向运动，最高速度为10mm/s、运动幅度为8mm。
图5示出了以这种方式获得的空化泡云的图像。云由白盒标识。伪影存在于图像的上方和下方；它们最可能是由于超声脉冲从水-模型界面的反射。
在静态条件下和在换能器移动的情况下，测量对于N＝6、10和14(再次p＝2和q＝N)获得的空化泡云的“被动”图像的对比度噪声比CNR。还使用主动成像技术进行类似的测量。本发明人已经观察到，在静态条件下，空化泡的主动成像给出了比被动成像略好的结果，但是在换能器移动的情况下相反。这证实了本发明中的被动方法的使用。
换能器移动时获得的结果在图6中示出，其中，参考标记P表示被动图像的CNR，A1-A5表示使用主动方法获得的图像的CNR：
-A1：通过在每个治疗脉冲之后发射平面超声波获得的主动成像。
-A2：超快序列，在每个治疗脉冲之后发射十一个定向发散波，合成单个图像。
-A3：超快序列，在每个治疗脉冲之后发射10个发散波(没有组合)，这导致每10毫秒10个图像。
-A4：超快序列，发射分布在三个治疗脉冲上的、一百十次定向发散波，其被组合成十个图像。每个图像是十一个定向发散波的组合的结果。该序列大约每35毫秒给出十个图像。
实际上，所观察到的是，换能器的移动引起了主动成像技术的性能的非常显著的退化(一个数量级)，而根据本发明获得的图像的CNR仅降低大约两倍。这可以通过以下方式来解释。当要治疗的区域正在(或者，等效地，换能器正在)移动时，来自组织的反射包含干扰从空化泡反射的有用分量的高频非相干时空分量。然而，在主动成像中，整个观察区域中的组织暴露于超声并因此生成回波。相反，在被动成像中，使用聚焦超声，这降低了由组织引起的寄生影响。
参考文献：
[1]B.Arnal,J.Baranger,C.Demene,M.Tanter and M.Pernot(2016).“运动器官中的体内实时空化成像(In vivo real-time cavitation imaging in moving organs)”。Phys Med Biol.2017Feb 7；62(3):843-857.
[2]Christian Coviello,Richard Kozick,James Choi,Miklós
Figure BDA0002795687750000101
CarlJensen,Penny Probert Smith,and Constantin-C.Coussios(2015).“超声治疗监测中利用最优波束形成的被动声标测(Passive acoustic mapping utilizing optimalbeamforming in ultrasound therapy monitoring)”。The Journal of the AcousticalSociety of America 137,2573.
[3]Kevin J.Haworth,Kenneth B.Bader,Kyle T.Rich,Christy K.Holland,Member and T.Douglas Mast(2017).“定量频域被动空化成像(Quantitative Frequency-Domain Passive Cavitation Imaging)”。“关于超声、铁电和频率控制的IEEE处理(IEEEtransactions on ultrasonics,ferroelectrics,and frequency control)”,vol.64,no.1.
[4]Bader Kenneth,Haworth Kevin,D.Maxwell Adam,Holland Christy(2017).“被动空化成像在体外组织摧毁术液化分类中的应用(Post Hoc Analysis of PassiveCavitation Imaging for Classification of Histotripsy-Induced Liquefaction inVitro)”。IEEE Transactions on Medical Imaging.
[5]P.Boulos,F.Varray,A.Poizat,J.C.Bera,C.Cachard(2015).“使用开放式超声系统和时间反转重建的被动空化成像(Passive cavitation imaging using an openultrasonic system and time reversal reconstruction)”。22ème Congrès
Figure BDA0002795687750000102
deMécanique Lyon.

Claims (13)

1.一种用于提供超声治疗的设备，包括：
-治疗超声换能器，其适于生成聚焦的超声波；
-成像超声换能器，其与所述治疗超声换能器相关联；以及
-电子系统，其被配置为：
-控制所述治疗超声换能器，以便发射超声波能量的脉冲串，所述超声波能量的脉冲串的持续时间适于在所述治疗超声换能器的焦斑位于人体或动物体的待治疗区域内时导致在所述焦斑中产生空化泡云；
-控制所述成像超声换能器，以在由所述治疗超声换能器发射的两个超声波脉冲之间发射指向所述待治疗区域的超声波，以采集所述超声波的回波并对其进行处理以重建所述待治疗区域的至少一个图像；
-获取由所述治疗超声换能器发射的N>1个超声波脉冲的多个回波信号，所述信号由所述成像超声换能器捕获；
-使用波束成形算法来处理由所述治疗超声换能器发射的N个超声波脉冲的所述回波信号，以便形成相应的回波图像，并且使用时空滤波来提取所述回波图像的代表来自所述空化泡云的超声波脉冲的反向散射的分量、将之与代表来自所述待治疗区域的组织的反向散射的分量分离，以重建所述空化泡云的图像；以及
-在所述待治疗区域的图像上叠加显示所述空化泡云的图像。
2.根据权利要求1所述的用于提供超声治疗的设备，其中所述电子系统还被配置为：在所述空化泡云的图像的重建过程中，引入对超声波在所述治疗超声换能器与所述焦斑之间、以及在所述焦斑与所述成像超声换能器之间的传播时间差的补偿。
3.根据权利要求1或2所述的用于提供超声治疗的设备，其中所述电子系统被配置为：控制所述治疗超声换能器以便发射具有落在0.1 µs与50 µs之间的持续时间的超声波脉冲串。
4.根据权利要求1或2所述的用于提供超声治疗的设备，其中所述电子系统被配置为：控制所述治疗超声换能器以便发射具有落在0.5 µs与20 µs之间的持续时间的超声波脉冲串。
5.根据权利要求1或2所述的用于提供超声治疗的设备，其中N落在2与10000之间。
6.根据权利要求1或2所述的用于提供超声治疗的设备，其中N落在2与1000之间。
7.根据权利要求1或2所述的用于提供超声治疗的设备，其中所述电子系统被配置为：借助于并行波束成形算法来重建所述空化泡云的图像。
8.根据权利要求1或2所述的用于提供超声治疗的设备，其中，所述电子系统和所述治疗超声换能器被配置为：以落在1 Hz和1000 Hz之间的重复率发射持续时间落在0.1 µs和50 µs之间、中心频率落在100 kHz和5 MHz之间的超声波脉冲。
9.根据权利要求1或2所述的用于提供超声治疗的设备，其中所述电子系统和所述成像超声换能器被配置为：发射中心频率落在2 MHz与15 MHz之间的超声波脉冲。
10.根据权利要求1或2所述的用于提供超声治疗的设备，其中所述电子系统被配置为：使用奇异值分解来执行所述时空滤波。
11.根据权利要求1或2所述的用于提供超声治疗的设备，其中所述治疗超声换能器和所述成像超声换能器是同轴布置的。
12.根据权利要求1或2所述的用于提供超声治疗的设备，还包括用于在所述治疗区域中移动由所述治疗超声换能器发射的超声波脉冲的焦斑的装置。
13.根据权利要求1或2所述的用于提供超声治疗的设备，其中所述电子系统还被配置为：根据所述空化泡云的图像来调节由所述治疗超声换能器发射的超声波脉冲的功率水平。
CN201980035031.1A 2018-05-25 2019-05-24 包括用于使空化泡成像的装置的超声处理设备 Active CN112566694B (zh)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
FR1854428 2018-05-25
FR1854428A FR3081334B1 (fr) 2018-05-25 2018-05-25 Appareil de traitement par ultrasons comportant des moyens d'imagerie des bulles de cavitation
PCT/EP2019/063455 WO2019224350A1 (fr) 2018-05-25 2019-05-24 Appareil de traitement par ultrasons comportant des moyens d'imagerie des bulles de cavitation

Publications (2)

Publication Number Publication Date
CN112566694A CN112566694A (zh) 2021-03-26
CN112566694B true CN112566694B (zh) 2022-09-27

Family

ID=63637977

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201980035031.1A Active CN112566694B (zh) 2018-05-25 2019-05-24 包括用于使空化泡成像的装置的超声处理设备

Country Status (9)

Country Link
US (1) US11571230B2 (zh)
EP (1) EP3801761B1 (zh)
JP (1) JP7358391B2 (zh)
KR (1) KR20210027249A (zh)
CN (1) CN112566694B (zh)
CA (1) CA3101381A1 (zh)
FR (1) FR3081334B1 (zh)
IL (1) IL278882A (zh)
WO (1) WO2019224350A1 (zh)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US11813484B2 (en) 2018-11-28 2023-11-14 Histosonics, Inc. Histotripsy systems and methods
EP4096782A4 (en) 2020-01-28 2024-02-14 Univ Michigan Regents SYSTEMS AND METHODS FOR IMMUNOSENSITIZATION BY HISTOTRIPSY
CN112451865A (zh) * 2020-11-24 2021-03-09 北京小超科技有限公司 超声治疗监测装置及方法
JP7137682B1 (ja) 2021-11-29 2022-09-14 ソニア・セラピューティクス株式会社 超音波治療装置
JP7122483B1 (ja) 2022-01-18 2022-08-19 ソニア・セラピューティクス株式会社 Ｈｉｆｕ照射装置用の超音波イメージングプローブおよび超音波画像表示装置
CN115500865A (zh) * 2022-10-08 2022-12-23 珠海医凯电子科技有限公司 用于被动声成像的快速延迟乘和自适应波束形成方法

Family Cites Families (20)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE69634714T2 (de) * 1995-03-31 2006-01-19 Kabushiki Kaisha Toshiba, Kawasaki Therapeutisches Ultraschallgerät
CN101152646B (zh) * 2006-09-27 2012-07-04 香港理工大学 柔性超声换能器阵列及其应用装置
WO2008062342A2 (en) * 2006-11-20 2008-05-29 Koninklijke Philips Electronics, N.V. Control and display of ultrasonic microbubble cavitation
GB0820377D0 (en) * 2008-11-07 2008-12-17 Isis Innovation Mapping and characterization of cavitation activity
CN101623203B (zh) * 2009-08-07 2011-02-09 西安交通大学 一种瞬态物理过程中多模式和多参量同步检测成像监控系统及监控方法
AU2010289769B2 (en) * 2009-08-26 2016-06-30 Histosonics, Inc. Micromanipulator control arm for therapeutic and imaging ultrasound transducers
CN103028203B (zh) * 2011-10-09 2016-08-31 北京汇福康医疗技术股份有限公司 超声减脂治疗头及其工作方法
CN103284753B (zh) * 2012-02-22 2015-12-09 香港理工大学 超声波成像系统及成像方法
US10624612B2 (en) * 2014-06-05 2020-04-21 Chikayoshi Sumi Beamforming method, measurement and imaging instruments, and communication instruments
WO2016030889A1 (en) * 2014-08-27 2016-03-03 Fus Mobile Inc. Handheld devices for projecting focused ultrasound and related methods
CN104622504A (zh) * 2015-02-12 2015-05-20 黄晶 一种相控阵高强度聚焦超声检测、治疗系统
JP6943768B2 (ja) * 2015-03-30 2021-10-06 コーニンクレッカ フィリップス エヌ ヴェＫｏｎｉｎｋｌｉｊｋｅ Ｐｈｉｌｉｐｓ Ｎ．Ｖ． 超音波血栓溶解処置およびモニタリングのための超音波トランスデューサ・アレイ
KR102574559B1 (ko) * 2015-04-02 2023-09-05 카르디아웨이브 판막증을 치료하기 위한 방법 및 장치
CN104815399B (zh) * 2015-04-03 2018-04-17 西安交通大学 基于六轴机械臂的高强度聚焦超声治疗引导和控制系统及方法
CN105054972A (zh) * 2015-07-20 2015-11-18 东南大学 一种超声诊疗系统及利用该系统治疗肿瘤的方法
EP3236467A1 (en) * 2016-04-22 2017-10-25 Cardiawave Ultrasound imaging and therapy device
CN106037815A (zh) * 2016-05-17 2016-10-26 西安交通大学 一种用于热凝固监测的超声回波统计参量成像系统与方法
CN106063712B (zh) * 2016-05-24 2019-06-28 黄晶 一种超声诊疗方法和系统
CN106419955B (zh) * 2016-09-07 2019-08-13 苏州国科昂卓医疗科技有限公司 超声波束合成方法的应用及剪切波超声弹性成像方法
US11369810B2 (en) * 2016-12-19 2022-06-28 Michalakis Averkiou Method and apparatus for ultrasonic mediation of drug delivery using microbubbles

Also Published As

Publication number Publication date
IL278882A (en) 2021-01-31
US20210212709A1 (en) 2021-07-15
FR3081334A1 (fr) 2019-11-29
WO2019224350A1 (fr) 2019-11-28
US11571230B2 (en) 2023-02-07
CN112566694A (zh) 2021-03-26
JP7358391B2 (ja) 2023-10-10
EP3801761B1 (fr) 2022-07-06
CA3101381A1 (en) 2019-11-28
KR20210027249A (ko) 2021-03-10
JP2021523776A (ja) 2021-09-09
FR3081334B1 (fr) 2020-05-01
EP3801761A1 (fr) 2021-04-14

Similar Documents

Publication Publication Date Title
CN112566694B (zh) 包括用于使空化泡成像的装置的超声处理设备
OˈReilly et al. A super‐resolution ultrasound method for brain vascular mapping
EP3328285B1 (en) A method and system for correcting fat-induced aberrations
US6334846B1 (en) Ultrasound therapeutic apparatus
US6280402B1 (en) Ultrasound therapeutic apparatus
JP5707148B2 (ja) 医用画像診断装置及び医用画像処理装置
Soulioti et al. Super-resolution imaging through the human skull
Hall et al. A real-time measure of cavitation induced tissue disruption by ultrasound imaging backscatter reduction
US20100249597A1 (en) Method and system for imaging vessels
KR20180013956A (ko) 단일 추적 위치 전단파 탄성 이미징을 위한 방법, 시스템 및 컴퓨터 프로그램 제품
US20160174945A1 (en) Image processing apparatus, medical image apparatus and image fusion method for the medical image
JP2007181698A (ja) 病変を検出するための超音波診断システム
KR20160062874A (ko) 프로브, 초음파 영상 장치 및 초음파 영상 장치의 제어 방법
Escudero et al. 2D and 3D real-time passive cavitation imaging of pulsed cavitation ultrasound therapy in moving tissues
JPWO2007032134A1 (ja) 超音波診断装置
US10238369B2 (en) Real time ultrasound thermal dose monitoring system for tumor ablation therapy
JP5121384B2 (ja) 超音波診断装置
CN107073287B (zh) 放射治疗操作的超声波图像引导
KR20200043048A (ko) 초음파 영상 장치 및 그 제어 방법
Fahey et al. ARFI imaging of thermal lesions in ex vivo and in vivo soft tissues
Zheng et al. A targeting method based on acoustic backscatter for treatment planning in tissue ablation using focused ultrasound
JP5159885B2 (ja) 超音波診断装置
Shiozaki et al. 1Pa5-6 Ultrasound Imaging of Cavitation Bubbles by Triplet Pulse Sequence with Reduction of Therapeutic Ultrasound Noise
Ebbini et al. Quadratic b-mode and pulse inversion imaging of thermally-induced lesions in vivo
Escudero New echocardiographic imaging modalities for guiding and monitoring of therapy by focused ultrasound of calcified aortic stenosis

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