CN113925598A - 利用超声成像系统对组织进行热监测的方法和系统 - Google Patents

利用超声成像系统对组织进行热监测的方法和系统 Download PDF

Info

Publication number
CN113925598A
CN113925598A CN202110744846.6A CN202110744846A CN113925598A CN 113925598 A CN113925598 A CN 113925598A CN 202110744846 A CN202110744846 A CN 202110744846A CN 113925598 A CN113925598 A CN 113925598A
Authority
CN
China
Prior art keywords
signal
tissue
interest
region
roi
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.)
Pending
Application number
CN202110744846.6A
Other languages
English (en)
Inventor
J·齐格勒
A·G·伊蓝斯玛利奎茨
E·J·戈麦斯阿泰德
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.)
GE Precision Healthcare LLC
Original Assignee
GE Precision Healthcare LLC
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 GE Precision Healthcare LLC filed Critical GE Precision Healthcare LLC
Publication of CN113925598A publication Critical patent/CN113925598A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00—Image analysis
    • G06T7/0002—Inspection of images, e.g. flaw detection
    • G06T7/0012—Biomedical image inspection
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B18/00—Surgical instruments, devices or methods for transferring non-mechanical forms of energy to or from the body
    • A61B18/04—Surgical instruments, devices or methods for transferring non-mechanical forms of energy to or from the body by heating
    • A61B18/12—Surgical instruments, devices or methods for transferring non-mechanical forms of energy to or from the body by heating by passing a current through the tissue to be heated, e.g. high-frequency current
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/01—Measuring temperature of body parts ; Diagnostic temperature sensing, e.g. for malignant or inflamed tissue
    • A61B5/015—By temperature mapping of body part
    • 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
    • A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/46—Ultrasonic, sonic or infrasonic diagnostic devices with special arrangements for interfacing with the operator or the patient
    • A61B8/467—Ultrasonic, sonic or infrasonic diagnostic devices with special arrangements for interfacing with the operator or the patient characterised by special input means
    • A61B8/469—Ultrasonic, sonic or infrasonic diagnostic devices with special arrangements for interfacing with the operator or the patient characterised by special input means for selection of a region of interest
    • 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/5215—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data
    • 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/5215—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data
    • A61B8/5223—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data for extracting a diagnostic or physiological parameter from medical diagnostic data
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00—Image analysis
    • G06T7/0002—Inspection of images, e.g. flaw detection
    • G06T7/0012—Biomedical image inspection
    • G06T7/0014—Biomedical image inspection using an image reference approach
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V10/00—Arrangements for image or video recognition or understanding
    • G06V10/20—Image preprocessing
    • G06V10/25—Determination of region of interest [ROI] or a volume of interest [VOI]
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V10/00—Arrangements for image or video recognition or understanding
    • G06V10/40—Extraction of image or video features
    • G06V10/42—Global feature extraction by analysis of the whole pattern, e.g. using frequency domain transformations or autocorrelation
    • G06V10/431—Frequency domain transformation; Autocorrelation
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V10/00—Arrangements for image or video recognition or understanding
    • G06V10/40—Extraction of image or video features
    • G06V10/52—Scale-space analysis, e.g. wavelet analysis
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B18/00—Surgical instruments, devices or methods for transferring non-mechanical forms of energy to or from the body
    • A61B2018/00571—Surgical instruments, devices or methods for transferring non-mechanical forms of energy to or from the body for achieving a particular surgical effect
    • A61B2018/00577—Ablation
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B18/00—Surgical instruments, devices or methods for transferring non-mechanical forms of energy to or from the body
    • A61B2018/00636—Sensing and controlling the application of energy
    • A61B2018/00773—Sensed parameters
    • A61B2018/00791—Temperature
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2218/00—Aspects of pattern recognition specially adapted for signal processing
    • G06F2218/12—Classification; Matching
    • G06F2218/16—Classification; Matching by matching signal segments
    • G06F2218/20—Classification; Matching by matching signal segments by applying autoregressive analysis
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/10—Image acquisition modality
    • G06T2207/10132—Ultrasound image
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/20—Special algorithmic details
    • G06T2207/20048—Transform domain processing
    • G06T2207/20064—Wavelet transform [DWT]
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/20—Special algorithmic details
    • G06T2207/20092—Interactive image processing based on input by user
    • G06T2207/20104—Interactive definition of region of interest [ROI]
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/30—Subject of image; Context of image processing
    • G06T2207/30004—Biomedical image processing
    • G06T2207/30096—Tumor; Lesion
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V2201/00—Indexing scheme relating to image or video recognition or understanding
    • G06V2201/03—Recognition of patterns in medical or anatomical images
    • G06V2201/031—Recognition of patterns in medical or anatomical images of internal organs

Abstract

本发明提供了利用超声成像系统对组织进行热监测的各种方法和系统。在一个实施方案中，方法包括经由超声探头采集超声图像，选择超声图像中的至少一个感兴趣区域，从至少一个感兴趣区域提取特征，基于特征对至少一个感兴趣区域中的组织的热状态进行分类，以及经由显示装置输出组织的热状态。这样，超声成像可用于消融规程期间组织的实时热监测，从而提高消融规程的准确性和功效。

Description

