CN113226208B - 用于监测组织温度的方法和系统 - Google Patents

用于监测组织温度的方法和系统 Download PDF

Info

Publication number
CN113226208B
CN113226208B CN201980086822.7A CN201980086822A CN113226208B CN 113226208 B CN113226208 B CN 113226208B CN 201980086822 A CN201980086822 A CN 201980086822A CN 113226208 B CN113226208 B CN 113226208B
Authority
CN
China
Prior art keywords
interest
tissue
temperature
peak
acoustic signal
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
CN201980086822.7A
Other languages
English (en)
Other versions
CN113226208A (zh
Inventor
曹骏桓
保罗·A·皮科特
迈克尔·M·桑顿
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.)
Andra Life Sciences
Original Assignee
Andra Life Sciences
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 Andra Life Sciences filed Critical Andra Life Sciences
Publication of CN113226208A publication Critical patent/CN113226208A/zh
Application granted granted Critical
Publication of CN113226208B publication Critical patent/CN113226208B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/01Measuring temperature of body parts ; Diagnostic temperature sensing, e.g. for malignant or inflamed tissue
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/0093Detecting, measuring or recording by applying one single type of energy and measuring its conversion into another type of energy
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/0093Detecting, measuring or recording by applying one single type of energy and measuring its conversion into another type of energy
    • A61B5/0095Detecting, measuring or recording by applying one single type of energy and measuring its conversion into another type of energy by applying light and detecting acoustic waves, i.e. photoacoustic measurements
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/01Measuring temperature of body parts ; Diagnostic temperature sensing, e.g. for malignant or inflamed tissue
    • A61B5/015By temperature mapping of body part
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/05Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves 
    • A61B5/0507Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves  using microwaves or terahertz waves
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/48Other medical applications
    • A61B5/4836Diagnosis combined with treatment in closed-loop systems or methods
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/48Other medical applications
    • A61B5/4869Determining body composition
    • A61B5/4872Body fat
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/72Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7271Specific aspects of physiological measurement analysis
    • A61B5/7278Artificial waveform generation or derivation, e.g. synthesising signals from measured signals
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/74Details of notification to user or communication with user or patient ; user input means
    • A61B5/7475User input or interface means, e.g. keyboard, pointing device, joystick
    • A61B5/748Selection of a region of interest, e.g. using a graphics tablet
    • A61B5/7485Automatic selection of region of interest
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • 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/5238Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data for combining image data of patient, e.g. merging several images from different acquisition modes into one image
    • A61B8/5261Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data for combining image data of patient, e.g. merging several images from different acquisition modes into one image combining images from different diagnostic modalities, e.g. ultrasound and X-ray
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B18/00Surgical instruments, devices or methods for transferring non-mechanical forms of energy to or from the body
    • A61B18/02Surgical instruments, devices or methods for transferring non-mechanical forms of energy to or from the body by cooling, e.g. cryogenic techniques
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B18/00Surgical instruments, devices or methods for transferring non-mechanical forms of energy to or from the body
    • A61B18/04Surgical instruments, devices or methods for transferring non-mechanical forms of energy to or from the body by heating
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B18/00Surgical instruments, devices or methods for transferring non-mechanical forms of energy to or from the body
    • A61B2018/00636Sensing and controlling the application of energy
    • A61B2018/00773Sensed parameters
    • A61B2018/00791Temperature
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B2576/00Medical imaging apparatus involving image processing or analysis
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/20Measuring for diagnostic purposes; Identification of persons for measuring urological functions restricted to the evaluation of the urinary system
    • A61B5/201Assessing renal or kidney functions
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/42Detecting, measuring or recording for evaluating the gastrointestinal, the endocrine or the exocrine systems
    • A61B5/4222Evaluating particular parts, e.g. particular organs
    • A61B5/4244Evaluating particular parts, e.g. particular organs liver
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/45For evaluating or diagnosing the musculoskeletal system or teeth
    • A61B5/4519Muscles
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/48Other medical applications
    • A61B5/4887Locating particular structures in or on the body
    • A61B5/489Blood vessels
    • 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/899Combination of imaging systems with ancillary equipment

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Veterinary Medicine (AREA)
  • Public Health (AREA)
  • Pathology (AREA)
  • Biophysics (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Medical Informatics (AREA)
  • Molecular Biology (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • General Health & Medical Sciences (AREA)
  • Physiology (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Radiology & Medical Imaging (AREA)
  • Gastroenterology & Hepatology (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Dentistry (AREA)
  • Oral & Maxillofacial Surgery (AREA)
  • Orthopedic Medicine & Surgery (AREA)
  • Rheumatology (AREA)
  • Endocrinology (AREA)
  • Urology & Nephrology (AREA)
  • Vascular Medicine (AREA)
  • Artificial Intelligence (AREA)
  • Psychiatry (AREA)
  • Signal Processing (AREA)
  • Acoustics & Sound (AREA)
  • Ultra Sonic Daignosis Equipment (AREA)

Abstract

一种用于监测关注区域内的组织温度的热声成像系统和方法,该关注区域具有由至少一个边界隔开的关注对象和参考物。该系统和方法包括:热声成像系统以及一个或更多个处理器,所述热声成像系统具有:可调节射频(RF)施加器,其被配置为将RF能量脉冲发射到组织的关注区域中并加热其中的组织;声学接收器,其被配置为接收响应于关注区域中的组织的加热而生成的双极声学信号;所述一个或更多个处理器处理响应于RF能量脉冲在关注区域中生成的至少一个接收到的双极声学信号以确定其峰峰幅度以及使用至少一个双极声学信号的峰峰幅度计算在至少一个边界处的温度。

Description

用于监测组织温度的方法和系统
领域
本主题公开涉及热声成像,并且特别涉及一种用于监测组织温度的方法和系统。
背景
在热疗或冷冻治疗期间,需要监测组织中的温度分布,以确保热能在周围健康组织中的安全沉积并有效破坏肿瘤和异常细胞。为此,需要具有空间分辨率(约1毫米)和高温度灵敏度(1开氏度或更好)的实时温度监测。
最准确的温度监测是通过用热电偶或热敏电阻直接测量温度。然而,直接测量温度是侵入性的,通常不是优选的,并且根本不可行。
已经开发了几种非侵入式温度监测方法。例如,红外热成像已被用于监测组织温度。尽管红外热成像具有0.1摄氏度的准确度,但它仅限于表面温度。
超声波已被用于监测组织温度。虽然超声波具有良好的空间分辨率和较高的穿透深度,但超声波的温度灵敏度较低。
磁共振成像已被用于监测组织温度。虽然磁共振成像具有高分辨率和灵敏度的优点,但磁共振成像价格昂贵、体积大并且缓慢。
尽管用于监测组织温度的技术已经被考虑到,但仍需要改进。因此,至少一个目的是提供一种用于监测组织温度的新方法和系统。
概述
应当理解,提供本概述是为了以简化的形式介绍在下面的详细描述中进一步被描述的概念的选择。本概述无意用于限制要求保护的主题的范围。
因此,在一个方面,提供了一种用于估计组织温度的方法,该方法包括(i)使用射频(RF)施加器将一个或更多个RF能量脉冲引导到组织的关注区域中,该组织的关注区域包括由至少一个边界隔开的关注对象和至少一个参考物;(ii)使用声学接收器检测响应于RF能量脉冲而在关注区域中生成的至少一个双极声学信号,并处理该至少一个双极声学信号以确定该至少一个双极声学信号的峰峰幅度;以及(iii)使用该至少一个双极声学信号的峰峰幅度计算在至少一个边界处的温度。
在一个或更多个实施例中,该方法包括基于计算的在至少一个边界处的温度来估计关注对象的温度。
在一个或更多个实施例中,该方法包括使用声学接收器检测由第二关注对象响应于RF能量脉冲而生成的至少一个另外的双极声学信号,并处理该至少一个另外的双极声学信号以确定该至少一个另外的双极声学信号的峰峰幅度。
在一个或更多个实施例中,该方法包括比较双极声学信号的峰峰幅度、计算的在至少一个边界处的温度和另外的双极声学信号的峰峰幅度以估计第二关注对象的温度。
在一个或更多个实施例中,边界位于至少两个不同类型的组织之间的位置处。
在一个或更多个实施例中,两个不同类型的组织是以下中的一种:肌肉和脂肪;血管和脂肪;以及肝组织和肾组织。
在一个或更多个实施例中,该方法包括在治疗期间重复步骤(i)至(iii)以连续监测组织温度。
根据另一方面,提供了一种用于确定关注区域内的组织温度的系统,该关注区域包括由至少一个边界隔开的关注对象和参考物,该系统包括热声成像系统以及一个或更多个处理器,该热声成像系统包括:可调射频(RF)施加器,其被配置为将RF能量脉冲发射到组织的关注区域并加热在该组织的关注区域中的组织;以及声学接收器,其被配置为接收响应于关注区域中的组织的加热而生成的双极声学信号;所述一个或更多个处理器被配置为:处理响应于RF能量脉冲而在关注区域中生成的至少一个接收到的双极声学信号以确定该至少一个接收到的双极声学信号的峰峰幅度;以及使用至少一个双极声学信号的峰峰幅度计算在至少一个边界处的温度。
在一个或更多个实施例中,该一个或更多个处理器还被配置为基于计算的在至少一个边界处的温度来估计关注对象的温度。
在一个或更多个实施例中,该一个或更多个处理器还被配置为使用声学接收器检测由第二关注对象响应于RF能量脉冲而生成的至少一个另外的双极声学信号,并处理该至少一个另外的双极声学信号以确定该至少一个另外的双极声学信号的峰峰幅度。
在一个或更多个实施例中,该一个或更多个处理器还被配置为比较双极声学信号的峰峰幅度、计算的在至少一个边界处的温度和另外的双极声学信号的峰峰幅度以估计第二关注对象的温度。
根据另一方面,提供了一种监测组织温度的方法,该方法包括(i)使用射频(RF)施加器将一个或更多个RF能量脉冲引导到组织的关注区域中,该组织的关注区域包括由至少一个边界隔开的关注对象和至少一个参考物;(ii)使用声学接收器检测响应于RF能量脉冲而在关注区域中生成的至少一个双极声学信号,并处理该至少一个双极声学信号以确定该至少一个双极声学信号的峰峰幅度;(iii)使用该至少一个双极声学信号的峰峰幅度计算在至少一个边界处的温度;(iv)重复所述引导、检测和计算以监测在至少一个边界处的温度。
附图简述
现在将参考附图更全面地描述实施例,其中:
图1是成像系统的示意图;
图2是示出示例性双极声学信号的曲线图;
图3是示出示例性电场强度衰减曲线的曲线图;
图4是用于估计组织温度的方法的流程图;
图5是包含关注对象和参考的示例性关注区域;
图6是用于估计组织温度的另一方法的流程图;
图7示出可以使用图1的成像系统成像的人体的各个部分;以及
图8是示出Grüneisen参数如何随温度变化的曲线图。
具体实施方式
当结合附图阅读时,将更好地理解前述概述以及某些示例的以下详细描述。如本文所用,以单数形式引入且前面有词“一(a)”或“一(an)”的元素或特征应理解为不一定排除多个元素或特征。此外,对“一个示例”或“一个实施例”的引用并不旨在被解释为排除也包含所描述的元素或特征的另外的示例或实施例的存在。此外,除非有明确相反地指出,否则“包括(comprising)”或“具有(having)”或“包括(including)”具有特定特性的一个元素或特征或多个元素或特征的示例或实施例可以包括不具有该特性的另外的元素或特征。此外,应当理解,术语“包括(comprising)”、“具有(having)”、“包括(including)”是指“包括但不限于”并且术语“包括(comprising)”、“具有(having)”和“包括(including)”具有等效含义。
如本文所用,术语“和/或”可以包括一个或更多个相关联的所列元素或特征的任何和所有组合。
将要理解的是,当元素或特征被称为“在另一元素或特征上(on)”、“附接(attached)”到另一元素或特征、“连接(connected)”到另一元素或特征、与另一元素或特征“耦合(coupled)”、“接触(contacting)”另一元素或特征等时,该元素或特征可直接在另一元素或特征上、附接到另一元素或特征、连接到另一元素或特征、与另一元素或特征耦合或接触另一元素或特征,或者也可以存在中间元素。相比之下,当元素或特征被称为例如“直接在另一元素或特征上”、“直接附接”到另一元素或特征、“直接连接”到另一元素或特征、与另一元素或特征“直接耦合”或“直接接触”另一元素或特征时,不存在中间元素或特征。
将要理解的是,空间相对术语,诸如“在……下面(under)”、“在……下方(below)”、“下部(lower)”、“在……上方(over)”、“在……上面(above)”、“上部(upper)”、“前方(front)”、“后方(back)”等,在本文可能为了便于描述而被用于描述在附图中所示的元素或特征与另一元素或特征的关系。然而,除了图中描绘的取向之外,在使用或操作中空间相对术语还可以包括不同取向。
在下文中,描述了用于估计组织温度的方法和系统。总体上,该方法和系统利用RF施加器来获得在受试者的关注区域(ROI)内的组织的热声数据。分析热声数据以估计组织温度。
现在转向图1,示出了示例性成像系统并且通常由附图标记20标识。可以看出,成像系统20包括通信地耦合到超声成像系统24和热声成像系统26的编程计算设备22。超声成像系统24和热声成像系统26被配置为分别获得与受试者S相关联的组织的关注区域ROI的超声图像数据和热声图像数据。
在该实施例中的编程计算设备22是个人计算机或其它合适的处理设备,例如,其包括:包括一个或更多个处理器的处理单元、系统存储器(易失性和/或非易失性存储器)、其它不可移动或可移动存储器(例如,硬盘驱动器、RAM、ROM、EEPROM、CD-ROM、DVD、闪存等)以及将各种计算机部件耦合到处理单元的系统总线。计算设备22还可以包括使用以太网、Wi-Fi和/或其它合适的网络格式的联网能力,以能够连接到共享或远程驱动器、一个或更多个联网计算机或其它联网设备。一个或更多个输入设备,诸如鼠标和键盘(未示出),耦合到计算设备22,以用于接收操作者输入。显示设备(未示出),诸如一个或更多个计算机屏幕或监视器,耦合到计算设备22,以用于显示一个或更多个生成的图像,该一个或更多个生成的图像基于从超声成像系统24接收的超声图像数据和/或从热声成像系统26接收的热声图像数据。
超声成像系统24包括超声换能器28形式的声学接收器,该超声换能器安置一个或更多个超声换能器阵列30,该一个或更多个超声换能器阵列30被配置为将声波发射到受试者S的关注区域ROI中。被引导到受试者的关注区域ROI内的声波在关注区域ROI产生回声,不同的组织反射不同程度的声音。由一个或更多个超声换能器阵列30接收的回声在作为超声图像数据传送到计算设备22进行进一步处理和呈现为可由操作者解释的超声图像之前由超声成像系统24进行处理。在该实施例中,超声成像系统24利用假设1540m/s的标称声速的B模式(B-mode)超声成像技术。由于超声成像系统在本领域中是已知的,因此本文将不再进一步描述超声成像系统24的另外的细节。
热声成像系统26包括以热声换能器32形式的声学接收器。热声换能器32安置一个或更多个热声换能器阵列34以及射频(RF)施加器36。然而,应当理解的是,RF施加器36可以与热声换能器32分开安置。RF施加器36被配置为发射短脉冲的RF能量,该短脉冲的RF能量被引导到受试者的关注区域ROI内的组织中。在该实施例中,RF施加器36具有大约10MHz到100GHz之间的频率并且具有大约0.1纳秒到10纳秒之间的脉冲持续时间。传递到关注区域ROI内的组织的RF能量脉冲加热组织,从而诱发被热声换能器32检测到的声压波(acousticpressure waves)。被热声换能器32检测到的声压波被处理并作为热声图像数据传送到计算设备22,以用于进一步处理和呈现为可由操作者解释的热声图像。
在该实施例中,超声换能器28和热声换能器32机械互连,使得在一个或更多个超声换能器阵列30、一个或更多个热声阵列34和RF施加器36之间的空间关系是已知的。使用一个或更多个超声换能器阵列30、一个或更多个热声换能器阵列34和RF施加器36的中心线来设置空间关系。每条中心线被定义为相应换能器阵列的区域的正中央(mid-point)。
在该实施例中,一个或更多个超声换能器阵列30与一个或更多个热声换能器阵列34之间的空间关系是使得一个或更多个热声换能器阵列34的中心线以相对于一个或更多个超声换能器阵列30的中心线(也称为轴向轴或超声换能器阵列束轴)的已知角度α设置。一个或更多个热声换能器阵列34和RF施加器36之间的空间关系是使得RF施加器36的中心线与一个或更多个热声换能器阵列34的中心线间隔开并且与之大体上平行。
成像系统20利用一个或更多个超声换能器阵列30与一个或更多个热声换能器阵列34之间的已知空间关系来增加热声成像的精度和准确度。
超声换能器28的一个或更多个超声换能器阵列30的坐标系和热声换能器32的一个或更多个热声换能器阵列34的坐标系由计算设备22映射,使得获取的超声和热声图像可以被记录(registered)。可替代地,热声成像系统26可通过将一个或更多个超声换能器阵列30与超声换能器28断开连接并将一个或更多个超声换能器阵列30连接到热声换能器32来利用超声换能器28的一个或更多个超声换能器阵列30。如将要理解的,通过这样做,不需要在一个或更多个超声换能器阵列30与一个或更多个热声换能器阵列34之间的坐标映射。
在包括在脂肪或脂肪组织与软组织或瘦组织之间的边界的关注区域的热声成像期间,生成由热声换能器32接收的双极声学信号。这是由于如下事实:软组织或瘦组织比脂肪或脂肪组织吸收更多的热量,导致软组织或瘦组织迅速膨胀越过边界并进入脂肪或脂肪组织,该脂肪或脂肪组织膨胀较少,并且然后迅速收缩。双极声学信号的强度或峰峰值取决于脂肪或脂肪组织和软组织或瘦组织的相对吸收特性。双极声学信号的强度或峰峰值还取决于脂肪或脂肪组织与软组织或瘦组织之间的边界的温度。
示例性双极声学信号200、205和210在图2中被示出。双极声学信号200、205和210是响应于组织的关注区域ROI的热声成像而生成的,该组织的关注区域ROI包括由边界215分开的第一组织220和不同类型的第二组织225。虚线230指示与边界215对应的时间点。双极声学信号200、205和210的峰峰值的差异表示第一组织220在收缩之前膨胀到边界215和第二组织225中的程度。随着边界215的温度增加,第一组织220膨胀进入边界215和进入第二组织225的量增加。因此,每个双极声学信号200、205和210的峰峰幅度与在边界215处的温度成比例。可以看出,双极声学信号200的峰峰值大于双极声学信号205、210的峰峰值,并且双极声学信号205的峰峰值大于双极声学信号210的峰峰值。因此,当生成双极声学信号200时边界215的温度大于当生成双极信号205时边界215的温度。类似地,当生成双极声学信号205时边界215的温度大于当生成双极信号210时边界215的温度。
示例性电场强度衰减曲线300和305在图3中被示出。每条电场强度衰减曲线300、305分别表示组织310、315随着距热声成像系统26的RF施加器36的距离变化而变化的电场强度衰减。相比于与电场强度衰减曲线305相关联的组织315,与电场强度衰减曲线300相关联的组织310具有不同的温度和Grüneisen参数。将要理解的是,温度对热声双极信号具有正面或负面影响,这取决于相应类型的组织的电场强度和Grüneisen参数之间的温度关系。
成像系统20利用热声双极信号和温度之间的关系来估计在热声成像期间的组织温度。
现在转向图4,示出了在热声成像期间确定组织温度的方法。在该方法期间,最初,定位要被成像的受试者S体内的关注区域ROI,该关注区域包含由至少一个边界分隔的关注对象和参考物(步骤410)。在该实施例中,使用超声成像系统24定位关注区域ROI。具体地,由超声成像系统24获得的超声图像数据被传送到计算设备22。超声图像数据由计算设备22处理并且重构建的超声图像被呈现在显示设备上。操作者在受试者的身体上移动超声换能器28,直到关注区域被定位为止。当定位关注区域时,计算设备22将与超声换能器28的一个或更多个换能器阵列30的中心线的角度相关联的信息叠加(overlay)在显示设备上的重构建的超声图像之上。该信息被用于向操作者提供反馈以确保超声换能器28的轴向轴大体上垂直于在关注对象和参考物之间的边界。
图5示出包含关注对象505和参考物510的示例性关注区域500。在该实施例中,关注对象505是受试者的肝脏并且参考物510是受试者的肾脏。在图5中还示出了RF施加器36和热声换能器32。
然后在重构建的超声图像中标识在关注对象和参考物之间的至少一个边界(步骤420)。在该实施例中,至少一个边界是由操作者使用耦合到计算设备22的输入设备(诸如鼠标)来标识的。具体地,操作者绘制(draw)包围关注对象505的至少一部分、参考物510的至少一部分以及在关注对象和参考物的部分之间所标识的边界的框。计算设备22经由显示设备向操作者提供反馈以指示在框与边界之间的大致角度以确保框大体上垂直于边界。
图5中示出了示例性框515。可以看出,框515包围关注对象505(肝脏)的一部分、参考物510(肾脏))的一部分以及在关注对象505和参考物510之间的边界520。在肝脏和肾脏彼此密切相关的特定位置525处选择边界520。
RF施加器36被调节以生成短的RF能量脉冲。RF能量脉冲被引导到关注区域500中以将能量传递到关注区域ROI内的关注对象505和参考物510。作为响应,生成由热声换能器32检测的双极声学信号(步骤430)。
由于在超声换能器28的一个或更多个换能器阵列30的中心线与热声换能器32的一个或更多个换能器阵列34的中心线之间的角度α是已知的,操作者能够调节热声换能器32相对于受试者的身体的位置,使得热声成像系统26能够以期望的成像角σ获得关注区域的热声图像数据。期望的成像角σ使得热声换能器32的一个或更多个换能器阵列34的中心线延伸通过在关注对象505和参考物510之间的边界520。
双极声学信号进而被传送到计算设备22以进行处理(步骤440)。在该实施例中,计算设备22被编程为确定双极声学信号的峰峰幅度。
边界520的温度使用双极声学信号的峰峰幅度确定(步骤450)。在本实施例中,边界520的温度基于以下确定。
声压波的压力上升与称为Grüneisen参数的无量纲参数(dimensionlessparameter)和局部注量(local fluence)成比例。局部注量取决于组织参数,诸如,例如吸收系数、散射系数和各向异性因子。局部注量不随温度显著变化。
Grüneisen参数取决于等温压缩性、体积膨胀的热系数、质量密度和组织的恒定体积下的比热容。如图8中所示,Grüneisen参数随温度显著变化。因此,根据声压波生成的双极声学信号的强度或峰峰值取决于组织温度。
使用Grüneisen参数和双极声学信号的强度或峰峰值,可以确定边界的温度。在该实施例中,使用以下方程式计算边界的温度。
由热源H(r,t)产生的热声压力p(r,t)遵循以下方程式:
其中,r是空间位置向量,β是等压体积膨胀系数,c是声速,以及Cp是比热容。求解关于声压波p(r,t)的方程式1,产生以下前向问题:
热源H(r,t)被建模为两个因子的乘积,该两个因子是能量吸收的空间分布A(r)和时间辐射函数I(t)。能量吸收的空间分布A(r)是基于被成像的组织的特性确定的。由于热声换能器阵列30具有有限带宽,因此接收到的热声数据pd(r,t)是声压波p(r,t)和热声换能器阵列30的脉冲响应h(t)的卷积结果,如在方程式3中所示:
pd(r,t)=p(r,t)*th(t) (3)
其中*t表示一维时间卷积。
将要理解的是,对于传统热声成像,目标是通过反转该前向问题来恢复空间吸收分布A(r)。因此,辐照函数被建模为在给定时间点在整个场中均匀的时间函数。
由于用于接收热声数据的热声换能器阵列的带宽有限,准确地恢复吸收分布并非易事。因此,从热声数据中提取诸如组织温度等定量信息需要超越传统重构建方法的复杂方法。
当用RF能量脉冲加热关注区域时,每单位体积的功率沉积A(r)被表达为:
A(r)=ωε0ε″rE2(r) (4)
其中,ω是弧度频率,ε0是真空介电常数,ε″r是组织的相对介电常数的虚部,以及E(r)是电场强度。从组织获得的热声数据S(r)的强度是沉积的能量和组织的Grüneisen参数Γ的乘积:
S(r)=ΓA(r)=Γωε0ε″rE2(r) (5)
在介电损耗介质中,电场强度在其通过介质传播时衰减。衰减量由各种因子确定,诸如由例如关注区域的特性和RF施加器36的特性确定。电场的空间分布是:
E(r)=E0EA(r) (6)
其中,E0是前向到关注区域的最大电场强度,以及EA(r)是给定空间上电场的衰减。对于简单的1D情况,衰减EA(r)可以用指数形式表达:
EA(d)=e-ηd (7)
其中,η是关注区域的电场吸收系数,以及d是关注区域与施加器的距离。
如上所述,双极信号的强度或峰峰幅度表示在关注对象和参考物之间的吸收特性差异。此外,在边界处的热声数据的相位指示哪个组织(关注对象或参考物)具有更高或更低的吸收系数。在边界位置rl处测量的热声信号的强度或峰峰幅度Sl在方程式8.0中进行表达:
Sl=μ(Γ1ε″r,12ε″r,2)ωε0El 2 (8.0)
其中,下标1和2表示位于边界rl的每一侧的两个不同组织,以及El表示在边界处的入射电场强度,以及μ是校准标度常数。给定系统的校准标度常数将取决于各种因素,诸如换能器特性和施加的信号处理技术。不失一般性,对于以下推导,校准标度常数的值可以被假设为1。
如方程式8.0中所示,所获取的热声数据的强度是根据组织特性和电场强度确定的。
由于参考物的特性是已知的,为了确定关注对象的温度,仅需要在边界处的电场强度。换句话说,由于温度较高的组织与温度较低的组织相比具有不同的介电特性和热特性,因此可以推导出在参考物和关注对象之间的边界处的温度。
组织的Grüneisen参数Γ是温度的函数。在该实施例中,假设在组织的Grüneisen参数Γ与组织的温度之间的关系是线性的:
Γ(T)=aT+b (8.1)
其中,a是线性关系的斜率,T是应用于被建模的温度范围的线性部分的温度偏移(例如对应于20至80摄氏度的线性温度范围),而b是在T=0时的值。要理解的是,在其它实施例中可以使用高阶函数或其它类型的数学关系。
使用方程式8.1的线性关系,方程式8.0被重写为:
Sl=((a1T+b1)ε″r,1-(a2T+b2)ε″r,2)ωε0El 2 (8.2)
求解温度T生成方程式8.3:
T=(((Sl/(ωε0El 2))-b1ε″r,1-b2ε″r,2)/(a1ε″r,1+a2ε″r,2) (8.3)
将要理解的是,组织的介电特性也是温度相关的。对于每种组织类型,在介电特性和温度之间的关系是唯一的。上述方程式假设组织的介电特性是固定的。为了更准确地确定温度,每种组织类型的介电常数被建模为温度的函数。因此,方程式(8.2)被重写为:
Sl=(Γ1(T)ε″r,1(T)-Γ2(T)ε″r,2(T))ωε0El 2 (8.4)
Sl=((a1T+b1)ε″r,1(T)-(a2T+b2)ε″r,2(T))ωε0El 2 (8.5)
通过重新排列项,以下温度函数被定义:
g(T)≡((a1T+b1)ε″r,1(T)-(a2T+b2)ε″r,2(T))=Sl/(ωε0El 2) (8.6)
因此,对于不同的温度,函数g(T)根据模型和实验制表列出(tabulated)。
将要理解的是,在边界处的入射电场El可以使用方程式(6)来确定:
其中,E0是参考物开始时的电场强度,η参考物是参考物的衰减系数,d参考物是参考物的厚度。在该实施例中,电场强度E0通过有限差分时域(FDTD)方法建模并且基于在皮肤处的测量结果来推断出。
组织的电场衰减也取决于温度。因此,方程式9.0可以被重写为温度关系,如下所示:
使用上述方程式,以下温度函数被定义:
对于不同的温度,函数h(T)可以根据模型和实验制表列出。h(T)的每个分量都针对不同的温度进行测量或建模。使用函数h(T)确定在边界处的温度。
假设关注对象和参考物在边界处具有相同的温度。因此,一旦确定了边界520的温度,关注对象和参考物的温度就被估计出。可以重复步骤430至450以在医疗程序期间(诸如例如在热疗期间)连续监测关注对象的温度。
在另一实施例中,边界、关注对象和参考物的温度可以被用于估计第二关注对象的温度(在单独的实施例中,第二关注对象可对应于整个关注区域或受医疗程序影响的区域)。示例性方法600在图6中被示出。该方法跟随上述步骤450,其中确定边界、关注对象和参考物的温度(步骤610)。使用超声成像系统24定位第二关注对象(步骤620)。在该实施例中,第二关注对象是肿瘤。RF能量脉冲被引导以将能量传递到第二关注对象。作为响应,生成由热声换能器32检测的双极声学信号(步骤630)。双极声学信号被传送到计算设备22以进行处理(步骤640)。在该实施例中,计算设备22被编程为确定双极声学信号的峰峰幅度。
通过比较在步骤440期间确定的双极声学信号的峰峰幅度、在步骤450期间确定的边界、关注对象和参考物的温度、以及在步骤640期间确定的双极声学信号的峰峰幅度,可以估计出第二关注对象的温度(步骤650)。可以重复步骤610至650以在医疗程序期间(诸如,例如在热疗期间)连续监测第二关注对象的温度。
虽然在实施例中,关注对象被描述为肝脏并且参考物被描述为肾脏,但是本领域技术人员将理解,可以获得身体的其它部分的热声数据。如图7中所示,可以使用上述系统和方法成像的身体的各个部位包括外膜/心包脂肪组织1301、肝脏1302、皮下脂肪组织1303、内脏脂肪组织1304、皮下臀股脂肪组织1305、血管周围脂肪组织1306、心肌脂肪1307、胰腺脂肪1308、肾窦脂肪1309和肌肉脂肪1310。
将要理解的是,在一些实施例中,可以根据热声数据调节来校正热声数据。例如,热声信号以声压波的形式在空间中传播。在超声换能器阵列处接收的信号可以根据方程式10进行表达:
ps(t)=∫sp(r,t)dS (10)
其中,S是超声换能器阵列的表面积。热声换能器阵列的特性及其相对于关注区域的定位二者都可能对热声数据产生影响。
尽管在上述实施例中,参考物被描述为由操作者选择,但是本领域技术人员将理解替代方案是可用的。例如,在另一个实施例中,可以使用计算设备基于在关注区域内特定类型的组织的已知几何形状和/或已知超声特性执行的算法来自动定义参考物。此外,在参考物和关注对象之间的边界可以使用基于超声分割或热声数据分析的算法来自动定义。将要理解的是,操作者定义的方法和自动方法二者可以结合。
尽管在上述实施例中使用了双极信号的强度或峰峰幅度,但是本领域技术人员将理解可以使用其它度量,诸如例如简单峰(最大值)、p-范数、双极信号下的面积等。
将要理解的是,可以使用存储在存储器中的图像实时或离线地执行上述实施例。
尽管热声成像系统被描述为包括被配置为生成短脉冲的RF电磁辐射的RF源,但是本领域技术人员将理解,在其它实施例中,热声成像系统可以包括具有在400nm到10pm之间的波长和在10皮秒到10微秒之间的脉冲持续时间的可见光源或红外辐射源。
尽管在上述实施例中,热声成像系统和超声成像系统被描述为使用一个或更多个超声换能器阵列,但是本领域技术人员将理解,替代方案是可用的。例如,可以使用单个换能器元件、具有线性或弯曲的一维阵列的超声换能器阵列或二维超声换能器阵列。此外,凝胶状材料或水囊可以被用于将一个或更多个超声换能器阵列与关注区域对接。
尽管在上述实施例中,使用从单个关注区域获得的热声数据来估计温度,但是本领域技术人员将理解,可以分析和组合多个关注区域。
尽管在上述实施例中,血管被描述为由操作者手动标识,但是本领域技术人员将理解,可以以其它方式标识血管。例如,在另一实施例中,可以使用自动或半自动算法来标识一个或更多个血管。在其它实施例中,多普勒成像方法可以被用于标识血管。
本领域技术人员将理解,上述超声图像数据和热声数据可以是一维、二维或三维的。在实施例中,超声图像数据可以处于与热声数据不同的维度。例如,超声图像数据可以是二维的,而热声数据可以是一维的。此外,可以使用不同的视场。
在另一个实施例中,不同类型或模型的换能器阵列可以与热声和超声成像系统一起使用。在该实施例中,可以使用变换来将热声吸收图像映射到超声图像。在另一个实施例中,在换能器阵列几何形状的知识不容易获得的情况下,可以使用模型参考点(phantomreference points)将热声吸收图像映射到超声图像。在该实施例中,可以使用变换将来自热声吸收图像的已知模型参考点映射到超声图像上的模型参考点。
尽管超声成像系统被描述为使用B模式超声成像技术,但也可以使用其它技术,诸如例如功率多普勒图像、连续波多普勒图像等。
本领域技术人员将理解,可以评估其它关注对象并且可以使用其它参考物,诸如,例如心脏、肾脏、肺、食道、胸腺、乳房、前列腺、脑、肌肉、神经组织、上皮组织、膀胱、胆囊、肠、肝、胰腺、脾、胃、睾丸、卵巢、子宫、皮肤和脂肪组织。
尽管在上述实施例中,获得的是关注区域的热声数据,但是本领域技术人员将理解,可以获得大于关注区域的区域的热声数据。
使用前述说明,本发明可以通过使用标准编程和/或工程技术来产生编程软件、固件、硬件或它们的任何组合而被实现为机器、过程或制品。
具有计算机可读指令的任何结果程序可以被存储在一个或更多个计算机可用介质中,诸如存储器设备或传输设备,从而制造根据本发明的计算机程序产品或制品。因此,功能可以作为计算机程序被赋予在物理设备上,该计算机程序作为由处理器执行的指令存在于任何计算机可读介质上(诸如在任何存储器设备上或在任何传输设备中)。
存储器设备的示例包括硬盘驱动器、软盘、光盘、磁带、诸如闪存、RAM、ROM、PROMS的半导体存储器等。网络的示例包括但不限于互联网、内联网、基于电话/调制解调器的网络通信、硬接线/有线通信网络、蜂窝通信、无线电波通信、卫星通信以及其它固定或移动网络系统/通信链路。
体现本发明的机器可以涉及一个或更多个处理系统,包括例如计算机处理单元(CPU)或处理器、存储器/存储设备、通信链路、通信/传输设备、服务器、I/O设备,或一个或更多个处理系统的任何子部件或单独部分,包括软件、固件、硬件或它们的任何组合或子组合,它们体现了权利要求中阐述的本发明。
使用本文提供的描述,本领域技术人员将能够容易地将如所述创建的软件与适当的或专用计算机硬件组合以创建体现本发明的计算机系统和/或计算机子组件,并且创建用于执行本发明方法的计算机系统和/或计算机子部件。
尽管上面已经参考附图描述了实施例,但是本领域技术人员将理解,可以在不脱离由所附权利要求限定的范围的情况下做出变化和修改。

Claims (12)

1.一种用于估计组织温度的方法,所述方法包括:
(i)使用射频(RF)施加器将一个或更多个RF能量脉冲引导到组织的关注区域中,所述组织的关注区域包括由至少一个边界隔开的关注对象和至少一个参考物;其中,所述至少一个边界位于至少两个不同类型的组织之间的位置处,进一步地其中,所述两个不同类型的组织是所述关注对象和所述至少一个参考物;
(ii)使用声学接收器检测响应于所述RF能量脉冲而在所述关注区域中生成的至少一个双极声学信号,并处理所述至少一个双极声学信号以确定所述至少一个双极声学信号的峰峰幅度;
(iii)使用所述至少一个双极声学信号的所述峰峰幅度计算在所述至少一个边界处的温度;以及
(iv)基于所计算的在所述至少一个边界处的温度来估计所述关注对象的温度。
2.根据权利要求1所述的方法,还包括使用所述声学接收器检测由第二关注对象响应于所述RF能量脉冲而生成的至少一个另外的双极声学信号,并且处理所述至少一个另外的双极声学信号以确定所述至少一个另外的双极声学信号的峰峰幅度。
3.根据权利要求2所述的方法,包括比较所述双极声学信号的峰峰幅度、所计算的在所述至少一个边界处的温度和所述另外的双极声学信号的峰峰幅度,以估计所述第二关注对象的温度。
4.根据权利要求1所述的方法,其中,所述两个不同类型的组织是以下之一:
肌肉和脂肪;
血管和脂肪;以及
肝组织和肾组织。
5.根据权利要求1所述的方法,还包括在治疗期间重复步骤(i)至(iii)以连续监测组织温度。
6.一种用于确定在关注区域内的组织温度的系统,所述关注区域包括由至少一个边界隔开的关注对象和参考物;其中,所述至少一个边界位于至少两个不同类型的组织之间的位置处,进一步地其中,所述两个不同类型的组织是所述关注对象和所述参考物,所述系统包括:
热声成像系统,所述热声成像系统包括:可调节射频(RF)施加器,所述可调节射频(RF)施加器被配置为将RF能量脉冲发射到组织的关注区域中并加热在所述组织的关注区域中的组织;以及声学接收器,所述声学接收器被配置为接收响应于在所述关注区域中的组织加热而生成的双极声学信号;以及
一个或更多个处理器,所述一个或更多个处理器被配置为:
处理接收到的响应于所述RF能量脉冲而在所述关注区域中生成的至少一个双极声学信号以确定所述至少一个双极声学信号的峰峰幅度;
使用所述至少一个双极声学信号的所述峰峰幅度计算在所述至少一个边界处的温度;以及
基于所计算的在所述至少一个边界处的温度来估计所述关注对象的温度。
7.根据权利要求6所述的系统,其中,所述一个或更多个处理器还被配置为:
使用所述声学接收器,检测由第二关注对象响应于所述RF能量脉冲而生成的至少一个另外的双极声学信号,并处理所述至少一个另外的双极声学信号以确定所述至少一个另外的双极声学信号的峰峰幅度。
8.根据权利要求7所述的系统,其中,所述一个或更多个处理器还被配置为:
比较所述双极声学信号的峰峰幅度、所计算的在所述至少一个边界处的温度以及所述另外的双极声学信号的峰峰幅度,以估计所述第二关注对象的温度。
9.一种监测组织温度的方法,所述方法包括:
(i)使用射频(RF)施加器将一个或更多个RF能量脉冲引导到组织的关注区域中,所述组织的关注区域包括由至少一个边界隔开的关注对象和至少一个参考物;其中,所述至少一个边界位于至少两个不同类型的组织之间的位置处,进一步地其中,所述两个不同类型的组织是所述关注对象和所述至少一个参考物;
(ii)使用声学接收器,检测响应于所述RF能量脉冲而在所述关注区域中生成的至少一个双极声学信号,并处理所述至少一个双极声学信号以确定所述至少一个双极声学信号的峰峰幅度;
(iii)使用所述至少一个双极声学信号的峰峰幅度计算在所述至少一个边界处的温度;以及
(iv)重复所述引导、检测和计算,以监测在所述至少一个边界处的温度。
10.根据权利要求9所述的方法,还包括:
基于所计算的在所述至少一个边界处的温度估计所述关注对象的温度。
11.根据权利要求9所述的方法,还包括:
使用所述声学接收器,检测由第二关注对象响应于所述RF能量脉冲而生成的至少一个另外的双极声学信号,并处理所述至少一个另外的双极声学信号以确定所述至少一个另外的双极声学信号的峰峰幅度。
12.根据权利要求11所述的方法,包括比较所述双极声学信号的峰峰幅度、所计算的在所述至少一个边界处的温度和所述另外的双极声学信号的峰峰幅度,以估计所述第二关注对象的温度。
CN201980086822.7A 2018-12-27 2019-12-04 用于监测组织温度的方法和系统 Active CN113226208B (zh)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
US16/233,935 US11141067B2 (en) 2018-12-27 2018-12-27 Method and system for monitoring tissue temperature
US16/233,935 2018-12-27
PCT/US2019/064447 WO2020139527A1 (en) 2018-12-27 2019-12-04 Method and system for monitoring tissue temperature

Publications (2)

Publication Number Publication Date
CN113226208A CN113226208A (zh) 2021-08-06
CN113226208B true CN113226208B (zh) 2024-03-26

Family

ID=71123694

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201980086822.7A Active CN113226208B (zh) 2018-12-27 2019-12-04 用于监测组织温度的方法和系统

Country Status (4)

Country Link
US (1) US11141067B2 (zh)
EP (1) EP3902490B1 (zh)
CN (1) CN113226208B (zh)
WO (1) WO2020139527A1 (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107861036A (zh) * 2017-11-10 2018-03-30 国网新疆电力有限公司培训中心 多用绝缘杆检测装置
US11067543B2 (en) * 2019-10-03 2021-07-20 Endra Life Sciences Inc. Method and system for determining a parameter of a material of interest

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105208959A (zh) * 2013-04-18 2015-12-30 皇家飞利浦有限公司 用于沉积能量的治疗系统
US9888880B1 (en) * 2017-08-01 2018-02-13 Endra Life Sciences Inc. Method and system for estimating fractional fat content of an object

Family Cites Families (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7048756B2 (en) * 2002-01-18 2006-05-23 Apasara Medical Corporation System, method and apparatus for evaluating tissue temperature
WO2012088149A2 (en) * 2010-12-20 2012-06-28 Virginia Tech Intellectual Properties, Inc. High-frequency electroporation for cancer therapy
WO2017192558A1 (en) * 2016-05-02 2017-11-09 Arizona Board Of Regents On Behalf Of The University Of Arizona Thermoacoustic image-guided microwave therapy system
EP3463071A4 (en) * 2016-06-05 2019-12-25 Endra, Inc. METHOD AND SYSTEM FOR ESTIMATING THE FACTIONAL CONTENT OF AN OBJECT
WO2018200865A1 (en) 2017-04-27 2018-11-01 Epix Therapeutics, Inc. Determining nature of contact between catheter tip and tissue
US10631734B1 (en) * 2018-12-27 2020-04-28 Endra Life Sciences Inc. Method and system for monitoring tissue temperature
US11478153B2 (en) * 2018-12-27 2022-10-25 Endra Life Sciences Inc. System for monitoring tissue temperature

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105208959A (zh) * 2013-04-18 2015-12-30 皇家飞利浦有限公司 用于沉积能量的治疗系统
US9888880B1 (en) * 2017-08-01 2018-02-13 Endra Life Sciences Inc. Method and system for estimating fractional fat content of an object

Also Published As

Publication number Publication date
CN113226208A (zh) 2021-08-06
US11141067B2 (en) 2021-10-12
WO2020139527A1 (en) 2020-07-02
EP3902490A4 (en) 2022-09-14
EP3902490B1 (en) 2024-02-14
EP3902490A1 (en) 2021-11-03
US20200205668A1 (en) 2020-07-02
EP3902490C0 (en) 2024-02-14

Similar Documents

Publication Publication Date Title
EP3668382B1 (en) Method for estimating fractional fat content of an object
EP3661416B1 (en) Method and system for estimating fractional fat content of an object
CN113624360B (zh) 用于监测组织温度的方法和系统
CN113226208B (zh) 用于监测组织温度的方法和系统
CN113227782B (zh) 用于估计感兴趣对象的分数脂肪含量的方法和系统
EP4037554B1 (en) Method and system for determining a material type of an object of interest
US11478153B2 (en) System for monitoring tissue temperature
EP4021284B1 (en) Method and system for determining at least one parameter of interest of a material
EP4037551B1 (en) Method and system for determining a parameter of a material of interest
Kim et al. A novel 3D ultrasound thermometry method for HIFU ablation using an ultrasound element
US20230363738A1 (en) Ablation monitoring and control using three-dimensional echo decorrelation imaging
WO2021026478A1 (en) Method and system for identifying a material of interest

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