CN113631953A - 使用基波和谐波信号的定量超声 - Google Patents

使用基波和谐波信号的定量超声 Download PDF

Info

Publication number
CN113631953A
CN113631953A CN202080022296.0A CN202080022296A CN113631953A CN 113631953 A CN113631953 A CN 113631953A CN 202080022296 A CN202080022296 A CN 202080022296A CN 113631953 A CN113631953 A CN 113631953A
Authority
CN
China
Prior art keywords
power spectrum
echo signal
signal power
frequency band
reference phantom
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
CN202080022296.0A
Other languages
English (en)
Inventor
Y·拉拜德
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.)
Siemens Medical Solutions USA Inc
Original Assignee
Siemens Medical Solutions USA Inc
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 Siemens Medical Solutions USA Inc filed Critical Siemens Medical Solutions USA Inc
Publication of CN113631953A publication Critical patent/CN113631953A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/52—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
    • G01S7/52017—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 particularly adapted to short-range imaging
    • G01S7/52023—Details of receivers
    • G01S7/52036—Details of receivers using analysis of echo signal for target characterisation
    • G01S7/52038—Details of receivers using analysis of echo signal for target characterisation involving non-linear properties of the propagation medium or of the reflective target
    • 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
    • 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
    • G01—MEASURING; TESTING
    • G01S—RADIO 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/00—Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
    • G01S15/88—Sonar systems specially adapted for specific applications
    • G01S15/89—Sonar systems specially adapted for specific applications for mapping or imaging
    • G01S15/8906—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques
    • G01S15/8909—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques using a static transducer configuration
    • G01S15/8915—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques using a static transducer configuration using a transducer array
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/52—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
    • G01S7/52004—Means for monitoring or calibrating
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/52—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
    • G01S7/52017—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 particularly adapted to short-range imaging
    • G01S7/52023—Details of receivers
    • G01S7/52036—Details of receivers using analysis of echo signal for target characterisation
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/52—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
    • G01S7/52017—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 particularly adapted to short-range imaging
    • G01S7/52085—Details related to the ultrasound signal acquisition, e.g. scan sequences
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/44—Constructional features of the ultrasonic, sonic or infrasonic diagnostic device
    • A61B8/4483—Constructional features of the ultrasonic, sonic or infrasonic diagnostic device characterised by features of the ultrasound transducer
    • A61B8/4494—Constructional features of the ultrasonic, sonic or infrasonic diagnostic device characterised by features of the ultrasound transducer characterised by the arrangement of the transducer elements
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/48—Diagnostic techniques
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/56—Details of data transmission or power supply

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • General Physics & Mathematics (AREA)
  • Biophysics (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Veterinary Medicine (AREA)
  • Public Health (AREA)
  • General Health & Medical Sciences (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Pathology (AREA)
  • Radiology & Medical Imaging (AREA)
  • Biomedical Technology (AREA)
  • Animal Behavior & Ethology (AREA)
  • Medical Informatics (AREA)
  • Molecular Biology (AREA)
  • Surgery (AREA)
  • Acoustics & Sound (AREA)
  • Physiology (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Nonlinear Science (AREA)
  • Ultra Sonic Daignosis Equipment (AREA)
  • Gynecology & Obstetrics (AREA)

Abstract

一种***和方法，包括：存储基波频带的参考体模的回波信号功率谱和谐波频带的参考体模的回波信号功率谱；获取基波频带的组织区域的回波信号功率谱和谐波频带的组织区域的回波信号功率谱；基于基波频带的组织区域的回波信号功率谱和基波频带的参考体模的回波信号功率谱来确定第一反向散射系数；基于谐波频带的组织区域的回波信号功率谱和谐波频带的参考体模的回波信号功率谱来确定第二反向散射系数；以及基于第一反向散射系数和第二反向散射系数来确定组织区域的非线性度。

Description

使用基波和谐波信号的定量超声
背景技术
常规的超声成像***通过检测体积内的声阻抗不连续性来创建该体积的内部图像（即，B模式图像）。更具体地，常规的超声成像涉及将超声波束发射到体积中，并且检测从该体积内的声阻抗不连续性反射的返回信号。由于不同的材料通常表现出不同的声阻抗，因此检测到的声阻抗不连续性表示不同材料在体积内的位置。
上面描述的B模式图像描绘了组织中的定性特征，但是不直接提供关于组织的定量信息。另一方面，定量超声***可以确定感兴趣区域（ROI）的衰减系数（AC）和反向散射系数（BSC）。AC是组织中超声能量损失的量度，并且BSC是从组织返回的超声能量的量度。诸如这些的定量超声值有助于表征组织的属性，诸如硬度和脂肪分数（fat fraction）。附加的定量超声量度将进一步辅助这些表征。
AC和BSC的确定需要针对获取超声***的***效应进行补偿。这些效应包括发射脉冲功率、换能器灵敏度、波束聚焦图样（pattern）、以及深度相关的接收器增益。常规地，补偿包括将基波频带中的组织样本的回波信号功率谱除以来自相同深度的基波频带中的良好表征的参考体模（phantom）的回波信号功率谱。所得的商是归一化谱，其取决于组织样本和参考体模的衰减和反向散射属性。由于参考体模的属性是已知的，因此可以从归一化谱中导出组织样本的AC和BSC。
前述确定的一个缺点是需要在临床扫描时获取参考体模数据。该附加获取阻碍了工作流程和患者吞吐量。此外，这种确定的准确性在许多场景中已经被视为是缺乏的。前述确定也未能提供其他潜在有用的定量量度，诸如组织非线性度。
附图说明
根据考虑如附图中所图示的以下说明书，实施例的构造和使用将变得明显，在附图中，相同的参考标号标示相同的部件，并且其中：
图1是根据一些实施例的定量超声***的框图；
图2是根据一些实施例的用于确定定量超声值的过程的流程图；
图3是根据一些实施例的用于获取参考体模校准数据的超声***的框图；
图4是根据一些实施例的用于确定参考体模校准数据的数据的表格表示；
图5是根据一些实施例的包括定量超声值的超声图像；以及
图6是根据一些实施例的超声***的框图。
具体实施方式
提供以下描述是为了使得本领域的任何人员能够制作和使用所描述的实施例，并且阐述了设想用于实行所描述的实施例的最佳模式。然而，各种修改对于本领域技术人员来说将仍然是明显的。
一些实施例提供了对定量超声值的高效且准确的确定。更具体地，一些实施例提供了一种基于来自谐波频带的回波信号功率谱来确定AC和BSC的创造性***。使用谐波频带在临床上更合期望。使用谐波信号的临床益处是由于减少的混响、杂波噪声和相位像差而改进的QUS估计的重复性和再现性。
这种实施例可以利用谐波频带中的良好表征的参考体模的预先存储的回波信号功率谱（或可以从其中计算功率谱的RF或IQ信号数据），从而消除了在临床数据获取时或接近临床数据获取时获取校准数据的需要。
一些实施例有利地确定组织非线性度的定量值。相比之下，上述常规***假设参考体模的非线性度基本上等于ROI中的组织的非线性度。
图1图示了根据一些实施例的实现方式。***100包括超声单元110、超声换能器120和显示器130。体积140可以包括人体，但是实施例不限于此。超声换能器100可以包括任何合适的超声换能器，诸如但不限于相控阵、线性或凸形超声换能器。
一般地，单元110的处理单元112可以执行程序代码来控制换能器120，以将超声波束发射到体积140中，并且从其中接收声学射频信号。单元110的处理单元112可以执行程序代码，以基于接收到的信号来生成图像和/或确定定量超声值。图像和/或所确定的值可以在显示器130上显示给技术人员。
根据一些实施例，技术人员对***100进行操作，以从体积140的ROI获取谐波频带的回波信号功率谱。获取可以包括：获取RF或IQ信号，并且从其中计算谐波频带的回波信号功率谱。使用第一扫描设置来执行该获取，第一扫描设置可以包括频率、F数/孔径大小、聚焦和变迹（apodization）函数参数的特定值。换能器120与发射单元116结合操作以将超声波束发射到ROI中，并且接收器单元118与换能器120结合操作以根据扫描设置从ROI接收谐波频带中的反射信号。
接下来，从存储装置114中检索对应于第一扫描设置的谐波校准数据。谐波校准数据可以包括使用与第一扫描设置相同或基本相似的扫描设置从参考体模获取的谐波频带的回波信号功率谱（或可以从其中导出回波信号功率谱的数据）。在下面详细描述了谐波校准数据的获取和存储。
通过将从体积140获取的回波信号功率谱除以所存储的回波信号功率谱，来对从体积140获取的回波信号功率谱进行归一化。结果是归一化谱，其仅取决于ROI组织和参考体模的衰减和反向散射属性。由于参考体模的衰减和反向散射属性是已知的，因此可以从归一化谱中导出组织的AC和BSC，如将在下面详细描述的那样。
根据一些实施例，还从体积140的ROI获取基波频带的回波信号功率谱。使用第二扫描设置来执行该获取，第二扫描设置可以不同于或者并非不同于上面提及的第一扫描设置。从存储装置114中检索对应于第二扫描设置的基波校准数据。基波校准数据包括使用与第二扫描设置相同或基本相似的扫描设置从参考体模获取的基波频带的回波信号功率谱（或者，同样地，可以从其中导出回波功率信号谱的RF或IQ数据）。
通过将从体积140获取的基波频带回波信号功率谱除以基波校准数据，来对从体积140获取的基波频带回波信号功率谱进行归一化，从而得到另一个归一化谱，其仅取决于ROI组织和参考体模的衰减和反向散射属性。也可以从该归一化谱中导出组织的AC和BSC。
根据一些实施例，可以基于从谐波频带数据中导出的BSC和从基波频带数据中导出的BSC来确定ROI的非线性度。这种确定提供了相比于常规***改进的ROI表征，常规***假设参考体模和ROI的非线性度是等效的。
图2是根据一些实施例的用于确定定量超声值的过程200的流程图。过程200可以由***100的元件执行，但是实施例不限于此。过程200和本文中提及的所有其他过程可以体现在从一个或多个非暂时性计算机可读介质（诸如，硬盘、易失性或非易失性随机存取存储器、DVD-ROM、闪存驱动器和磁带）读取的处理器可执行程序代码中，并且然后以压缩、未编译和/或加密的格式来存储。在一些实施例中，可以使用硬连线电路来代替程序代码、或者与程序代码组合地使用硬连线电路，以用于实现根据一些实施例的过程。因此，实施例不限于硬件和软件的任何特定组合。
初始地，在S210处，扫描参考体模以获取RF或IQ数据，并且从其中确定基波频带和谐波频带的回波信号功率谱。使用特定的扫描设置来扫描参考体模。在一些实施例中，使用其他扫描设置在S210处获取基波频带和谐波频带的附加回波信号功率谱。
图3图示了根据一些实施例的用于执行S210的***300。***300可以由超声***提供者、体模提供者或其他非临床实体来操作。超声单元310的处理单元312执行存储装置314的扫描程序，以控制换能器320向体模340发射信号和从体模340接收信号。体模340表示预期的ROI（例如，成年男性躯干），并且超声单元310和换能器320可以包括预期将在扫描预期ROI中使用的超声单元和换能器的生产等效的版本。存储装置314存储所获取的功率谱。
接下来，在S220处，将基波频带和谐波频带的所获取的回波信号功率谱存储在超声***中。在一些实施例中，所存储的数据包括可以从其中导出基波频带和谐波频带的回波信号功率谱的数据。在一些示例中，该谱被存储为要运送给客户的超声***的上述基波和谐波校准数据。根据一些实施例，该谱与用于获取相应谱的扫描设置相关联地存储。
图4是根据一些实施例的在S220处存储的数据的表格表示。表格400将每个所获取的功率谱与用于获取该功率谱的扫描设置和参考体模相关联。如所示的，在S210处，可以扫描一个以上类型的参考体模。每个参考体模/扫描设置对与基波频带功率谱数据和谐波频带功率谱数据相关联。谱数据列的值可以由包括对应谱数据的文件的文件名组成。
在一些实施例中，S210和S220在远在对其中获取的数据进行临床使用之前（例如，在将超声***运送到诊所之前）的时间段期间被执行。因此，S220与S230之间的虚线箭头指示可能比过程200的其他相邻步骤之间的时间显著更长的时间流逝。
在S230处，例如在临床设置中，对其中存储了该谱的超声***进行操作，以从ROI获取基波频带和谐波频带的回波信号功率谱。该获取使用第一扫描设置，第一扫描设置用于在S210处生成对应参考体模的校准数据。在一些实施例中，第一扫描设置被设置为超声***的默认扫描设置。
根据S230的一些实施例，超声***以3 MHz的中心频率以及从2至4 MHz的频率带宽来发射信号。然后，在2-4 MHz范围内接收基波信号。为了获得谐波频带中的信号，以1.5MHz的中心频率来发射信号，并且该信号具有1-2 MHz之间的带宽。所得的接收到的谐波信号可以表现出例如中心频率的两倍或者在2-4 MHz之间的频带。因此，QUS值的后续计算与2-4 MHz之间的频率相关联。
接下来，在S240处，确定对应于第一扫描设置的所存储的校准数据。所确定的校准数据由使用对应于ROI的参考体模和第一扫描设置所获取的基波频带和谐波频带的回波信号功率谱组成。例如，S240可以包括：在表格400的一行内标识适当的参考体模和扫描设置，并且获取在表格400的同一行内标识的所存储的基波和谐波校准数据文件。
在S250处，确定第一AC和第一BSC。该确定基于在S230处获取的基波频带的回波信号功率谱和在S240处确定的基波频带的校准回波信号功率谱。实施例不限于以下对S250处的确定的描述。
初始地，将ROI的每个射频回波线划分成若干个重叠的时间门控窗口。傅里叶变换被应用于每个窗口，并且对应于相同深度的窗口的功率谱被求平均。对参考体模的对应ROI执行相同的过程。在标准脉冲回波成像中，统计上均质的组织中的窗口化区域的基波频带中的测量功率谱由等式（1）给出：
Figure 432430DEST_PATH_IMAGE002
（1）。
下标s表示样本（即，ROI的组织）。从换能器表面到ROI内特定时间门控窗口的中心的距离由z来标示。频率由f来标示。T(f)表示发射脉冲的传递函数。E tx (f)和E rx (f)分别表示换能器电-声和声-电传递函数。D(f,z)标示衍射效应，衍射效应与换能器几何形状以及发射和接收聚焦相关。
Figure DEST_PATH_IMAGE004A
和BSC s (f)分别是样本的频率相关AC和BSC值。
类似地，来自参考体模的反向散射信号的功率谱是：
Figure 975670DEST_PATH_IMAGE006
（2）。
将样本的功率谱除以参考体模的功率谱得出：
Figure 755407DEST_PATH_IMAGE008
（3）。
补偿参考体模的已知衰减和反向散射属性，等式（3）变为：
Figure 150616DEST_PATH_IMAGE010
（4）。
计算自然对数得出：
Figure 964988DEST_PATH_IMAGE012
（5）。
然后，可以从比对深度z来拟合等式（5）的线的斜率（−4
Figure DEST_PATH_IMAGE013
）和截距（ln(BSC s (f))）中导出衰减系数
Figure DEST_PATH_IMAGE013A
（np/cm）和反向散射系数BSC s （1/cm-str）。
在S260处，确定第二AC和第二BSC。S260处的确定基于在S230处获取的谐波频带的回波信号功率谱和在S240处确定的谐波频带的校准回波信号功率谱。实施例不限于以下对S260处的确定的描述。
使用谐波频带来估计AC和BSC需要新模型，该新模型计及组织非线性度。来自基波频率f下的平面波P 0的二次谐波压力由下式给出：
Figure DEST_PATH_IMAGE015
（6），
其中
Figure DEST_PATH_IMAGE017
和
Figure DEST_PATH_IMAGE019
分别是基波和谐波信号的衰减系数（np/cm），并且K是与非线性度参数B/A成比例的常数。
可以使用指数函数的泰勒级数来进一步简化等式（6）中的比率项：
Figure DEST_PATH_IMAGE021
，
其中最后的步骤使用了常见假设
Figure DEST_PATH_IMAGE023
和泰勒级数近似。
在使用临床换能器阵列的聚焦发射的情况下，P 0(f)被建模为：
Figure DEST_PATH_IMAGE025
（8），
其中T(f)是发射脉冲的传递函数，E tx (f)是发射时的电子-机械（electro-mechanical）传递函数，并且D tx (f)是发射时的衍射图样。使用等式（7）和（8）：
Figure DEST_PATH_IMAGE027
（9）。
接收到的二次谐波信号的功率谱由下式给出：
Figure DEST_PATH_IMAGE029
（10），
其中E rx (f)是接收时的电子-机械传递函数，并且D rx (f)是接收时的衍射图样，BSC(2f)是二次谐波频率下的反向散射系数（np/cm-str），并且
Figure DEST_PATH_IMAGE031
计及返回谐波信号的衰减。对等式（9）和（10）进行组合：
Figure DEST_PATH_IMAGE033
（11）。
在组织谐波成像（THI）模式中，对来自180°异相（180° out of phase）的两个脉冲的射频信号进行求和，以获得谐波信号。
将ROI的每个射频回波线划分成若干个重叠的时间门控窗口。傅里叶变换被应用于每个窗口，并且对应于相同深度的窗口的功率谱被求平均。对参考体模的对应ROI执行相同的过程。
ROI中的样本和参考的功率谱由等式（11）给出。假设括号中的项对于样本和参考两者是相同的。通过计算来自样本和参考的功率谱的比率，我们获得了类似于等式（3）的等式：
Figure DEST_PATH_IMAGE035
（12）。
补偿参考体模的已知衰减、反向散射和非线性度，等式（12）变为：
Figure DEST_PATH_IMAGE037
（13）。
计算自然对数得出：
Figure DEST_PATH_IMAGE039
（14）。
可以从比对深度z来拟合等式（14）的线的斜率中导出频率2f下的衰减系数
Figure DEST_PATH_IMAGE041
（dB/cm）。比对深度z来拟合等式（14）的线的截距得出了由非线性度项
Figure DEST_PATH_IMAGE043
偏置的ln(BSC s (2f))。
返回到过程200，在S270处，基于在S250处确定的第一BSC和在S260处确定的第二BSC来确定ROI的非线性度。继续上面的示例，可以通过基于等式（5）使用基波频带来测量ln(BSC s (2f))（其中基波频率现在是2f）并且将结果代入等式（14）中，从而估计出等式（14）中的非线性度项
Figure DEST_PATH_IMAGE043A
。
在S280处，生成并且显示ROI的图像。如本领域中已知的，可以基于在S230处获取的谱中的一个或其两者来生成该图像。该图像还可以指示非线性度以及基于接收到的信号所确定的任何其他值。
图5示出了根据一些实施例的在S280处生成并显示的图像500。如所示的，剪切波速（Vs）、弹性（E）、AC、BSC、以超声方式导出的脂肪分数（UDFF）以及非线性度（K）的所确定的值与超声B模式图像数据同时显示。这些定量测量结果可以通过移除对B模式图像的定性解释以及通过减少***相关因素来改进医学超声的诊断能力。
基于基波带和谐波带的对AC的估计应当是等效的。因此，图像500可以显示这两者的平均值或估计值中的任一个。在一些实施例中，可以通过确定并显示两个AC估计的加权平均值来减少可变性。
在BSC的情况下，基于谐波带在S260处确定的BSC值由非线性度项所偏置。因此，在S280处显示的BSC可以是基于基波带在S250处确定的BSC值。如果假设非线性度项是可忽略不计的，那么可以显示这两个BSC值的平均值。
图6是根据一些实施例的超声成像***600的框图。***600可以实现本文中描述的一个或多个过程。***600是相控阵超声成像***，但是实施例不限于此。典型的相控阵***利用64到256个接收信道和相当数量的发射信道。为了清楚起见，图6图示了单个发射和接收信道。
***600包括换能器元件605和发射/接收开关610。换能器元件605可以包括压电或电容膜元件的1维、1.25维、1.5维、1.75维或2维阵列的元件。操作发射/接收开关610以允许经由元件605来发射超声能量（例如，响应于跨元件605的电压施加），或者允许接收由元件605响应于接收到的超声能量（即，回波）而生成的电压。
发射波束形成器615结合数模转换器620和高压发射器625而可操作，以生成用于多个信道的波形，其中每个波形可以表现出不同的幅度、延迟和/或相位。接收波束形成器630从多个信道接收信号，每个信道可以经受放大635、滤波640、模数转换645、延迟和/或相位旋转器以及一个或多个求和器。接收波束形成器630可以由硬件或软件来配置，以响应于每个发射波束来应用相对延迟、相位和/或变迹以形成一个或多个接收波束。接收波束形成器630可以提供本领域已知的动态接收聚焦、以及固定聚焦接收。
由接收波束形成器630形成的接收波束表示发射波束和接收波束已经穿过的材料。接收波束被输出到处理器650以供处理。例如，处理器650可以基于接收波束来生成图像。
处理器650可以执行存储在存储器660中的处理器可执行程序代码，以执行和/或控制***600的其他组件来执行本文中描述的过程。处理器650可以包括B模式检测器、多普勒检测器、脉冲波多普勒检测器、相关性处理器、傅里叶变换处理器、专用集成电路、通用处理器、控制处理器、图像处理器、现场可编程门阵列、数字信号处理器、模拟电路、数字电路、它们的组合、或者用于基于波束形成的超声样本来生成数据（例如，图像数据）的其他当前已知或以后开发的设备。
存储器660可以包括非暂时性计算机可读存储介质，诸如随机存取存储器和/或非易失性存储器（例如，闪速存储器、硬盘存储器）。存储器660可以存储程序代码、校准数据、B模式图像和/或任何其他合适的数据。显示器655可以包括阴极射线管显示器、液晶显示器、发光二极管显示器、等离子显示器、或用于显示图像和/或测量值的其他类型的显示器。
本领域技术人员将领会，在不脱离权利要求的范围和精神的情况下，可以配置上述实施例的各种改编和修改。因此，要理解的是，权利要求可以不同于本文中具体描述的那样来实践。

Claims (18)

1.一种超声***，包括：
存储器，其存储基波频带的参考体模的回波信号功率谱和谐波频带的参考体模的回波信号功率谱；以及
处理单元，用于：
控制所述超声***来获取基波频带的组织区域的回波信号功率谱和谐波频带的组织区域的回波信号功率谱；
基于基波频带的组织区域的回波信号功率谱和基波频带的参考体模的回波信号功率谱来确定第一反向散射系数；
基于谐波频带的组织区域的回波信号功率谱和谐波频带的参考体模的回波信号功率谱来确定第二反向散射系数；以及
基于第一反向散射系数和第二反向散射系数来确定组织区域的非线性度。
2.根据权利要求1所述的超声***，其中，使用基本相似的扫描设置来获取基波频带的参考体模的回波信号功率谱和基波频带的组织区域的回波信号功率谱。
3.根据权利要求2所述的超声***，其中，使用基本相似的扫描设置来获取谐波频带的参考体模的回波信号功率谱和谐波频带的组织区域的回波信号功率谱。
4.根据权利要求1所述的***，所述存储器存储基波频带的第二参考体模的回波信号功率谱和谐波频带的第二参考体模的回波信号功率谱，以及
所述处理单元进一步基于组织区域与参考体模之间的相关性来确定使用基波频带的参考体模的回波信号功率谱以用于确定第一反向散射系数、以及使用谐波频带的参考体模的回波信号功率谱以用于确定第二反向散射系数。
5.根据权利要求4所述的***，已经使用第一扫描设置获取了基波频带的参考体模的回波信号功率谱和谐波频带的参考体模的回波信号功率谱，
所述存储器存储基波频带的参考体模的第二回波信号功率谱和谐波频带的参考体模的第二回波信号功率谱，第二回波信号功率谱是使用第二扫描设置来获取的，
使用第三扫描设置来获取基波频带的组织区域的回波信号功率谱和谐波频带的组织区域的回波信号功率谱，以及
其中所述处理单元用于基于第一扫描设置与第三扫描设置之间的相关性来确定使用基波频带的参考体模的回波信号功率谱和谐波频带的参考体模的回波信号功率谱，以用于确定第一反向散射系数和确定第二反向散射系数。
6.根据权利要求1所述的***，所述处理单元进一步用于：
基于谐波频带的组织区域的回波信号功率谱、谐波频带的参考体模的回波信号功率谱、以及参考体模的衰减系数来确定组织区域的衰减系数。
7.根据权利要求6所述的***，进一步包括：
显示器，用于同时显示第二反向散射系数、组织区域的衰减系数、非线性度、以及组织区域的B模式图像。
8.一种方法，包括：
在超声***中存储基波频带的参考体模的回波信号功率谱和谐波频带的参考体模的回波信号功率谱；
控制所述超声***来获取基波频带的组织区域的回波信号功率谱和谐波频带的组织区域的回波信号功率谱；
基于基波频带的组织区域的回波信号功率谱和基波频带的参考体模的回波信号功率谱来确定第一反向散射系数；
基于谐波频带的组织区域的回波信号功率谱和谐波频带的参考体模的回波信号功率谱来确定第二反向散射系数；以及
基于第一反向散射系数和第二反向散射系数来确定组织区域的非线性度。
9.根据权利要求8所述的方法，其中，使用基本相似的扫描设置来获取基波频带的参考体模的回波信号功率谱和基波频带的组织区域的回波信号功率谱。
10.根据权利要求9所述的方法，其中，使用基本相似的扫描设置来获取谐波频带的参考体模的回波信号功率谱和谐波频带的组织区域的回波信号功率谱。
11.根据权利要求8所述的方法，进一步包括：
存储基波频带的第二参考体模的回波信号功率谱和谐波频带的第二参考体模的回波信号功率谱，以及
基于组织区域与参考体模之间的相关性来确定使用基波频带的参考体模的回波信号功率谱以用于确定第一反向散射系数、以及使用谐波频带的参考体模的回波信号功率谱以用于确定第二反向散射系数。
12.根据权利要求11所述的方法，已经使用第一扫描设置获取了基波频带的参考体模的回波信号功率谱和谐波频带的参考体模的回波信号功率谱，进一步包括：
在所述超声***中存储基波频带的参考体模的第二回波信号功率谱和谐波频带的参考体模的第二回波信号功率谱，第二回波信号功率谱是使用第二扫描设置来获取的，
其中使用第三扫描设置来获取基波频带的组织区域的回波信号功率谱和谐波频带的组织区域的回波信号功率谱，以及
进一步包括：基于第一扫描设置与第三扫描设置之间的相关性来确定使用基波频带的参考体模的回波信号功率谱以用于确定第一反向散射系数、以及使用谐波频带的参考体模的回波信号功率谱以用于确定第二反向散射系数。
13.根据权利要求8所述的方法，进一步包括：
基于谐波频带的组织区域的回波信号功率谱、谐波频带的参考体模的回波信号功率谱、以及参考体模的衰减系数来确定组织区域的衰减系数。
14.根据权利要求13所述的方法，进一步包括：
同时显示第二反向散射系数、组织区域的衰减系数、非线性度、以及组织区域的B模式图像。
15.一种超声成像***，用于：
存储基波频带的参考体模的回波信号功率谱和谐波频带的参考体模的回波信号功率谱；
获取基波频带的组织区域的回波信号功率谱和谐波频带的组织区域的回波信号功率谱；
基于基波频带的组织区域的回波信号功率谱和基波频带的参考体模的回波信号功率谱来确定第一反向散射系数；
基于谐波频带的组织区域的回波信号功率谱和谐波频带的参考体模的回波信号功率谱来确定第二反向散射系数；以及
基于第一反向散射系数和第二反向散射系数来确定组织区域的非线性度。
16.根据权利要求15所述的***，其中，使用基本相似的扫描设置来获取基波频带的参考体模的回波信号功率谱和基波频带的组织区域的回波信号功率谱。
17.根据权利要求15所述的***，进一步用于：
存储基波频带的第二参考体模的回波信号功率谱和谐波频带的第二参考体模的回波信号功率谱，以及
基于组织区域与参考体模之间的相关性来确定使用基波频带的参考体模的回波信号功率谱以用于确定第一反向散射系数、以及使用谐波频带的参考体模的回波信号功率谱以用于确定第二反向散射系数。
18.根据权利要求17所述的***，已经使用第一扫描设置获取了基波频带的参考体模的回波信号功率谱和谐波频带的参考体模的回波信号功率谱，所述***用于：
存储基波频带的参考体模的第二回波信号功率谱和谐波频带的参考体模的第二回波信号功率谱，第二回波信号功率谱是使用第二扫描设置来获取的，
其中使用第三扫描设置来获取基波频带的组织区域的回波信号功率谱和谐波频带的组织区域的回波信号功率谱，以及
所述***用于基于第一扫描设置与第三扫描设置之间的相关性来确定使用基波频带的参考体模的回波信号功率谱以用于确定第一反向散射系数、以及使用谐波频带的参考体模的回波信号功率谱以用于确定第二反向散射系数。
CN202080022296.0A 2019-03-20 2020-02-06 使用基波和谐波信号的定量超声 Pending CN113631953A (zh)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
US16/358,808 US11331078B2 (en) 2019-03-20 2019-03-20 Quantitative ultrasound using fundamental and harmonic signals
US16/358808 2019-03-20
PCT/US2020/016912 WO2020190410A1 (en) 2019-03-20 2020-02-06 Quantitative ultrasound using fundamental and harmonic signals

Publications (1)

Publication Number Publication Date
CN113631953A true CN113631953A (zh) 2021-11-09

Family

ID=69771130

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202080022296.0A Pending CN113631953A (zh) 2019-03-20 2020-02-06 使用基波和谐波信号的定量超声

Country Status (5)

Country Link
US (2) US11331078B2 (zh)
KR (1) KR102628879B1 (zh)
CN (1) CN113631953A (zh)
DE (1) DE112020001333T5 (zh)
WO (1) WO2020190410A1 (zh)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10682098B2 (en) * 2018-03-22 2020-06-16 Shenzhen Mindray Bio-Medical Electronics Co., Ltd. Predictive use of quantitative imaging
US11457896B2 (en) * 2020-09-23 2022-10-04 GE Precision Healthcare LLC Ultrasound imaging system and method for generating an enhanced image to enhance a shadow region
US11700839B2 (en) 2021-09-01 2023-07-18 X. Development Calibration target for ultrasonic removal of ectoparasites from fish
US11864535B2 (en) 2021-12-21 2024-01-09 X Development Llc Mount for a calibration target for ultrasonic removal of ectoparasites from fish

Citations (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4470303A (en) * 1982-09-20 1984-09-11 General Electric Company Quantitative volume backscatter imaging
US4936308A (en) * 1988-05-27 1990-06-26 Agency Of Industrial Science & Technology Method and apparatus for measuring acoustic characteristics and temperature
CN1249164A (zh) * 1998-09-30 2000-04-05 柯坚 人体组织超声衰减成像技术
WO2004093671A2 (en) * 2003-04-22 2004-11-04 Wisconsin Alumni Research Foundation Parametric ultrasound imaging using angular compounding
US20050277835A1 (en) * 2003-05-30 2005-12-15 Angelsen Bjorn A Ultrasound imaging by nonlinear low frequency manipulation of high frequency scattering and propagation properties
WO2007088508A2 (en) * 2006-02-03 2007-08-09 Koninklijke Philips Electronics N.V. Ultrasonic method and apparatus for measuring or detecting flow behavior of a non-sinusoidal periodicity
US20120232388A1 (en) * 2011-03-10 2012-09-13 University Of Washington Through Its Center For Commercialization Ultrasound systems and methods for real-time noninvasive spatial temperature estimation
CA2912791A1 (en) * 2013-05-24 2014-11-27 Sunnybrook Research Institute System and method for classifying and characterizing tissues using first-order and second-order statistics of quantitative ultrasound parametric maps
CN106037796A (zh) * 2015-04-16 2016-10-26 美国西门子医疗解决公司 量化的粘弹性超声成像
CN106841385A (zh) * 2017-01-15 2017-06-13 长沙理工大学 基于声‑超声的聚丙烯生产管道粉末粘附状态的检测方法
CN106846458A (zh) * 2016-12-15 2017-06-13 哈尔滨工业大学（威海） 基于3d打印的超声波立体模型构建方法及装置
US20170273667A1 (en) * 2016-03-22 2017-09-28 Siemens Medical Solutions Usa, Inc. Relative backscatter coefficient in medical diagnostic ultrasound
US20180120272A1 (en) * 2016-11-02 2018-05-03 Wyatt Technology Corporation Method to eliminate periodic noise from data collected with a chromatography system
CN108680928A (zh) * 2018-03-30 2018-10-19 南京牧镭激光科技有限公司 消光系数的测量方法、雷达及存储装置
CN108685596A (zh) * 2017-04-06 2018-10-23 美国西门子医疗解决公司 利用超声医疗成像的组织性质估计

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP6532802B2 (ja) 2014-12-15 2019-06-19 キヤノンメディカルシステムズ株式会社 超音波診断装置

Patent Citations (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4470303A (en) * 1982-09-20 1984-09-11 General Electric Company Quantitative volume backscatter imaging
US4936308A (en) * 1988-05-27 1990-06-26 Agency Of Industrial Science & Technology Method and apparatus for measuring acoustic characteristics and temperature
CN1249164A (zh) * 1998-09-30 2000-04-05 柯坚 人体组织超声衰减成像技术
WO2004093671A2 (en) * 2003-04-22 2004-11-04 Wisconsin Alumni Research Foundation Parametric ultrasound imaging using angular compounding
US20050277835A1 (en) * 2003-05-30 2005-12-15 Angelsen Bjorn A Ultrasound imaging by nonlinear low frequency manipulation of high frequency scattering and propagation properties
WO2007088508A2 (en) * 2006-02-03 2007-08-09 Koninklijke Philips Electronics N.V. Ultrasonic method and apparatus for measuring or detecting flow behavior of a non-sinusoidal periodicity
US20120232388A1 (en) * 2011-03-10 2012-09-13 University Of Washington Through Its Center For Commercialization Ultrasound systems and methods for real-time noninvasive spatial temperature estimation
CA2912791A1 (en) * 2013-05-24 2014-11-27 Sunnybrook Research Institute System and method for classifying and characterizing tissues using first-order and second-order statistics of quantitative ultrasound parametric maps
CN106037796A (zh) * 2015-04-16 2016-10-26 美国西门子医疗解决公司 量化的粘弹性超声成像
US20170273667A1 (en) * 2016-03-22 2017-09-28 Siemens Medical Solutions Usa, Inc. Relative backscatter coefficient in medical diagnostic ultrasound
US20180120272A1 (en) * 2016-11-02 2018-05-03 Wyatt Technology Corporation Method to eliminate periodic noise from data collected with a chromatography system
CN106846458A (zh) * 2016-12-15 2017-06-13 哈尔滨工业大学（威海） 基于3d打印的超声波立体模型构建方法及装置
CN106841385A (zh) * 2017-01-15 2017-06-13 长沙理工大学 基于声‑超声的聚丙烯生产管道粉末粘附状态的检测方法
CN108685596A (zh) * 2017-04-06 2018-10-23 美国西门子医疗解决公司 利用超声医疗成像的组织性质估计
CN108680928A (zh) * 2018-03-30 2018-10-19 南京牧镭激光科技有限公司 消光系数的测量方法、雷达及存储装置

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
KAJOLI BANERJEE KRISHNAN: "A Two-Parameter Model for Ultrasonic Tissue Characterization with Harmonic Imaging", ARXIV MEDICAL PHYSICS, 31 December 2017 (2017-12-31), pages 1 - 26 *
MICHAEL F. INSANA等: "PARAMETRIC ULTRASOUND IMAGING FROM BACKSCATTER COEFFICIENT MEASUREMENTS: IMAGE FORMATION AND INTERPRETATION", 《ULTRASONIC IMAGING》, 11 April 2004 (2004-04-11), pages 245 - 267 *
李硕;杨春兰;盛磊;朱皓;白燕萍;吴水才;: "利用超声射频实现肿瘤热疗实时无创测温方法", 北京生物医学工程, no. 06, 15 December 2009 (2009-12-15) *
李莉, 吴朝霞, 程敬之: "定量背向散射成像研究", 生物医学工程学杂志, no. 01, 25 March 1999 (1999-03-25) *
王玲;涂沛;徐铭恩;: "激光辐照下组织光热响应的实时监测研究", 中国激光, no. 01, 10 January 2015 (2015-01-10) *

Also Published As

Publication number Publication date
US20200297322A1 (en) 2020-09-24
DE112020001333T5 (de) 2021-12-02
US11617567B2 (en) 2023-04-04
US20220240900A1 (en) 2022-08-04
WO2020190410A1 (en) 2020-09-24
KR102628879B1 (ko) 2024-01-23
KR20210137561A (ko) 2021-11-17
US11331078B2 (en) 2022-05-17

Similar Documents

Publication Publication Date Title
KR102628879B1 (ko) 기본 및 고조파 신호들을 사용하는 정량적 초음파
CN107616809B (zh) 医学诊断超声中的组织表征
US8197408B2 (en) Sparse tissue property measurements in medical ultrasound imaging
US20060241456A1 (en) Ultrasonic imaging apparatus and ultrasonic imaging method
CN108652660B (zh) 用于医学诊断超声中的衰减估计的衍射校正
US9168027B2 (en) Adaptive acoustic pressure estimation in medical ultrasound
US5052394A (en) Method and apparatus for ultrasonic beam compensation
US20060079780A1 (en) Ultrasonic imaging apparatus
JP5642910B1 (ja) 超音波観測装置、超音波観測装置の作動方法および超音波観測装置の作動プログラム
US10799208B2 (en) Compressional sound speed imaging using ultrasound
EP3023058B1 (en) Ultrasound observation apparatus, method for operating ultrasound observation apparatus, and program for operating ultrasound observation apparatus
JPS60122367A (ja) 超音波測定方法およびその装置
US20230404537A1 (en) Ultrasound medical imaging with optimized speed of sound based on fat fraction
KR20190087041A (ko) 초음파 영상장치 및 그 제어방법
US7481769B2 (en) Ultrasonic diagnosing apparatus
KR20120067535A (ko) 미드 포인트 알고리즘에 기초하여 ｈｐｒｆ 도플러 영상을 제공하는 초음파 시스템 및 방법
KR102220822B1 (ko) Arfi 이미징을 위한 교정
US8398548B2 (en) Ultrasound diagnostic apparatus and ultrasound diagnostic method
JPH0556466B2 (zh)
JP2001000434A (ja) 物質をイメージングする方法及びイメージング・システム
Ehab et al. Ultrafast Phase Aberration Correction in Ultrasound Imaging Using a Simple Model for fat layer
JPH0390139A (ja) 血流のエコーを基準とした超音波計測装置
Ehab Using a Simple T1Vodel forFat Layer
Le Croissette Signal Processing In Ultrasound
JPH0425017B2 (zh)

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