利用超声成像系统对组织进行热监测的方法和系统
技术领域
本文公开的主题的实施方案涉及超声成像，并且更具体地涉及经由超声成像的热监测。
背景技术
癌症是一组涉及异常细胞生长的疾病，有可能侵入或者扩散到身体的其他部位。癌症通常表现为肿瘤形式的组织异常生长，该肿瘤可能局限于患者身体的特定区域(例如，与特定身体部位或器官相关联)或者可能扩散到全身。良性和恶性肿瘤通常均经由外科介入来治疗和移除，因为通常外科手术提供完全移除和治愈的最大机会，尤其是在癌症尚未扩散到身体的其他部位的情况下。
在一些情况下，对于此种侵入式规程，肿瘤的位置可能是不可触及的。在此类情况下，微创热疗法是固体内部恶性肿瘤的潜在疗法。例如，射频消融(RFA)可用于通过经皮施加直接施加到肿瘤的高频能量来通过热破坏癌或病理组织。更具体地，通过在组织内产生交变电场，在邻近电场源的区域(通常为经由RFA针定位在该区域中的电极)中搅动离子，从而增加该区域内的温度，导致病理组织的破坏或消融。RFA是针对肝细胞癌(HCC，成人中最常见类型的原发性肝癌之一)的主要治疗，尤其是当HCC不可切除时。
发明内容
在一个实施方案中，方法包括经由超声探头采集超声图像，选择超声图像中的至少一个感兴趣区域，从至少一个感兴趣区域提取特征，基于特征对至少一个感兴趣区域中的组织的热状态进行分类，以及经由显示装置输出组织的热状态。这样，超声成像可用于消融规程期间组织的实时热监测，从而提高消融规程的准确性和功效。
应当理解，提供上面的简要描述来以简化的形式介绍在具体实施方式中进一步描述的精选概念。这并不意味着识别所要求保护的主题的关键或必要特征，该主题的范围由具体实施方式后的权利要求书唯一地限定。此外，所要求保护的主题不限于解决上文或本公开的任何部分中提到的任何缺点的实施方式。
附图说明
通过参考附图阅读以下对非限制性实施方案的描述将更好地理解本发明，其中以下：
图1示出了根据一个实施方案的示出示例性超声成像系统的框图；
图2示出了根据一个实施方案的示出用于在超声成像期间进行热监测的示例性方法的高级流程图；
图3示出了根据一个实施方案的示出用于从超声图像中的感兴趣区域提取纹理特征的高级流程图；
图4示出了根据一个实施方案的示出利用超声成像系统进行热监测的示例性方法的框图；
图5示出了根据一个实施方案的示出用于从超声图像中的感兴趣区域提取纹理特征的示例性方法的框图；
图6示出了根据一个实施方案的描绘不同热状态下组织的超声图像的感兴趣区域的一组图像；
图7示出了根据一个实施方案的描绘用于图像到信号转换的示例性像素轨线的一组示意图；
图8示出了根据一个实施方案的示出第一热状态的具有不同像素轨线的单个感兴趣区域构造的示例性信号的一组曲线图；
图9示出了根据一个实施方案的示出第二热状态的具有不同像素轨线的单个感兴趣区域构造的示例性信号的一组曲线图；
图10示出了根据一个实施方案的示出信号偏置校正的一组曲线图；
图11示出了根据一个实施方案的示出第一热状态的纹理信号到伪频带的示例性分解的一组曲线图；
图12示出了根据一个实施方案的示出第二热状态的纹理信号到伪频带的示例性分解的一组曲线图；
图13示出了根据一个实施方案的示出通过分解纹理信号的自回归建模获得的示例性功率谱密度估计的一组曲线图；
图14示出了根据一个实施方案的针对多个热状态计算的一组示例性相关系数标测图；
图15示出了根据一个实施方案的示出不同热状态的经排序的相关系数的示例性相异性度量的曲线图；
图16示出了根据一个实施方案的示出不同热状态的从相异性度量中提取的第一示例性特征的一组曲线图；
图17示出了根据一个实施方案的示出不同热状态的从相异性度量中提取的第二示例性特征的一组曲线图；以及
图18示出了根据一个实施方案的示出不同热状态的从经排序的相关系数的相异性度量中提取的第三示例性特征和第四示例性特征的一组曲线图。
具体实施方式
以下描述涉及超声成像的各种实施方案。具体地，提供了利用超声成像系统对组织进行热监测的方法和系统。在图1中示出了根据本发明技术可用于获取图像和执行热监测的超声成像系统的示例。如上文所讨论，射频消融(RFA)是主要应用于HCC的治疗方法。在此类热规程期间，观察消融进展(即，监测和确定组织的温度范围)是主要问题。达到高于60℃的温度确保病理组织中的细胞被成功并且不可逆地破坏。因此，热监测用于指示某个区段中是否已达到消融温度。目前，磁共振成像(MRI)是在热消融期间用于实时和非侵入性温度监测的黄金标准。然而，使用MRI进行热监测的主要缺点在于需要专用的MRI兼容消融系统、设备和MRI时间的成本高、MRI孔中的空间有限以及可用性低。相比之下，超声成像系统(诸如图1所示的超声成像系统)是便携式的低成本成像模态，并且可广泛获得。此种系统可用于例如RFA中的针引导，但是先前的系统不能在RFA规程期间测量和显示组织温度。利用超声成像系统进行热监测的方法，诸如图2所示的方法，包括从超声图像中的感兴趣区域(ROI)提取特征，以及基于所提取的特征来评估组织的状态。例如，可训练一个或多个机器学习算法以基于所提取的特征对ROI中所描绘的组织的温度进行分类。一种用于提取特征的方法，诸如图3所示的方法，包括将ROI从二维图像转换成一维信号，将每个信号分解成多个频带，用自回归(AR)模型生成每个分解信号的估计值，测量相同频带中估计值的相关性，基于相关性确定相异性的量度，以及从这些相异性的量度中提取特征。如图4和图5所示，用于提取特征的方法，包括将ROI的改变纹理(例如，由于组织的受热)与ROI的初始纹理(例如，消融规程开始时组织的状态)进行比较，例如，如图6所示。因此，特征包括由热引起的图像纹理之间的相异性的定量量度，其可因此与组织的相对温度升高相关联(例如，通过一个或多个机器学习算法的有效训练)。与直接评估组织的图像纹理相比，此类特征可为分类组织的热状态提供更高的精度。为了从组织的图像导出特征，可通过使用相应的多个像素遍历轨线将图像转换成多个信号，如图7所示。初始纹理(如图8所示)以及改变纹理(如图9所示)的所得信号尽管从同一图像获得，但仍表现出不同的信号动态。此外，在执行如图10所示的偏置校正之后，每个信号可分解成如图11和图12所示的多个伪频带。使用基于频率的信号的自回归建模产生功率谱密度(PSD)的AR估计值，例如，如图13所示，其作为信号相对于频率的功率含量的量度，提供了对组织纹理的稳健表征。此外，本文提供的方法包括计算初始纹理和改变纹理的此类AR估计值之间的相关系数，从而产生由热引起的组织纹理的相对变化的稳健量度，如图14所示。然后将这些量度组织成相异性度量，如图15所示，从而提供了用于提取上文提及的特征的动态纹理偏移的合适表示。具体地，如图16至图18所示，特征可来源于相异性度量，以提供组织纹理之间随时间推移的相异性的相对简单的定量量度，这继而允许在采集超声图像时对组织的热状态进行稳健分类。这样，超声成像系统可通过处理如本文所述的超声图像来快速并且自动地对组织的温度进行分类。此外，为了简化分类，超声成像系统可相对于临界消融温度对组织纹理的状态进行分类。此外，考虑到不同患者中的不同组织可在不同温度下实现消融，本文提供的用于对组织的热状态进行分类的方法通常可区分对组织的可逆和不可逆损伤，从而针对患者而不是在预定的临界消融温度下对组织纹理的分类进行个性化。这样，执行消融规程的人可使用超声成像系统对组织进行稳健、实时的热监测，以确保消融规程的准确性和功效。
图1示出了根据本公开的一个实施方案的超声成像系统100的示意图。超声成像系统100包括发射波束形成器101和发射器102，该发射器驱动换能器阵列(本文中称为探头106)内的元件(例如，换能器元件)104，以将脉冲超声信号(本文中称为发射脉冲)发射到身体(未示出)中。根据一个实施方案，探头106可以是一维换能器阵列探头。然而，在一些实施方案中，探头106可以是二维矩阵换能器阵列探头。如以下进一步解释的，换能器元件104可以由压电材料构成。当向压电晶体施加电压时，晶体物理地膨胀和收缩，从而发射超声球波。这样，换能器元件104可将电子发射信号转换为声学发射波束。
在探头106的元件104将脉冲超声信号发射到(患者的)身体中之后，脉冲超声信号从身体内部的结构(如血细胞或肌肉组织)向后散射，以产生返回到元件104的回波。回波被元件104转换成电信号或超声数据，并且电信号被接收器108接收。表示所接收的回波的电信号穿过输出超声数据的接收波束形成器110。另外，换能器元件104可根据所接收的回波产生一个或多个超声脉冲以形成一个或多个发射波束。
根据一些实施方案，探头106可包含电子电路来执行发射波束形成和/或接收波束形成的全部或部分。例如，发射波束形成器101、发射器102、接收器108和接收波束形成器110的全部或一部分可位于探头106内。在本公开中，术语“扫描”或“扫描中”可也用于指通过传输和接收超声信号的过程来采集数据。本公开中，术语“数据”可以用于指称用超声成像系统来采集的一个或多个数据集。在一个实施方案中，可使用经由超声系统100获取的数据来训练机器学习模型。用户界面115可用于控制超声成像系统100的操作，包括用于控制患者数据(例如，患者病史)的输入、用于改变扫描或显示参数、用于启动探头复极化序列等。用户界面115可包括以下项中的一者或多者：旋转元件、鼠标、键盘、轨迹球、链接到特定动作的硬键、可被配置为控制不同功能的软键，以及显示在显示设备118上的图形用户界面。
超声成像系统100还包括处理器116，该处理器用以控制发射波束形成器101、发射器102、接收器108和接收波束形成器110。处理器116与探头106进行电子通信(例如，通信地连接)。出于本公开的目的，术语“电子通信”可被定义为包括有线通信和无线通信两者。处理器116可根据存储在处理器116的存储器和/或存储器120上的指令来控制探头106采集数据。处理器116控制元件104中的哪些是活动的以及从探头106发射的波束的形状。处理器116还与显示设备118进行电子通信，并且处理器116可将数据(例如，超声数据)处理成图像以用于在显示设备118上显示。处理器116可以包括根据一个实施方案的中央处理器(CPU)。根据其他实施方案，处理器116可包括能够执行处理功能的其他电子部件，诸如数字信号处理器、现场可编程门阵列(FPGA)或图形板。根据其他实施方案，处理器116可包括能够执行处理功能的多个电子部件。例如，处理器116可包括从电子部件的列表中选择的两个或更多个电子部件，这些电子部件包括：中央处理器、数字信号处理器、现场可编程门阵列和图形板。根据另一个实施方案，处理器116还可包括解调RF数据并且生成原始数据的复合解调器(未示出)。在另一个实施方案中，解调可以在处理链中较早地执行。处理器116适于根据数据上的多个可选超声模态来执行一个或多个处理操作。在一个示例中，可在扫描会话期间实时处理数据，因为回波信号被接收器108接收并且被传输至处理器116。出于本公开的目的，术语“实时”被定义为包括在没有任何有意延迟的情况下执行的过程。例如，实施方案可以7至20帧/秒的实时速率获取图像。超声成像系统100能够以显著更快的速率获取一个或多个平面的2D数据。然而，应当理解，实时帧速率可取决于采集用于显示的每帧数据所花费的时间长度。因此，当采集相对大量的数据时，实时帧速率可能较慢。因此，一些实施方案可具有显著快于20帧/秒的实时帧速率，而其他实施方案可具有低于7帧/秒的实时帧速率。数据可在扫描会话期间临时存储在缓冲器(未示出)中，并且在实时或离线操作中以不太实时的方式处理。本发明的一些实施方案可包括多个处理器(未示出)，以处理根据上文所述的示例性实施方案的由处理器116处理的处理任务。例如，在显示图像之前，可利用第一处理器来解调和抽取RF信号，同时可使用第二处理器来进一步处理数据(例如，通过如本文进一步描述的那样扩充数据)。应当理解，其他实施方案可使用不同的处理器布置方式。
超声成像系统100可以例如10Hz至30Hz的帧速率(例如，每秒10帧至30帧)连续采集数据。根据数据生成的图像可以在显示设备118上以相似的帧速率刷新。其他实施方案能够以不同速率获取并且显示数据。例如，根据帧的大小和预期的应用，一些实施方案可以小于10Hz或大于30Hz的帧速率采集数据。包括存储器120，用于存储经处理的采集数据的帧。在示例性实施方案中，存储器120具有足够的容量来存储至少几秒钟的超声数据帧。数据帧的存储方式便于根据其采集顺序或时间进行检索。存储器120可包括任何已知的数据存储介质。
在本发明的各种实施方案中，处理器116可通过不同的模式相关模块(例如，B模式、彩色多普勒、M模式、彩色M模式、频谱多普勒、弹性成像、TVI、应变、应变速率等)来处理数据，以形成2D或3D数据。例如，一个或多个模块可生成B模式、彩色多普勒、M模式、彩色M模式、频谱多普勒、弹性成像、TVI、应变、应变速率以及它们的组合，等等。作为一个示例，一个或多个模块可处理彩色多普勒数据，其可包括传统彩色血流多普勒、功率多普勒、HD流，等等。图像线和/或帧存储在存储器中，并且可包括指示图像线和/或帧存储在存储器中的时间的定时信息。这些模块可包括例如扫描转换模块，用于执行扫描转换操作，以将所采集的图像从波束空间坐标转换为显示空间坐标。可以提供视频处理器模块，该视频处理器模块从存储器读取所采集的图像并且在对患者执行规程(例如，超声成像)时实时显示图像。视频处理器模块可包括单独的图像存储器，并且超声图像可被写入图像存储器以便由显示设备118读取和显示。
在本公开的各种实施方案中，超声成像系统100的一个或多个部件可以包括在便携手持式超声成像设备中。例如，显示设备118和用户界面115可以集成到手持式超声成像设备的外部表面中，该手持式超声成像设备可还包括处理器116和存储器120。探头106可以包括与手持式超声成像设备进行电子通信以收集原始超声数据的手持式探头。发射波束形成器101、发射器102、接收器108和接收波束形成器110可以包括在超声成像系统100的相同或不同部分中。例如，发射波束形成器101、发射器102、接收器108和接收波束形成器110可以包括在手持式超声成像设备、探头以及它们的组合中。
在执行二维超声扫描之后，生成包含扫描线及其样本的数据块。在应用后端滤波器之后，执行称为扫描转换的过程，以将二维数据块变换为具有附加扫描信息(诸如深度、每条扫描线的角度等)的可显示的位图图像。在扫描转换期间，应用内插技术来填充所得图像中的丢失的孔(即，像素)。出现这些丢失的像素是因为二维块的每个元素通常应该覆盖所得图像中的许多像素。例如，在当前超声成像系统中，应用了双三次插值，其利用了二维块的相邻元素。因此，如果与位图图像的尺寸相比，二维块相对较小，则扫描转换后的图像将包括分辨率差或分辨率低的区域，特别是对于深度较大的区域。
图2示出了根据一个实施方案的示出用于在超声成像期间进行热监测的示例性方法200的高级流程图。具体地，方法200涉及处理超声图像以确定组织随时间推移的热状态，例如在消融规程期间。参照图1的系统和部件描述方法200，但是应当理解，在不脱离本公开的范围的情况下，方法200可用其他系统和部件来实现。方法200可被实现为非暂态存储器诸如存储器120中的可执行指令，该可执行指令在由处理器诸如处理器116执行时使得处理器执行下文所述的动作。
方法200在205处开始。在205处，方法200开始患者的超声扫描会话以提供对患者组织的热监测。超声扫描会话可在射频消融(RFA)规程期间执行，以为用户在组织内定位射频消融针提供图像引导。在这种示例中，超声扫描会话可包括控制超声换能器或超声探头106，例如以在正在执行RFA规程的区域执行标准B模式成像。在其他示例中，超声扫描会话可主要涉及经由聚焦超声执行热消融。在此类示例中，超声扫描会话可包括控制超声探头106经由聚焦超声执行感兴趣消融区域(ROI)的热消融，此外还采集用于组织的热监测的B模式图像。另选地，单独的超声成像系统可用于热消融和热监测。此外，可执行超声扫描会话以在经由其他消融模式(诸如微波治疗、激光消融和其他热递送方法)消融组织期间提供组织的热监测。在其他示例中，可在没有并发消融规程的情况下执行超声扫描会话。在任何情况下，在205处发起的超声扫描会话包括采集用于组织热监测的B模式超声图像。
为此，在210处，方法200采集超声图像。超声图像可包括B模式或灰度超声图像。在215处，方法200选择超声图像中的感兴趣区域(ROI)。例如，方法200选择多个ROI，每个ROI包括分布在整个超声图像中的图像贴片。作为例示性并且非限制性示例，每个ROI可包括超声图像内的20*20像素贴片。为了包括关于组织纹理的更多信息，作为另一个例示性和非限制性示例，每个ROI可包括超声图像中的31*31像素贴片。应当理解，每个ROI的尺寸可大于或小于31*31像素，例如，取决于超声图像的分辨率和视场。
至少一个ROI可被定位成与消融ROI对准。例如，如果消融方法是RFA，则多个ROI中的至少一个ROI可定位在超声图像内的RFA针的尖端处。除此之外或另选地，多个ROI中的至少一个ROI可以消融目标对象(诸如病变)的位置为中心。多个ROI中的剩余ROI可彼此相距最小距离分布在整个超声图像中，例如以确保具有不同热变化率的多种ROI。ROI的位置可自动地、手动地或它们的组合来选择。例如，用户可经由超声成像系统100的用户界面115选择消融ROI或消融目标诸如病灶的位置。然后，方法200可在消融目标的所选择的位置处选择用于热监测的ROI，并且相对于以消融目标为中心的ROI位置自动选择附加ROI的位置。在其他示例中，用户可经由用户界面115选择用于热监测的每个ROI的位置。在其他示例中，方法200可例如基于预先确定的ROI位置，或者通过经由图像处理和/或机器学习方法识别感兴趣区域，在没有用户输入的情况下自动选择每个ROI的位置。例如，可训练卷积神经网络以利用B模式超声图像作为输入自动确定用于热监测的ROI的多个位置。这样，ROI可自动选择，例如，具有用于热监测的期望位置(例如，远离骨和组织之间的边界，例如或器官之间的边界，这可能在热监测中潜在地引入错误)。
在一些示例中，所选择的ROI可一致地应用于在超声扫描会话期间采集的每个超声图像，以在所选择的ROI内执行热监测。为此，所选择的ROI的位置可相对于超声图像内描绘的解剖结构固定或设置，使得不管超声探头106在超声扫描会话期间相对于组织轻微运动或组织本身轻微运动(例如，由于呼吸或其他患者运动)，热监测都可一致地应用于每个ROI。
在选择用于热监测的ROI之后，方法200继续到220。在220处，方法200从ROI中提取特征。具体地，方法200从ROI中提取用于热监测的特征或热特征。方法200将特征或热特征与例如被配置用于ROI热监测的机器学习一起使用。因此，每个特征包括ROI的可测量特性或特征，其可响应于ROI中组织热状态的变化而动态地改变。例如，当组织被加热到更高的温度时，如超声图像中所描绘的组织的纹理变化。因此，提取的特征量化超声图像中的组织纹理。一种用于提取特征以研究由热引起的组织纹理的动态变化的方法，诸如用于本文参照图3和图4进一步描述的ROI中提取特征的方法，该方法包括将ROI的图像数据转换成信号，以及经由基于频率的自回归模型方法从信号中提取特征。
在从每个ROI中提取热特征之后，方法200继续到225。在225处，方法200基于所提取的特征评估ROI中的组织的状态。方法200可使用所提取的特征作为一个或多个机器学习算法的输入，该一个或多个机器学习算法被配置为能够基于所提取的特征评估组织随时间推移的状态。作为例示性和非限制性示例，方法200可使用包括随机森林分类器(RFC)、支持矢量机(SVM)和人工神经网络(ANN)中的一者或多者的机器学习算法，基于所提取的特征评估ROI中组织的状态。
基于所提取的特征评估ROI中组织的状态可包括对每个ROI中组织的热状态进行分类。例如，可训练一个或多个机器学习算法以对每个ROI中组织的温度进行分类。因此，一个或多个机器学习算法可基于输入到一个或多个机器学习算法的所提取的特征输出每个ROI中组织的温度分类。例如，温度分类可包括预定义温度仓内的温度分类。例如，每个ROI中的温度可被分类为高于或低于一个或多个温度阈值。作为例示性和非限制性示例，每个ROI中的温度可被分类为例如高于或低于60℃，或者另一个温度阈值。一个或多个温度阈值可基于消融动态确定，例如，使得一个或多个阈值对应于用于消融的一个或多个临界温度，诸如组织被消融规程永久影响或损坏的温度。此种临界温度可针对例如不同类型的组织或不同类型的患者而变化，因此在一些示例中，可提供至少两个温度阈值以限定消融范围。例如，消融可发生在55℃和65℃之间，因此每个ROI中的温度可被分类为低于第一温度阈值(例如，55℃)、高于第二温度阈值(例如，高于65℃)或在其间(例如，在55℃和65℃之间)。然而，由于消融规程的长度不足以在组织达到第一温度阈值后供组织达到第二温度阈值，因此在一些示例中，可仅训练一个或多个机器学习算法以相对于第一温度阈值或较低温度阈值对温度进行分类。又如，预定义温度仓可包括用于每个温度登记的仓，使得可更精确地对组织的温度进行分类。
另选地，可训练一个或多个机器学习算法以对纹理状态的动态变化进行分类。例如，一个或多个机器学习算法可将组织的状态分类为可逆(例如，组织尚未消融)或不可逆(例如，组织消融)，并且除此之外可将组织的状态分类为凝结。这样，机器学习算法可识别组织的状态是否指示消融完成，这可针对不同组织和/或不同患者在不同温度下发生。
在评估每个ROI中组织的状态之后，方法200继续到230。在230处，方法200输出在225处确定的组织的状态。例如，方法200可经由显示设备118在225处显示评估的结果。如果组织的状态经由如上文所述的温度分类来表征，则方法200可经由显示装置118显示至少一个ROI的温度分类。在一些示例中，温度分类可显示为例如与超声图像相邻的文本，或者叠加在ROI位置处的超声图像上。除此之外或另选地，温度分类可相对于超声图像的显示以图形方式显示(例如，作为热图)，以帮助用户可视化超声图像中组织的温度。通过以这种方式输出组织的状态，超声成像系统100的用户可在消融规程期间实时理解组织的状态。继而，用户可响应于组织的当前状态来调整消融规程(例如，通过结束规程或调整消融规程的参数)。
除了将组织的状态输出到显示设备118以显示给超声成像系统100的用户之外，方法200还可将组织的状态输出到存储器120，用于存储和后续查看。
在230处输出组织的状态之后，方法200继续到235。在235处，方法200确定超声扫描会话是否完成。方法200可确定例如超声扫描会话相对于例如经由用户界面115从用户接收输入命令是完成的，以结束超声扫描会话。在不存在用于结束超声扫描会话的这种命令的情况下，方法200可推测扫描未完成。如果超声扫描会话未完成(“否”)，则方法200返回到210。因此，对于在超声扫描会话期间采集的每个超声图像，方法200可继续采集超声图像，从所采集的超声图像中的ROI中提取特征，并且基于所提取的特征来评估组织的状态。应当理解，可保持在215处在超声图像中选择的ROI用于后续的超声图像。此外，从ROI提取的特征以及因此组织的状态可相对于在超声扫描会话期间采集的初始超声图像来确定。
最后，再次参考235，一旦扫描完成(“是”)，则方法200继续到240。在240处，方法200结束超声扫描会话。然后，方法200返回。因此，提供了一种用于组织的实时热监测的方法。该方法可被执行用于消融规程期间的热监测，例如，或者用于患者组织的一般非侵入性热监测。
图3示出了根据一个实施方案的示出用于从超声图像中的感兴趣区域提取纹理特征的示例性方法300的高级流程图。具体地，方法300涉及从在超声扫描会话期间采集的超声图像中的ROI提取特征，并且因此可包括上文所述的方法200的子例程。参照图1的系统和部件以及图2的方法200描述方法300，但是应当理解，在不脱离本公开的范围的情况下，方法300可用其他系统、部件和方法来实现。方法300可被实现为非暂态存储器诸如存储器120中的可执行指令，该可执行指令在由处理器诸如处理器116执行时使得处理器执行下文所述的动作。
方法300在305处开始。在305处，方法300将ROI转换成信号。具体地，对于每个ROI，方法300将二维ROI的所有像素值(例如，每个像素的灰度值)级联成单个一维矢量。在一些示例中，方法300使用多个轨线在转换期间遍历每个ROI的像素，使得每个ROI被转换成多个信号。例如，方法300可使用螺旋轨线(例如，从图像中心开始并向外螺旋)、反向螺旋轨线(例如，螺旋轨线的转置轨线)、之字形轨线(例如，沿着ROI的列遍历)和转置之字形轨线(例如，沿着ROI的行遍历)。本文参照图7进一步描述了在将ROI转换成信号时用于遍历像素的此类示例性轨线。
因此，作为例示性和非限制性示例，方法300可将每个ROI转换成四个对应信号。由此针对ROI获得的四个信号以四种不同的方式表示ROI的组织纹理。由于像素遍历轨线的多样性，单个ROI的信号的不同动态是可观察的，并且随着组织纹理变化(例如，由于消融规程)的时间推移，此类动态的演变，即使对于相同的轨线，也可如本文进一步所述的进行量化。
在频域中更容易观察到信号动态。因此，在310处继续，方法300将每个信号分解成多个频带，每个频带表示纹理的不同动态或方面。在一个示例中，方法300使用离散小波变换(DWT)分解每个信号或纹理信号。每个信号的信号趋势影响小波变换信号两端的奇异点，这可能对频带分解造成负面影响。因此，在分解之前，方法300首先通过计算每个信号的移动平均强度(例如，灰度值)并且从信号中减去移动平均强度来对每个信号执行偏置移除(即信号趋势的移除)。本文参照图10进一步描述了信号趋势或偏置移除的示例性信号。
方法300使用DWT中的适当标度选择将移除了偏置的每个信号分解成低频带(LF)、中频带(MF)、高频带(HF)和总频带(TB)。当从空间序列而不是时间序列生成信号时，方法300使用伪频率范围以将小波变换应用于信号。作为例示性的和非限制性的示例，可限定纹理信号中从0Hz到100Hz的频率范围，并且可为小波变换选择200Hz的采样频率。标度选择可针对被成像的组织的类型而不同，因为DWT的标度包含可特定于被成像的组织的类型的动态。用于DWT的特定标度可通过分析不同类型的组织的ROI数据的光谱学来确定，例如在整个消融序列上具有1至10的标度，这揭示了组织类型的基于频率的纹理动态。因此，可基于每个标度中观察到的频率或伪频率以及期望的频带选择标度以与DWT一起使用。例如，对于肝脏组织，标度1可涵盖具有显著变化的宽频率范围，这对于后续分析而言将是复杂的，因此标度1可被省略。同时，标度2和标度3可表现出期望的频率动态以解决例如12Hz至35Hz之间的高频带。类似地，标度4和标度5可表现出期望的频率动态以解决在3.5Hz和11Hz之间的中间频带，而标度6和标度7可表现出期望的频率动态以解决在0Hz和3.3Hz之间的低频带。由于低频信息或没有可观察到的频率观察，因此可忽略或省略高于7的标度。此外，总频带可通过2至7的所有所选择的标度的组合来选择。
因此，方法300使用离散的小波变换，其中为被成像的组织选择合适的标度参数，以将每个信号分解成四个窄带信号。应当理解，在不脱离本公开的范围的情况下，可使用除DWT之外的频率分解方法。如前所述，针对图像中的每个ROI获得四个纹理信号。在DWT之后，其中每个纹理信号被分解成四个信号，图像中的每个ROI有十六个信号。
在315处继续，方法300利用自回归模型生成每个分解信号的估计值。作为例示性和非限制性示例，方法300使用参数自回归(AR)模型与Yule-Walker方法对十六个窄带纹理信号进行建模。具体地，方法300通过由对应的估计AR参数计算功率谱密度(PSD)来对每个窄带信号的周期图进行建模。针对每个窄带信号单独地设置AR模型的顺序，使得用于给定窄带纹理信号的建模PSD(即，由对应的估计AR参数计算的PSD)最佳地匹配直接由窄带纹理信号估计的周期图。作为例示性和非限制性示例，AR顺序可包括LF信号的35、MF信号的25、HF信号的20和TB信号的30。AR顺序可例如基于对不同组织的样本组织纹理的研究预先确定，或者可在窄带纹理信号的建模期间通过方法300动态地确定。
所得的十六个AR PSD估计值以简化方式表示ROI内的组织纹理。在两个不同温度下的纹理之间，在AR PSD估计值中可容易地观察到每个频带中的频率偏移，并且导致此类频率偏移的纹理变化可归因于在消融规程期间施加到组织的热量。本文参照图13进一步描述了示例性周期图和AR PSD估计值。
在利用AR模型生成每个分解信号的估计值之后，方法300继续到320。在320处，方法300计算相同频带中的估计值和ROI到信号转换的相同像素轨线的相关系数。即，方法300计算每个估计值与初始估计值的相关系数。具体地，方法300可计算相同频带的一对AR PSD估计值的皮尔逊相关系数，包括频带的第一AR PSD估计值。因此，对于所采集的第一超声图像，针对第一AR PSD估计值与其自身计算给定频带和像素轨线的相关系数，由于精确的相关性，导致相关系数为1。所得相关系数可存储在相关系数标测图或矩阵中，例如如本文参照图14进一步所述。
彼此类似的AR PSD估计值产生较高的相关系数，而彼此不同的AR PSD估计值产生较低的相关系数。虽然最初所有AR PSD估计值是类似的并且因此高度相关，但是随着ROI的温度升高并且组织纹理响应于升高的温度而变化，在各种频带和转换类型中开始出现相异性。例如，AR PSD估计值中的相异性可开始出现在转置螺旋信号的窄带信号上，而之字形信号和转置之字形信号的窄带中的相异性在较高温度下增加。
由于每种转换方法的窄带中的相异性已经表现出作为温度增加的函数的某些行为，因此相关系数在一些示例中可用作将纹理变化与温度水平相关联的特征。然而，更稳健的方法可包括组合来自相关系数标测图或矩阵的所有相异性。为此，在325处继续，方法300确定相关系数的相异性度量。为了确定相关系数的相异性度量，方法300对每个ROI的相关系数进行排序，并且将经排序的系数组合成包括相异性度量的相异性矢量。例如，对于所采集的第一超声图像，相异性矢量或相异性度量包括所有值等于1的直线，因为相关系数指示ARP SD估计值与其自身的精确相关性。随着时间的推移，当温度升高并且各种相关系数远离一演变(即，随着ARP SD估计值在某些频带和转换类型中变得更相异)，相异性矢量表现出指示相异性增加的曲率。本文参照图15进一步描述了示例性相异性矢量。
由于随着时间和温度的推移，相异性度量易于区分，因此可容易地用某些特征来描述相异性度量。因此，在330处，方法300从相异性度量中提取特征。作为例示性和非限制性示例，方法300针对每个相异性度量计算多个特征。特征包括相异性度量的平均值和中值。相异性度量的平均值和中值均随温度升高而降低。此外，方法300针对每个相异性度量执行线性回归，并且确定线性回归的所得斜率和截距。因此，特征可包括线性回归的斜率和截距。随着温度升高，斜率升高并且截距减小。因此，方法300提取每个相异性度量的特征，包括相异性度量的中值、相异性度量的平均值、相异性度量的线性回归斜率以及相异性度量的线性回归截距。如上文参照图2所述，这些特征可用于对每个ROI的温度范围进行分类。然后，方法300返回。
因此，提供了一种用于从图像提取特征以评估组织的热状态的方法。该方法包括使用多个像素遍历轨线将ROI从二维图像转换成一维信号，将每个信号分解成多个频带，利用自回归模型生成每个分解信号的估计值，计算相同频带和轨线类型中的估计值的相关系数，确定相关系数的相异性度量，以及由相异性度量计算特征。应当理解，虽然本文所述的方法涉及选择性地处理超声图像以确定超声图像中描绘的组织的热状态，但是该方法可适于原始超声数据，诸如经由超声探头106采集的RF数据，而不是超声图像。
为了示出方法200的示例性实施方式，图4示出了示出利用超声成像系统对组织进行热监测的示例性方法400的框图。方法400可被实现为超声成像系统的非暂态存储器(诸如超声成像系统100的存储器120)中的可执行指令，并且可由超声成像系统的处理器(诸如超声成像系统100的处理器116)执行以执行本文所述的动作。
方法400采集参考超声图像帧402，并且选择参考纹理ROI(ROIT1)412。参考超声图像帧402包括处于未消融状态(即，在执行消融之前)的组织纹理结构。方法400采集新超声图像帧404，其包括与参考超声图像帧402相同的组织纹理结构但处于消融状态(即，在消融规程期间)。方法400在与参考纹理ROIT1 412相同的位置处选择新纹理ROI(ROITn)414。参考超声图像帧402和新超声图像帧404在具有相同解剖结构的相同图像平面中采集。
方法400计算相异性度量420，其包括ROI 412和ROI 414的相异性的测量。方法400如下文进一步所述计算相异性度量420，例如通过将每个ROI转换成一维信号，执行一维信号的趋势移除，将校正的信号分解成伪频带，对伪频带中的特定范围建模，计算不同ROI的对应模型之间的相关性，以及构建相关性矢量(例如，上文所述的相异性矢量或相异性度量)作为输出。相异性度量420包括ROI 412和414中描绘的纹理之间的相异性量度。由于ROI412中的纹理描绘了处于未消融状态的组织，而ROI 414中的纹理描绘了消融规程期间的相同组织，因此增加的相异性指示ROI 414中的组织相对于ROI 412中的组织的热状态的变化。如本文所述，这种相异性的定量量度可与温度直接相关联，使得ROI 414中的组织的热状态可被稳健地分类。此外，如果相异性度量420指示ROI 412和ROI 414之间的相似性而不是相异性，则方法400也可确定组织的热状态未改变。也就是说，相似性指示没有向组织施加热量。
再次参考图4，在计算相异性度量420之后，方法400接着从相关矢量执行特征提取425，例如通过从相关矢量中提取特征，该相关性矢量包括但不限于以下中的一者或多者：平均相关性、中值相关性、相关矢量的线性回归的斜率和截距、能量、能量差以及与第一相关矢量的比率等等。从特征提取425提取的特征存储在特征缓冲器427中。此外，从特征提取430提取的特征用作新纹理ROITn 414的分类430的输入。具体地，分类430使用所提取的特征的特性来对新纹理ROITn 414的热状态435进行分类。因此，热状态435包括分配给新纹理ROITn 414的表示新纹理ROITn 414中所描绘的组织的热状态类别的值。方法400在440处显示热状态435。例如，方法400可将新纹理ROITn 414的热状态435的指定值转换成颜色，并且用该颜色增强新超声图像帧404中的对应ROI 414。
基于所提取的特征的新纹理ROITn 414的热状态435的分类430更稳健，并且与更原始的热监测方法相比进一步降低了计算复杂性，例如，基于对ROI 412和ROI 414中描绘的纹理的直接分析，或者通过将ROI 412和/或ROI 414直接输入到一个或多个机器学习算法中。此外，当相异性度量420表征ROI 412和ROI 414之间的相异性时，在特征提取425期间提取的特征提供用于评估和分类组织的热状态的个性化基础。
此外，在一些示例中，分类430可对新纹理ROITn 414中的组织的非热状态445进行分类。例如，气泡的膨胀、组织收缩或组织膨胀可导致观察到的ROI内的组织纹理的偏移。因此，用作参考纹理ROIT1 412的参考超声图像帧402的ROI可基于新纹理ROITn 414内的组织的非热状态445进行更新。因此，方法400确定更新的参考纹理ROI 450，并且用更新的参考纹理ROI 450替换参考纹理ROIT1 412以用于后续监测。
为了使分类430对非热状态445进行分类，将对应于存储在特征缓冲器427中的先前ROI或最后已知热状态的特征输入到分类430，使得分类430可基于来自特征缓冲器427的先前特征与在特征提取425期间提取的当前特征的比较来识别非热状态445。
因此，提供了一种用于利用超声成像系统对组织进行热监测的方法，该方法包括基于ROI中的组织纹理与来自参考ROI的组织纹理的比较对ROI中的组织的热状态进行分类，其中由于组织的非热状态的变化，参考ROI随时间推移进行更新。
为了示出方法300的示例性实施方式，图5示出了根据一个实施方案的示出用于从超声图像中的感兴趣区域提取纹理特征的示例性方法500的框图。如上文所讨论的，用于利用超声成像对组织进行热监测的方法包括比较感兴趣区域随时间推移的纹理。例如，方法500包括将第一时间T1的ROI中的初始纹理505与第n时间Tn的ROI中的改变纹理507进行比较，其中ROI中的初始纹理505描绘于第一时间T1采集的超声图像中，而ROI中的改变纹理507描绘于第n时间Tn采集的超声图像中。值n包括超声序列(例如，诸如上文参照图2所述的超声扫描会话)中的图像的索引数量，因此第n时间Tn对应于采集第Tn超声图像时的时间。
为了比较纹理505和507，方法500执行纹理505和507的ROI到信号转换510。如上文参照图3所述，方法500通过根据多个像素遍历轨线(包括螺旋轨线、转置螺旋(T.Spiral)轨线、之字形轨线和转置之字形(T.Zigzag)轨线)级联ROI像素值(例如，灰度值)，将每个二维ROI转换成一维信号。因此，对于初始纹理505，转换510产生包括初始纹理505的不同一维表示的四个信号515。类似地，对于改变纹理507，转换510产生包括改变纹理507的不同的一维表示的四个信号517。
方法500接着将四个信号515和四个信号517中的每个信号分解成不同的频带(例如，如上文所述的LF、MF、HF和TB)，例如使用小波变换诸如DWT。当每个信号分解成四个不同的信号时，方法500因此针对初始纹理505由四个信号515产生十六个信号525，并且针对改变纹理507由四个信号517产生十六个信号527。
然后，方法500对十六个信号525以及十六个信号527执行AR建模530，以估计例如每个信号的PSD。AR建模530因此产生包括十六个信号525的AR PSD估计值的十六个信号535，以及包括十六个信号527的AR PSD估计值的十六个信号537。方法500计算十六个信号535和十六个信号537的相关系数540，产生指示十六个信号535和十六个信号537之间的相关性的相关系数542。具体地，在一些示例中被组织为一组4*4值的相关系数542包括具有相同频带和相同像素遍历轨线的信号之间的相关性的量度。例如，对于频带(例如，LF、MF、HF、TB)的每个组合和像素遍历轨线(例如，螺旋、转置螺旋、之字形、转置之字形)，相关系数542包括：指示源于初始纹理505和改变纹理507的螺旋LF信号的AR PSD估计值之间的相关性的相关系数，指示源于纹理505和507的螺旋MF信号的AR PSD估计值之间的相关性的相关系数，以及指示源于纹理505和507的螺旋HF信号的AR PSD估计值之间的相关性的相关系数等等。
然后，方法500确定相关系数542的相异性度量550，并且从相异性度量550中提取特征560。在特征提取560期间获得的特征包括改变纹理507相对于初始纹理505的定量量度，这继而可与ROI中的组织的热状态的变化相关联。例如，如上所述，在特征提取560期间获得的特征可用作一个或多个机器学习算法的输入，以对时间Tn处的ROI的温度进行分类。
应当理解，在超声成像序列或超声扫描会话期间采集第一超声图像之后，改变纹理507可包括与初始纹理505相同的纹理。此外，在初始生成十六个信号535之后，方法500可将初始纹理505的十六个信号535存储在本地存储器中，使得可仅处理在初始超声图像之后采集的超声图像的改变纹理507，以获得十六个信号537。
为了示出上文所述的方法300、400和500，图6至图18描绘了方法300、400和500的示例性应用。具体地，图6示出了根据一个实施方案的一组超声图像600，该组超声图像描绘了处于不同热状态的组织的超声图像的感兴趣区域。该组超声图像600包括第一超声图像602和第二超声图像622，其中第二超声图像622描绘了在执行消融规程的给定持续时间之后被成像的受试者的相同解剖区域。图像602和图像622均描绘了组织内的RF消融针607。此外，第一超声图像602包括多个ROI 610，其包括邻近RF消融针607的端部定位的至少一个ROI 612。多个ROI 610中的剩余ROI远离ROI 612分布。
此外，第二超声图像622还包括多个ROI 630，其包括邻近RF消融针607的端部定位的至少一个ROI 632。第二超声图像622中的多个ROI 630对应于第一超声图像602中的多个ROI 610。
该组超声图像600还包括包含第一超声图像602中的ROI 612的超声图像618，以及包含第二超声图像622中的ROI 632的超声图像638。超声图像618描绘ROI 612中的病变，其经由超声图像638中的RFA被破坏，从而产生超声图像618和638中可观察到的不同纹理。由于图像618和638之间的图像纹理的这种变化由于RFA而发生，因此纹理的差异可如本文所述进行定量，以将纹理的变化与温度的变化相关联。
具体地，如上文所述，分别用于初始纹理和改变纹理的超声图像618和超声图像638可从图6中描绘的二维表示转换成多个一维表示。为此，图7示出了一组示意图，其描绘了用于示例性感兴趣区域701的图像到信号转换的示例性像素轨线700。像素轨线700包括螺旋轨线710，其中从图像中心开始并且在ROI 701中沿逆时针方向螺旋的像素值被级联以形成一维信号或矢量。像素轨线700还包括转置螺旋轨线720，在本文中也称为反向螺旋轨线720，其类似地开始于ROI 701的中心处。像素轨线700还包括之字形轨线730，该之字形轨线从ROI 701的角落开始并沿着ROI701的列遍历。类似地，像素轨线700包括转置之字形轨线740，在本文中也称为反之字形轨线740，其从ROI 701的角落开始并沿着ROI 701的遍历，如所描绘的。
因此，对于图6中描绘的多个ROI 610和ROI 630中的每个ROI，可根据像素轨线700将像素值级联成信号，从而为每个ROI产生四个信号。
例如，图8示出了示出由单个感兴趣区域(诸如根据像素轨线700的ROI 612)构造的示例性信号的一组曲线图800，而图9示出了示出由单个感兴趣区域(诸如根据像素轨线700的ROI 632)构造的示例性信号的一组曲线图900。具体地，该组曲线图800包括：描绘根据螺旋轨线710转换的ROI 612的螺旋信号812的第一曲线图810，描绘根据之字形轨线730转换的ROI 612的之字形信号822的第二曲线图820，描绘根据转置螺旋轨线720转换的ROI612的转置螺旋信号832的第三曲线图830，以及描绘根据转置之字形轨线740转换的ROI612的转置之字形信号842的第四曲线图840。类似地，该组曲线图900包括：描绘根据螺旋轨线710转换的ROI 632的螺旋信号912的第一曲线图910，描绘根据之字形轨线730转换的ROI632的之字形信号922的第二曲线图920，描绘根据转置螺旋轨线720转换的ROI 632的转置螺旋信号932第三曲线图930，以及描绘根据转置之字形轨线740转换的ROI 632的转置之字形信号942的第四曲线图940。
如所描绘的，尽管由相同的ROI 612转换，但该组曲线图800中的曲线图清楚地指示不同的信号动态。类似地，尽管由相同的ROI 632转换，但该组曲线图900中的曲线图也清楚地指示不同的动态。
在如上文所述将每个信号分解成不同的频率(或伪频率)带之前，例如通过计算每个信号的移动平均强度(例如，灰度值)并且从信号中减去移动平均强度来从每个信号中移除信号趋势或偏置。作为例示性和非限制性示例，图10示出了一组曲线图1000，其示出了曲线图810中描绘的ROI 612的螺旋信号812以及曲线图910中描绘的ROI 632的螺旋信号912的信号偏置校正。具体地，该组曲线图1000包括：描绘包括偏置被移除的螺旋信号812的校正的螺旋信号1014的曲线图1012，以及描绘包括偏置被移除的螺旋信号912的校正的螺旋信号1024的曲线图1022。如所描绘的，校正的螺旋信号1014和1024仍然分别表现出螺旋信号812和螺912的信号动态，同时还改善了信号对小波分解的适用性。应当理解，虽然描绘了螺旋信号812和912的偏置移除或校正，但是针对两个ROI的剩余信号执行偏置移除。
在从信号中移除偏置之后，执行小波分解(例如，DWT)以将校正的信号分解成多个频带。作为例示性和非限制性示例，图11示出了根据一个实施方案的示出第一热状态的纹理信号到伪频带的示例性分解的一组曲线图1100。具体地，该组曲线图1100示出了由ROI612的初始纹理获得的校正的螺旋信号1014的分解。该组曲线图1100包括：描绘低频(LF)带中的DWT螺旋信号1112的第一曲线图1110、描绘中频(MF)带中的DWT螺旋信号1122的第二曲线图1120、描绘高频(HF)带中的DWT螺旋信号1132的第三曲线图1130以及描绘总频带(TB)中的DWT螺旋信号1142的第四曲线图1140。如上文所讨论的，使用DWT的标度参数6和7获得DWT螺旋信号1112，使用DWT的标度参数4和5获得DWT螺旋信号1122，使用DWT的标度参数2和3获得DWT螺旋信号1132，并且使用DWT的标度参数2至7获得DWT螺旋信号1142。
类似地，图12示出了根据一个实施方案的示出第二热状态的纹理信号到伪频带的示例性分解的一组曲线图1200。具体地，该组曲线图1200示出了由ROI 632中的改变纹理获得的校正的螺旋信号1024的分解。该组曲线图1200包括：描绘低频(LF)带中的DWT螺旋信号1212的第一曲线图1210、描绘中频(MF)带中的DWT螺旋信号1222的第二曲线图1220、描绘高频(HF)带中的DWT螺旋信号1232的第三曲线图1230以及描绘总频带(TB)中的DWT螺旋信号1242的第四曲线图1240。与图11的DWT螺旋信号一致，使用DWT的标度参数6和7获得DWT螺旋信号1212，使用DWT的标度参数4和5获得DWT螺旋信号1222，使用DWT的标度参数2和3获得DWT螺旋信号1232，并且使用DWT的标度参数2至7获得DWT螺旋信号1242。
因此，对于每种转换类型(即，用于将ROI转换成信号的每个像素遍历轨线)，信号被分解成四个不同频带，从而产生给定的四种转换类型的十六个信号。分解后的十六个信号都以不同的表示定量地描绘ROI的纹理，从而允许清楚地限定纹理的不同方面，该纹理的不同方面可能无法通过简单直接观察ROI中的图像纹理来分辨。
在信号分解之后，利用上文所述的AR模型来估计分解的信号。作为例示性且非限制性的示例，图13示出了一组曲线图1300，其示出了通过图11和图12的分解纹理信号的自回归建模获得的示例性功率谱密度估计值。该组曲线图1300包括：描绘通过源于ROI 612的DWT螺旋信号的AR建模确定的并且在图11中描绘的示例性功率谱密度估计值的第一组曲线图1305；以及描绘通过源于ROI 632的DWT螺旋信号的AR建模确定的并且在图12中描绘的示例性功率谱密度估计值的第二组曲线图1355。
具体地，第一组曲线图1305包括描绘通过LF带中的DWT螺旋信号1112的AR建模确定的PSD估计值1312的曲线图1310。曲线图1310还描绘了针对DWT螺旋信号1112生成的周期图1314。第一组曲线图1305还包括：描绘MF带中的DWT螺旋信号1122的AR PSD估计值1322和周期图1324的曲线图1320、描绘HF带中的DWT螺旋信号1132的AR PSD估计值1332和周期图1334的曲线图1330以及描绘TB中的DWT螺旋信号1142的AR PSD估计值1342和周期图1344的曲线图1340。
类似地，第二组曲线图1355包括：描绘LF带中的DWT螺旋信号1212的AR PSD估计值1362和周期图1364的曲线图1360、描绘MF带中的DWT螺旋信号1222的AR PSD估计值1372和周期图1374的曲线图1370、描绘HF带中的DWT螺旋信号1232的AR PSD估计值1382和周期图1384的曲线图1380以及描绘TB中的DWT螺旋信号1242的AR PSD估计值1392和周期图1394的曲线图1390。
如上文所述，用于组织的热监测的方法包括计算初始纹理和改变纹理之间的信号的相关系数。作为例示性示例，图14示出了根据一个实施方案的针对多个热状态计算的一组示例相关系数(CC)标测图1400。具体地，该组CC标测图1400描绘了针对ROI 612随时间推移以及因此ROI 612中组织的多个温度计算的相关系数。该组CC标测图1400中的每个CC标测图描绘了示出从零到一的相关系数的值的颜色，如图例1480所示。
该组CC标测图1400包括第一CC标测图1410，该第一CC标测图描绘了在27.5℃的第一温度下在ROI 612的AR PSD估计值之间计算的相关系数。具体地，针对AR PSD估计值与其自身计算每个相关系数，因此每个相关系数为1。
该组CC标测图1400还包括在30℃的第二温度下相同ROI的第二CC标测图1420、在40℃的第三温度下ROI的第三CC标测图1430、在50℃的第四温度下ROI的第四CC标测图1440，在60℃的第五温度下的ROI的第五CC标测图1450、在70℃的第六温度下的ROI的第六CC标测图1460以及在80℃的第七温度下的ROI的第七CC标测图1470。
如CC标测图1400所描绘的，观察到随着温度升高，相异性增加。CC标测图1410和1420例如描绘了在转置螺旋(T螺旋)转换的窄带中对小温度升高(例如，从27.5℃到30℃)的敏感性。CC标测图1430和1440指示对之字形和转置之字形转换中的温度升高的轻微敏感性。如CC标测图1450和1460所描绘的，对于在较高温度下升高温度，之字形和转置之字形转换的窄带中的相异性显著更敏感。螺旋转换的窄带中的相异性也似乎随着高于70℃的温度而增加，如CC标测图1470所示。
虽然CC标测图1400用于在视觉上描绘源于单个ROI的十六个通道中不断增加的相异性，并且虽然它们也可用作将纹理变化与温度变化相关联的特征，但是更稳健的方法包括将CC标测图的相异性组合成例如相异性度量或相异性矢量。
作为示例性和非限制性示例，图15示出了根据一个实施方案的示出不同热状态的经排序的相关系数的示例性相异性矢量的曲线图1500。CC标测图1400中的每个CC标测图按CC值排序并且组合成图1500中的相异性矢量。具体地，曲线图1500包括：包括CC标测图1410的经排序的相关系数的相异性矢量1510、包括CC标测图1420的经排序的相关系数的相异性矢量1520、包括CC标测图1430的经排序的相关系数的相异性矢量1530、包括CC标测图1440的经排序的相关系数的相异性矢量1540、包括CC标测图1450的经排序的相关系数的相异性矢量1550、包括CC标测图1460的经排序的相关系数的相异性矢量1560以及包括CC标测图1470的经排序的相关系数的相异性矢量1570。
相异性矢量1510对应于27.5℃下与其自身相比的初始纹理，并且因此是所有相关系数等于1的直线。组织温度增加至30℃产生相异性矢量1520相对于相异性矢量1510的曲率。后续相异性矢量1530、1540、1550和1560的曲率随着温度升高而继续增加，如通过最低系数(即，相异性矢量中的第一排序系数)特别指出的。相异性矢量1570指示纹理相异性从70℃(相异性矢量1560中示出了相关系数)到80℃的显著偏移。
如曲线图1500所描绘的，不同温度的相异性矢量易于区分。此外，相异性矢量可使用多种特征表征，该多种特征可从相异性矢量中提取。
作为一个示例，图16示出了一组曲线图1600，其示出了从图15的相异性矢量中提取的第一示例性特征。具体地，第一示例性特征包括相异性矢量的平均相关系数。该组曲线图1600包括曲线图1605，该曲线图描绘了第一相异性矢量1510的平均相关系数1610、第二相异性矢量1520的平均相关系数1612、第三相异性矢量1530的平均相关系数1614、第四相异性矢量1540的平均相关系数1616、第五相异性矢量1550的平均相关系数1618以及第六相异性矢量1560的平均相关系数1620。类似地，该组曲线图1600包括条形图1635，该条形图描绘了根据其对应组织温度的多个相异性矢量的平均相关系数。初始时，初始纹理的平均相关系数为1，并且平均相关系数随着温度升高而减小。
中值相关系数是可从相异性矢量中提取的另一个特征。作为示例性示例，图17示出了一组曲线图1700，该组曲线图示出了从图15的相异性矢量中提取的第二示例性特征。具体地，第二示例性特征包括相异性矢量的中值相关系数。该组曲线图1700包括曲线图1705，其描绘了第一相异性矢量1510的中值相关系数1710、第二相异性矢量1520的中值相关系数1712、第三相异性矢量1530的中值相关系数1714，第四相异性矢量1540的中值相关系数1716、第五相异性矢量1550的中值相关系数1718以及第六相异性矢量1560的中值相关系数1720。类似地，该组曲线图1700包括条形图1735，其描绘了根据其对应组织温度的多个相异性矢量的中值相关系数。初始时，初始纹理的中值相关系数为1，并且中值相关系数随着温度升高而减小。如所描绘的，作为温度的函数的中值相关系数的偏移小于平均相关系数的对应偏移，并且结合中值相关系数和平均相关系数的两个特征提供了将ROI的图像纹理与温度相关联的简洁方式。
可通过确定相异性矢量的线性回归来获得附加的特征。作为示例性示例，图18示出了一组曲线图1800，其示出了从图15的相异性矢量中提取的第三示例性特征和第四示例性特征。具体地，该组曲线图1800包括曲线图1805，其描绘了每个相异性矢量的示例性线性回归，包括相异性矢量1510的第一线性回归1810、相异性矢量1520的第二线性回归1812、相异性矢量1530的第三线性回归1814、相异性矢量1540的第四线性回归1816、相异性矢量1550的第五线性回归1818和相异性矢量1560的第六线性回归1820。该组曲线图1800还包括描绘曲线图1805中的每个线性回归的斜率的曲线图1835，和描绘曲线图1805中的每个线性回归的截距的曲线图1845。如曲线图1835所描绘的，斜率随着温度升高而增大。类似地，如曲线图1845所描绘的，截距随着温度升高而减小。
因此，通过处理ROI以获得ROI的初始纹理和改变纹理之间的相异性矢量，可从相异性矢量中提取多个特征，包括相异性矢量的平均相关系数、相异性矢量的中值相关系数、相异性矢量的线性回归的斜率和相异性矢量的线性回归的截距。这四个特征可用于对每个改变纹理的温度范围进行分类，如上文所述。
可训练上文所述的一个或多个机器学习算法以基于此类特征对组织的温度进行分类。例如，用于训练一个或多个机器学习算法的训练数据集可包括四个提取的特征和每个特征组合的标签。标签可指示组织纹理的温度，其可在创建训练数据集期间通过使用光纤或用于测量ROI的精确温度的另一装置来测量。除此之外或另选地，标签可指示对应于温度范围的值。例如，对于低于60℃的温度，标签可包括0，并且对于等于或高于60℃的温度，标签可包括1。训练数据集可随机置换进行重新排列，以对行进行随机排列，并且还分成两部分(例如，包括75％的数据集的第一部分和包括25％的数据集的第二部分)用于训练和测试。如上所述，该一个或多个机器学习算法可包括随机森林分类器(RFC)、支持矢量机(SVM)、人工神经网络(ANN)和适用于分类的其他机器学习算法中的一者或多者。
本公开的技术效果包括采集超声图像。本公开的另一个技术效果包括对超声图像中描绘的组织的温度进行分类。本公开的又一个技术效果包括基于对超声图像的分析来显示组织的热状态。本公开的另一个技术效果包括相对于先前采集的超声图像对超声图像中的图像纹理进行定量评估。
在一个实施方案中，方法包括经由超声探头采集超声图像，选择超声图像中的至少一个感兴趣区域，从至少一个感兴趣区域提取特征，基于特征对至少一个感兴趣区域中的组织的热状态进行分类，以及经由显示装置输出组织的热状态。
在该方法的第一示例中，从至少一个感兴趣区域提取特征包括将至少一个感兴趣区域转换成至少一个信号，将至少一个信号分解成多个频带，用自回归模型生成每个分解信号的估计值，计算每个估计值与相同频带中的初始估计值的相关系数，基于每个估计值的相关系数确定相异性度量，以及从相异性度量中提取特征。在任选地包括第一示例的方法的第二示例中，特征包括相异性度量的中值相关系数、相异性度量的平均相关系数、相异性度量的线性回归的斜率和相异性度量的线性回归的截距中的一者或多者。在任选地包括第一示例和第二示例中的一个或多个的方法的第三示例中，将至少一个感兴趣区域转换成至少一个信号包括，对于至少一个感兴趣区域中的感兴趣区域，在沿多个像素轨线遍历感兴趣区域的像素时对感兴趣区域的像素值进行采样，以生成感兴趣区域的多个信号，多个信号中的每个信号对应于多个像素轨线中的不同像素轨线。在任选地包括第一示例至第三示例中的一个或多个的方法的第四示例中，将至少一个信号分解成多个频带包括使用离散小波变换将至少一个信号分解成多个频带。在任选地包括第一示例至第四示例中的一个或多个的方法的第五示例中，利用自回归模型生成每个分解信号的估计值包括，对每个分解信号执行自回归建模以获得每个分解信号的自回归参数，以及基于相应的自回归参数计算每个分解信号的功率谱密度，估计值包括功率谱密度。在任选地包括第一示例至第五示例中的一个或多个的方法的第六示例中，该方法还包括计算每个估计值与相同频带中的初始估计值和相同图像到信号转换类型的相关系数。在任选地包括第一示例至第六示例中的一个或多个的方法的第七示例中，基于每个估计值的相关系数确定相异性度量包括将估计值的相关系数分类为矢量，以及从至少一个感兴趣区域提取特征包括从矢量计算特征。在任选地包括第一示例至第七示例中的一个或多个的方法的第八示例中，基于特征对至少一个感兴趣区域中的组织的热状态进行分类包括将特征输入到一个或多个机器学习算法中，以及从一个或多个机器学习算法接收组织的热状态的分类。在任选地包括第一示例至第八示例中的一个或多个的方法的第九示例中，分类包括温度分类或组织损伤可逆性分类中的一个或多个。在任选地包括第一示例至第九示例中的一个或多个的方法的第十示例中，选择超声图像中的至少一个感兴趣区域包括选择超声图像中的消融目标处的第一感兴趣区域，以及选择超声图像中远离消融目标定位的多个附加的感兴趣区域。
在另一个实施方案中，方法包括利用超声探头采集超声图像，将超声图像中的感兴趣区域转换成多个信号，将多个信号中的每个信号分解成多个频带，利用多个频带中的每个频带的自回归模型生成每个分解信号的估计值，测量每个分解信号的估计值与从先前采集的超声图像获得的对应估计值的相关性，基于每个分解信号的相关性计算多个特征，利用基于多个特征的机器学习算法对感兴趣区域中的组织的热状态进行分类，以及经由显示装置输出组织的热状态。
在该方法的第一示例中，将感兴趣区域转换成多个信号包括利用多个像素遍历轨线对感兴趣区域的像素值进行采样以生成多个信号。在任选地包括第一示例的方法的第二示例中，将多个信号中的每个信号分解成多个频带包括对每个信号应用离散小波变换，以将每个信号分解成低频带信号、中频带信号、高频带信号和总频带信号。在任选地包括第一示例和第二示例中的一个或多个的方法的第三示例中，多个特征包括以下中的一者或多者：中值相关系数、平均相关系数、包括每个分解信号的相关的矢量的线性回归的斜率以及矢量的线性回归的截距。在任选地包括第一示例至第三示例中的一个或多个的第四示例中，利用基于多个特征的机器学习算法对感兴趣区域中的组织的热状态进行分类包括相对于临界消融温度对组织的温度进行分类。
在又一个实施方案中，系统包括超声探头、显示装置和配置有非暂态存储器中的可执行指令的处理器，该可执行指令在被执行时使得处理器：经由超声探头采集超声图像；选择超声图像中的至少一个感兴趣区域；从至少一个感兴趣区域提取特征；基于特征对至少一个感兴趣区域中的组织的热状态进行分类；以及经由显示装置输出组织的热状态。
在该系统的第一示例中，处理器还配置有非暂态存储器中的可执行指令，该可执行指令在被执行时使得处理器将至少一个感兴趣区域转换成多个信号；将多个信号中的每个信号分解成多个频带；利用自回归模型生成每个分解信号的估计值；计算每个估计值与从相同频带中的先前采集的超声图像获得的初始估计值的相关系数；基于每个估计值的相关系数确定相异性度量；以及从相异性度量中提取特征。在任选地包括第一示例的系统的第二示例中，特征包括以下中的一者或多者：平均相关系数、中值相关系数、相异性度量的线性回归的斜率和相异性度量的线性回归的截距。在任选地包括第一示例和第二示例中的一个或多个的系统的第三示例中，处理器进一步配置有非暂态存储器中的可执行指令，该可执行指令在被执行时使得处理器：利用基于多个特征的机器学习算法对感兴趣区域中的组织的热状态进行分类，该热状态包括组织的相对于临界消融温度的温度。
如本文所用，以单数形式列举并且以单词“一个”或“一种”开头的元件或步骤应当被理解为不排除多个所述元件或步骤，除非明确说明此类排除。此外，对本发明的“一个实施方案”的引用不旨在被解释为排除也包含所引用特征的附加实施方案的存在。此外，除非明确地相反说明，否则“包含”、“包括”或“具有”具有特定特性的元件或多个元件的实施方案可包括不具有该特性的附加此类元件。术语“包括”和“在...中”用作相应的术语“包含”和“其中”的简明语言等同形式。此外，术语“第一”、“第二”和“第三”等仅用作标记，而不旨在对其对象施加数字要求或特定位置次序。
该书面描述使用示例来公开本发明，包括最佳模式，并且还使相关领域中的普通技术人员能够实践本发明，包括制造和使用任何设备或系统以及执行任何包含的方法。本发明可取得专利权的范围由权利要求书限定，并且可包括本领域普通技术人员想到的其他示例。如果此类其他示例具有与权利要求书的字面语言没有区别的结构元件，或者如果它们包括与权利要求书的字面语言具有微小差别的等效结构元件，则此类其他示例旨在落入权利要求书的范围内。

Claims (17)

1.一种方法，所述方法包括：
经由超声探头采集超声图像；
选择所述超声图像中的至少一个感兴趣区域；
从所述至少一个感兴趣区域中提取特征；
基于所述特征对所述至少一个感兴趣区域中的组织的热状态进行分类；以及
经由显示装置输出所述组织的所述热状态。
2.根据权利要求1所述的方法，其中从所述至少一个感兴趣区域提取所述特征包括：
将所述至少一个感兴趣区域转换成至少一个信号；
将所述至少一个信号分解成多个频带；
利用自回归模型生成每个分解信号的估计值；
计算每个估计值与相同频带中的初始估计值的相关系数；
基于每个估计值的所述相关系数确定相异性度量；以及
从所述相异性度量中提取所述特征。
3.根据权利要求2所述的方法，其中所述特征包括以下中的一者或多者：所述相异性度量的中值相关系数、所述相异性度量的平均相关系数、所述相异性度量的线性回归的斜率以及所述相异性度量的线性回归的截距。
4.根据权利要求2所述的方法，其中将所述至少一个感兴趣区域转换成所述至少一个信号包括：对于所述至少一个感兴趣区域中的感兴趣区域，在沿多个像素轨线遍历所述感兴趣区域的像素时对所述感兴趣区域的像素值进行采样，以生成所述感兴趣区域的多个信号，所述多个信号中的每个信号对应于所述多个像素轨线中的不同像素轨线。
5.根据权利要求2所述的方法，其中将所述至少一个信号分解成多个频带包括使用离散小波变换将所述至少一个信号分解成所述多个频带。
6.根据权利要求2所述的方法，其中利用所述自回归模型生成每个分解信号的所述估计值包括对每个分解信号执行自回归建模，以获得每个分解信号的自回归参数，以及基于所述相应的自回归参数计算每个分解信号的功率谱密度，所述估计值包括所述功率谱密度。
7.根据权利要求2所述的方法，还包括计算每个估计值与相同频带中并且相同图像到信号转换类型的初始估计值的所述相关系数。
8.根据权利要求2所述的方法，其中基于每个估计值的所述相关系数确定所述相异性度量包括将所述估计值的所述相关系数排序为矢量，并且其中从所述至少一个感兴趣区域提取所述特征包括从所述矢量计算所述特征。
9.根据权利要求1所述的方法，其中基于所述特征对所述至少一个感兴趣区域中的所述组织的所述热状态进行分类包括将所述特征输入到一个或多个机器学习算法中，以及从所述一个或多个机器学习算法接收所述组织的所述热状态的分类。
10.根据权利要求9所述的方法，其中所述分类包括温度分类或组织损伤可逆性分类中的一者或多者。
11.根据权利要求1所述的方法，其中选择所述超声图像中的所述至少一个感兴趣区域包括选择所述超声图像中的消融目标处的第一感兴趣区域，以及选择远离所述消融目标定位的所述超声图像中的多个附加的感兴趣区域。
12.一种方法，所述方法包括：
利用超声探头采集超声图像；
将所述超声图像中的感兴趣区域转换成多个信号；
将所述多个信号中的每个信号分解成多个频带；
利用所述多个频带中的每个频带的自回归模型生成每个分解信号的估计值；
测量每个分解信号的估计值与从先前采集的超声图像获得的对应估计值的相关性；
基于每个分解信号的所述相关性计算多个特征；
利用基于所述多个特征的机器学习算法对所述感兴趣区域中的组织的热状态进行分类；以及
经由显示装置输出所述组织的所述热状态。
13.根据权利要求12所述的方法，其中将所述感兴趣区域转换成所述多个信号包括利用多个像素遍历轨线对所述感兴趣区域的像素值进行采样以生成所述多个信号。
14.根据权利要求12所述的方法，其中将所述多个信号中的每个信号分解成所述多个频带包括对每个信号应用离散小波变换，以将每个信号分解成低频带信号、中频带信号、高频带信号和总频带信号。
15.根据权利要求12所述的方法，其中所述多个特征包括以下中的一者或多者：中值相关系数、平均相关系数、包括每个分解信号的所述相关性的矢量的线性回归的斜率以及所述矢量的所述线性回归的截距。
16.根据权利要求12所述的方法，其中利用基于所述多个特征的所述机器学习算法对所述感兴趣区域中的所述组织的所述热状态进行分类包括相对于临界消融温度对所述组织的温度进行分类。
17.一种系统，所述系统包括：
超声探头；
显示设备；和
处理器，所述处理器配置有非暂态存储器中的可执行指令，所述可执行指令在被执行时使得所述处理器执行根据权利要求1至16中任一项所述的方法。
CN202110744846.6A 2020-07-13 2021-06-30 利用超声成像系统对组织进行热监测的方法和系统 Pending CN113925598A (zh)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
US16/927,685 2020-07-13
US16/927,685 US11439308B2 (en) 2020-07-13 2020-07-13 Methods and systems for thermal monitoring of tissue with an ultrasound imaging system

Publications (1)

Publication Number Publication Date
CN113925598A true CN113925598A (zh) 2022-01-14

Family

ID=79171864

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110744846.6A Pending CN113925598A (zh) 2020-07-13 2021-06-30 利用超声成像系统对组织进行热监测的方法和系统

Country Status (2)

Country Link
US (1) US11439308B2 (zh)
CN (1) CN113925598A (zh)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US11439308B2 (en) * 2020-07-13 2022-09-13 GE Precision Healthcare LLC Methods and systems for thermal monitoring of tissue with an ultrasound imaging system
US11527329B2 (en) * 2020-07-28 2022-12-13 Xifin, Inc. Automatically determining a medical recommendation for a patient based on multiple medical images from multiple different medical imaging modalities
US20220301127A1 (en) * 2021-03-18 2022-09-22 Applied Materials, Inc. Image processing pipeline for optimizing images in machine learning and other applications
CN114587565B (zh) * 2022-03-01 2023-09-29 河南中医药大学 一种射频消融中温度控制方法及系统

Citations (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1303525A (zh) * 1997-10-10 2001-07-11 美国3M公司 膜电极组合件
US20110060221A1 (en) * 2009-09-04 2011-03-10 Siemens Medical Solutions Usa, Inc. Temperature prediction using medical diagnostic ultrasound
CN103040488A (zh) * 2012-12-21 2013-04-17 深圳大学 一种实时超声弹性成像位移估计方法和系统
JP2013244219A (ja) * 2012-05-25 2013-12-09 Fujifilm Corp 超音波診断装置、超音波診断方法、及び超音波診断プログラム
US20140180273A1 (en) * 2012-12-21 2014-06-26 Volcano Corporation Device, System, and Method for Imaging and Tissue Characterization of Ablated Tissue
US8852103B2 (en) * 2011-10-17 2014-10-07 Butterfly Network, Inc. Transmissive imaging and related apparatus and methods
US20150272653A1 (en) * 2014-04-01 2015-10-01 Siemens Medical Solutions Usa, Inc. Thermal therapy ablation detection with ultrasound
CN105793729A (zh) * 2013-07-23 2016-07-20 明尼苏达大学评议会 使用多频率波形的超声图像形成和/或重建
CN107123159A (zh) * 2017-01-23 2017-09-01 上海联影医疗科技有限公司 血流状态分析系统及方法
CN109890296A (zh) * 2016-10-27 2019-06-14 皇家飞利浦有限公司 具有组织类型分析器的超声系统
US20190282298A1 (en) * 2009-10-12 2019-09-19 Kona Medical, Inc. Intravascular energy delivery
CN111272305A (zh) * 2020-01-19 2020-06-12 南京大学 一种基于非线性热膨胀评估温度变化的超声方法及系统

Family Cites Families (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8376946B2 (en) * 2002-05-16 2013-02-19 Barbara Ann Karamanos Cancer Institute Method and apparatus for combined diagnostic and therapeutic ultrasound system incorporating noninvasive thermometry, ablation control and automation
US7367944B2 (en) 2004-12-13 2008-05-06 Tel Hashomer Medical Research Infrastructure And Services Ltd. Method and system for monitoring ablation of tissues
US8016757B2 (en) * 2005-09-30 2011-09-13 University Of Washington Non-invasive temperature estimation technique for HIFU therapy monitoring using backscattered ultrasound
WO2009045885A2 (en) * 2007-10-02 2009-04-09 Board Of Regents, The University Of Texas System Real-time ultrasound monitoring of heat-induced tissue interactions
US8864669B2 (en) * 2008-12-29 2014-10-21 Perseus-Biomed Inc. Method and system for tissue imaging and analysis
EP2797516B9 (en) * 2011-12-29 2017-04-19 St. Jude Medical, Atrial Fibrillation Division, Inc. System and device for ultrasound thermography during a cardiac ablation procedure
US20180310987A1 (en) * 2017-04-27 2018-11-01 Biosense Webster (Israel) Ltd. Systems and processes for map-guided automatic cardiac ablation
US20190053790A1 (en) * 2017-08-17 2019-02-21 Contraline, Inc. Systems and methods for automated image recognition of implants and compositions with long-lasting echogenicity
WO2019161914A1 (en) 2018-02-23 2019-08-29 Brainlab Ag Image based ultrasound probe calibration
US11439308B2 (en) * 2020-07-13 2022-09-13 GE Precision Healthcare LLC Methods and systems for thermal monitoring of tissue with an ultrasound imaging system

Patent Citations (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1303525A (zh) * 1997-10-10 2001-07-11 美国3M公司 膜电极组合件
US20110060221A1 (en) * 2009-09-04 2011-03-10 Siemens Medical Solutions Usa, Inc. Temperature prediction using medical diagnostic ultrasound
US20190282298A1 (en) * 2009-10-12 2019-09-19 Kona Medical, Inc. Intravascular energy delivery
US8852103B2 (en) * 2011-10-17 2014-10-07 Butterfly Network, Inc. Transmissive imaging and related apparatus and methods
JP2013244219A (ja) * 2012-05-25 2013-12-09 Fujifilm Corp 超音波診断装置、超音波診断方法、及び超音波診断プログラム
CN103040488A (zh) * 2012-12-21 2013-04-17 深圳大学 一种实时超声弹性成像位移估计方法和系统
US20140180273A1 (en) * 2012-12-21 2014-06-26 Volcano Corporation Device, System, and Method for Imaging and Tissue Characterization of Ablated Tissue
CN105793729A (zh) * 2013-07-23 2016-07-20 明尼苏达大学评议会 使用多频率波形的超声图像形成和/或重建
US20150272653A1 (en) * 2014-04-01 2015-10-01 Siemens Medical Solutions Usa, Inc. Thermal therapy ablation detection with ultrasound
CN105877779A (zh) * 2014-04-01 2016-08-24 美国西门子医疗解决公司 利用超声进行热疗消融检测
CN109890296A (zh) * 2016-10-27 2019-06-14 皇家飞利浦有限公司 具有组织类型分析器的超声系统
CN107123159A (zh) * 2017-01-23 2017-09-01 上海联影医疗科技有限公司 血流状态分析系统及方法
CN111272305A (zh) * 2020-01-19 2020-06-12 南京大学 一种基于非线性热膨胀评估温度变化的超声方法及系统

Also Published As

Publication number Publication date
US20220007946A1 (en) 2022-01-13
US11439308B2 (en) 2022-09-13

Similar Documents

Publication Publication Date Title
US11439308B2 (en) Methods and systems for thermal monitoring of tissue with an ultrasound imaging system
JP5303147B2 (ja) 弾性画像を生成する超音波診断装置
JP5334413B2 (ja) 超音波診断装置
EP3213108B1 (en) Imaging methods and apparatuses for performing shear wave elastography imaging
KR101579740B1 (ko) 초음파 진단장치, 그에 따른 초음파 진단 방법 및 그에 따른 컴퓨터 판독 가능한 저장매체
US10832405B2 (en) Medical image processing apparatus with awareness of type of subject pattern
CN112641464A (zh) 用于启用上下文感知的超声扫描的方法和系统
CN105188555A (zh) 超声波诊断装置以及图像处理装置
KR20170086311A (ko) 의료 영상 장치 및 그 동작방법
JP5113548B2 (ja) 超音波画像処理装置
Hussain et al. Strain-initialized robust bone surface detection in 3-D ultrasound
EP2599444B1 (en) Method and device for determining the elastic modulus of a biological tissue
US11559280B2 (en) Ultrasound imaging system and method for determining acoustic contact
JP7336766B2 (ja) 超音波診断装置、超音波診断方法および超音波診断プログラム
CN114159099A (zh) 乳腺超声成像方法及设备
JP2022507770A (ja) 熱焼灼領域の時間的マッピング
CN117224161A (zh) 一种超声成像方法以及超声成像系统
CN111281423A (zh) 一种超声图像优化方法和超声成像设备
CN112294361A (zh) 一种超声成像设备、盆底的切面图像生成方法
US20230414204A1 (en) Ultrasound diagnostic apparatus and control method of ultrasound diagnostic apparatus
JP4966226B2 (ja) 超音波画像処理装置及び方法並びにプログラム
US20220240898A1 (en) Ultrasound diagnostic apparatus, medical image processing apparatus, and medical image processing method
JP5396054B2 (ja) 超音波診断装置
US11227392B2 (en) Ultrasound imaging system and method
US11810294B2 (en) Ultrasound imaging system and method for detecting acoustic shadowing

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