CN104013419A - 医疗超声波中的自适应声压估计 - Google Patents

医疗超声波中的自适应声压估计 Download PDF

Info

Publication number
CN104013419A
CN104013419A CN201410058929.XA CN201410058929A CN104013419A CN 104013419 A CN104013419 A CN 104013419A CN 201410058929 A CN201410058929 A CN 201410058929A CN 104013419 A CN104013419 A CN 104013419A
Authority
CN
China
Prior art keywords
acoustic pressure
response
function
tissue
frequency
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.)
Granted
Application number
CN201410058929.XA
Other languages
English (en)
Other versions
CN104013419B (zh
Inventor
W.G.维尔克宁
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 CN104013419A publication Critical patent/CN104013419A/zh
Application granted granted Critical
Publication of CN104013419B publication Critical patent/CN104013419B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/58Testing, adjusting or calibrating the diagnostic device
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61NELECTROTHERAPY; MAGNETOTHERAPY; RADIATION THERAPY; ULTRASOUND THERAPY
    • A61N7/00Ultrasound therapy
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/48Diagnostic techniques
    • A61B8/481Diagnostic techniques involving the use of contrast agent, e.g. microbubbles introduced into the bloodstream

Landscapes

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

Abstract

一种医疗超声波中提供的自适应声压估计(40)。测量(30)不同频率下的特定患者的组织的响应(例如,线性和非线性响应)。该响应被用来确定(38)声压。特定患者中的测量使估计适合于患者。确定(38)对于期望位置的声压以便将发射功率设置成期望水平,诸如用于具有高SNR但具有造影剂的有限毁坏的造影剂成像,或者提供期望的热剂量声学治疗。

Description

医疗超声波中的自适应声压估计
技术领域
本实施例涉及超声波的医疗使用。特别地,超声波的声压适合于诊断或治疗需要。
背景技术
超声波造影剂是稳定化微泡。造影剂是压力敏感的。高声压破坏了微泡。为了不破坏微泡,发射功率与典型B模式成像条件相比通常减小10 dB—30 dB。此外,微泡的非线性响应可弱于组织中的散射体的线性响应。因此,强反差图像的信噪比(SNR)是低的,因此期望使用不损坏微泡的最大功率。
为改善SNR或使泡毁坏最小化来调整发射功率与数据获取相干扰,因为只能在微泡已进入成像平面之后执行调整。在可重复较高质量获取之前,微泡已清除血流(例如,约15分钟)。此外,从患者安全角度和从成本角度出发,造影剂的多次注射是不期望的。
为了尝试在造影剂的注射之前设置发射功率,临床医生使用机械指数(MI)值作为声压的指示符。虽然MI与声压有关,但MI并未考虑患者的解剖体且并不提供关于成像平面中的声压分布的信息。有经验的临床医生可基于B模式图像的SNR来进行调整,但是这是臆测。
对于超声波治疗而言,对患者施加期望的热剂量。过多的热剂量可损害健康组织,并且过少的热剂量可能不会充分地治疗患者。然而,各种中间结构可导致意外的热剂量。类似于造影剂成像,预测实际压力的能力是有限的。
其它治疗应用可需要某些声压水平,例如以通过破坏外壳来从微泡释放药物或确立稳态空化(stable cavitation)。
发明内容
通过介绍,下面描述的优选实施例包括用于医疗超声波中的自适应声压估计的方法、系统、计算机可读介质以及指令。测量不同频率下的特定患者的组织的响应。该响应被用来确定声压。特定患者中的测量使估计适于患者。确定用于期望位置的声压以便将发射功率设置成期望水平,诸如用于具有高SNR但具有造影剂的有限毁坏的造影剂成像,或者提供期望的热剂量或声压用于治疗。
在第一方面中,提供了一种用于医疗超声波中的自适应声压估计的方法。针对多个位置在第一频率下获取第一响应。针对所述多个位置在第二频率下获取第二响应。估计位置中的一个或多个处的声压。作为第一和第二响应的函数来估计声压。将声学输出设置为声压的函数。
在第二方面中,非临时计算机可读存储介质具有存储在其中的数据,其表示可由已编程处理器执行用于医疗超声波中的自适应声压估计的指令。该存储介质包括用于确定对声能的第一和第二频带中的组织响应、作为第一和第二频带中的组织响应的函数来设置发射功率并使用发射功率来执行造影剂成像或声学治疗的指令。
在第三方面中,提供了一种用于医疗超声波中的自适应声压估计的超声波系统。发射放大器与换能器相连。处理器被配置成作为第一频率下的作为深度的函数的第一振幅轮廓(amplitude profile)和不同于第一频率的第二频率下的作为深度的函数的第二振幅轮廓的函数来确立发射放大器的发射功率。
下面结合优选实施例来讨论本发明的更多方面和优点。本发明由下面的权利要求限定,并且不应将本部分中的任何东西理解为对那些权利要求的限制。
附图说明
部件和图不一定按比例,而是将重点放在说明实施例的原理上。此外,在图中,遍及不同的视图,相似的参考数字指示对应的部分。
图1是用于医疗超声波中的自适应声压估计的方法的一个实施例的流程图;
图2是处于不同接收频率和发射功率水平的示例性振幅轮廓的图形表示;
图3是处于不同接收频率和发射功率水平但具有与图2不同的发射焦点深度的其它示例性振幅轮廓的图形表示;以及
图4是用于医疗超声波中的自适应声压估计的系统的一个实施例的框图。
具体实施方式
超声波对比成像需要声学输出功率的谨慎调整以避免毁坏试剂并使灵敏度最大化。低于最佳的声压可导致降低的灵敏度和/或穿透度。高于最佳的声压可导致造影剂的过多毁坏。类似地,超声波治疗受益于声能的精确设置以提供期望的剂量或压力。
生物组织中的超声波传播是非线性的,导致接收回波信号中的多个谱分量。通过分析对发射声压具有不同相关性(诸如线性的和二次的)的谱分量的振幅,可估计和调整目标器官中的声压以优化对比成像或声学治疗。
例如,可通过测量对于接收回波的基波和谐波(例如,二次)分量的强度轮廓(强度对比深度)来估计身体内部的声压。在成像平面或体积中的多个位置处确定基波和谐波响应。此确定可在造影剂到达所关心区域处之前小于一分钟或几分钟内发生,因此几乎不增加超声波检验的时间。根据响应来估计用于至少一个位置的绝对声压。响应于该压力估计而调整或者另外设置声学输出以实现最佳的强反差图像质量。
组织可以是‘暗的’,因为不存在到达那里的太多声音,或者因为组织并未散射太多。非线性传播并未过多地受到散射的影响。这就是为什么可从‘暗的’组织得到与基波相比强大的二次谐波的原因,表明声压实际上是高的。通过准确地预测压力,确定了关于组织的回声反射性的更多东西。
图1图示出用于医疗超声波中的自适应声压估计的方法的一个实施例。使用不同频率下的对组织的超声波的响应来估计声压。基于此声压来设置用于后续发射的发射功率。对声压的适应可在没有声压的直接估计的情况下基于不同频率下的响应来设置发射功率。
针对任何超声波发射来设置发射功率。例如,针对治疗来设置发射功率以更可能地提供期望的剂量。作为另一示例,针对造影剂成像来设置发射功率以在最小化或避免期望区域中的造影剂的毁坏的同时改善SNR。
用图4的系统来实施图1的方法,但可使用其它系统。医疗诊断成像系统执行该方法。替换地,超声波治疗系统执行该方法。在又一实施例中,与有或没有额外波束成形器电子装置的换能器相连的计算机执行该方法。
以所示的顺序或不同顺序执行该方法。例如,以任何顺序或同时地执行动作30、34和36。可以提供附加、不同或较少动作。例如,不执行动作34和36中的一者或两者。作为另一示例,不执行动作42。可不直接地执行动作40,诸如其中在没有实际声压的特定估计的情况下指示声压的信息被用来设置发射功率水平。在又一示例中,执行用于收集更多测量或值以估计声压的一个或多个动作和/或配置动作。
自动地或手动地发起发射功率设置特征。例如,用户按下按钮,按下脚踏板,选择菜单项或说话(例如通过语音激活的控制)。作为另一示例,用户选择造影剂成像软件或造影剂成像应用所特定的软件。作为又一示例,由连续活动软件来生成该请求。在该请求之前或之后将造影剂注射到患者中。响应于该请求,任何当前成像被中断或被交错有数据的获取用以设置发射功率水平。
在动作30中,确定组织对声能的响应。针对不同的频率确定响应。不同的频率是不同的中心频率和/或不同的频带。该频带可重叠或者可不重叠。例如二次谐波回波振幅与基波回波振幅之间的比提供组织内部的绝对基波声压的指示。
通过在动作32中获取处于不同频率的响应来执行动作30。在动作32中,获取超声波回波数据。例如,针对一个或多个位置执行用以满足检测技术的发射和接收(例如用于每个基波数据的多个脉冲、用于每个谐波数据的多个脉冲或用于B模式的单个脉冲)。在一个实施例中,使用四个或更多脉冲来从给四个或更多回波导出不包括二次谐波的基波和不包括一次谐波的二次谐波 。
位置可沿着扫描或波束线和/或分布在平面或体积中的不同扫描线之中。针对相同的扫描线,可使用相同的发射和对应的接收操作来获得不同位置处的响应。针对沿着分离的扫描线的位置,使用对不同扫描线的连续发射和接收操作。可替换地,执行宽波束或多波束发射和对应的多波束接收操作。
该位置可处于用于成像的采样密度或此类位置的子集。较少的位置可需要较少的处理和/或扫描时间,因此在一个实施例中,线密度是用于B模式组织成像的线密度的1/4至1倍。可以使用其它线或样本密度,诸如沿着用于平面或体积的一个扫描线执行测量。针对每个给定线,可使用任何采样密度。例如,使用与被用于B模式成像的相同的采样密度。获取沿着扫描线的来自各种深度的回波信号。
被用于采样的一个或多个扫描线和沿着该线的样本位置被选作具有已知反向散射特性、诸如充分发展斑点的区域。在肝脏中,大部分器官是软组织,因此具有已知的反向散射特性。该位置被选则作为软组织位置而不是大的血管或静脉位置。在替换实施例中,该位置与其它类型的组织或结构相关联。
也可考虑传播路径。使用仅传播通过显示出充分发展斑点的组织的波束或扫描线。不使用横断或通过血管、骨、硬组织或流体组织而到达采样位置处的波束。替换地,使用通过其它类型的组织的波束。在许多情况下,近场中的组织无关紧要,只要近场不包括骨即可且只要大体上或平均起来与正在测量的组织相同即可。沿着‘均匀的’一定范围的波束取得样本。
手动地或自动地选择该位置和/或扫描线。例如,用户指示将在B模式或B模式和流模式组合图像上使用的波束和/或样本位置。作为另一示例,针对遍及平面和/或体积的位置确定信噪比(SNR)。斑点与SNR中的(空间)方差相关联。在平均方差与标准偏差的比为约1.9的情况下,指示与软组织相关联的斑点区域。
在两个或更多不同频率下获取用于各位置的响应。在本文中所使用的一个示例中,使用基波或发射频率(例如,2 MHz)和基波频率的二次谐波(例如,4 MHz)。基波响应充当回声反射性的度量。
替代基波(一次谐波)和二次谐波响应,可使用两个或更多谐波的其它组合来估计绝对声压。可使用与非线性波方程式中的不同项有关的信号分量的任何组合,诸如二次谐波和三次基波。
为了将处于不同频率的响应隔离,处理表示声学回波的接收信号。使用相同的数据来确定不同频率下的响应。例如,用滤波器来将回波信号分解以便确定反射谐波的振幅。未滤波信号被用于基波,但可使用滤波。作为另一示例,采用已编码发射脉冲序列(例如,脉冲反相)。针对脉冲反相,两个另外相同的发射脉冲的振幅是[1, -1]。通过减去接收信号,可抵消来自偶次谐波的贡献,在基波频率下提供信息。其它奇次谐波下的信息是相对小的,因此可保留,但是在其它实施例中可被滤出。通过加上接收信号,抵消了来自奇次谐波的贡献。类似地,除二次谐波之外的其它偶次谐波下的信息是小的,因此可保留,但是可被滤出。任何附加滤波可改善谐波分量之间的分离并拒绝噪声。可使用具有不同相位、振幅或相位和振幅后面是在接收时的组合(例如,加权和滤波)的三个或更多发射脉冲的序列来将处于期望波段的信息隔离。
在替换实施例中,使用不同的发射和接收事件来获取不同的响应。例如,执行发射脉冲和对应的接收操作以便获得用于基波频率的各位置处的响应。执行分离或连续的发射和接收操作以便获得用于不同基波或谐波频率的各位置处的响应。
该响应是作为位置的函数的振幅。例如,该响应是作为深度的函数的振幅。使用检测之前或之后的波束成形样本。例如,使用强度或B模式检测之后的波束成形样本。可去除由波束成形器施加的任何深度增益修正或补偿。替换地,并未首先提供深度增益修正。在又一实施例中,一致地施加深度增益修正,并且在不去除增益的情况下使用所得到的信号。
用相同的发射焦点位置来获取作为深度的函数的振幅。使用动态接收聚焦,但可提供单个接收焦点。图2示出了在20mm处具有发射焦点的示例性响应。可使用其它发射焦点深度,诸如图3中所示的40 mm焦点。
针对两个或更多频带提供作为深度的函数的振幅。图2和3示出了用于基波和二次谐波频率的响应。
为了改善估计或设置的准确度,可在连续发射和接收事件中获取不同发射焦点处的响应。例如,获取图2和3两者的响应。改变发射焦点深度以收集不同深度处的测量以确认关于不同深度处的B/A(非线性参数,其中B和A是描述声音速度对局部压力的相关性的泰勒级数的前两个系数)、衰减以及绝对声压的假设。
响应是针对处于给定发射功率的发射。在替换实施例中,针对不同发射功率中的每一个获取不同频率下的响应。针对其中每个不同发射水平确定作为位置(例如,深度)的函数的振幅。图2和3示出了用于每个频率的两个响应的示例,其中,所述两个响应对应于两个发射功率。一个发射功率从另一个向下-3 dB,但是可使用更小或更大的差。可使用超过两个不同的发射功率水平。连续地获取用于不同发射功率水平的响应,诸如响应于不同发射事件。
图2和3图示出用于估计声压的某些现象。当将发射功率减小3 dB时,二次谐波的信号水平对于与二次谐波相比处于基波的很好地分离的响应而言下降6 dB。在没有脉冲反相和窄脉冲的情况下,当发射信号水平由于基波漏泄而减小3 dB时,二次谐波可仅下降4 dB。在一定深度处,平方的行为可能不再适用,因为二次谐波下降接近于本底噪声。通过在两个功率水平分析两个曲线,系统可识别SNR极限。
第二谐波水平快速地增加接近于焦点深度。二次谐波与基波之间的比(以dB为单位的差)随着发射功率的增加而增加。用于焦点深度之外的基波信号和二次谐波信号的斜率由于波束样式和衰减方面的差异而不同。图2和3的示例具有作为深度的函数的组织模拟材料的回波强度,但是在其它实施例中使用实际患者的响应。
在动作34中,确定衰减。可使用衰减来更准确地确定不同频率下的响应的比。可将该衰减用于估计声压。
该衰减可以是假设或预定值。例如,衰减可以是基于组织类型的值。用户指示组织的类型,诸如选择肝脏成像应用。可根据超声波信号来确定组织的类型并被用来选择预定衰减值。可使用用于通过各种类型的组织的波束路径的衰减。替换地,衰减是与一般软组织而不非特定类型的组织相关联的预定值。
在其它实施例中,根据获取的响应或其它超声波测量来确定衰减。例如,根据不同频率下的响应的斜率来确定衰减。确定最大值之前和之后的基波和谐波信号的斜率。振幅的最大值处于或接近于发射焦点。可忽略任何近场最大值。
在组织中,衰减与频率近似地成比例。谐波遭受较高衰减。超声波成像通常使用聚焦的波束。如果f数是低的,基波的压力朝着焦点增加并达到接近于焦点或在焦点处的最大值。随着基波振幅的增加,谐波的振幅由于该现象的非线性和累积性质而快速地增加。在已达到最大值之后,谐波的振幅快速地减小,因为来自基波的转换减弱,并且因为谐波经受较高的衰减。在沿着超声波波束的给定位置处,基波的振幅与发射波成比例,忽略了基波耗尽,并且二次谐波的振幅与求平方的发射信号的振幅近似地成比例。
给定高斯发射脉冲,响应的斜率的比指示衰减。该斜率来自远场区,诸如到最大值之后的响应的拟合线。可使用相对浅的发射焦点,诸如20mm焦点。在其它实施例中可使用处于最大值或其它离散深度的近场或比。可计算不同深度处的衰减。可使用衰减的其它计算,诸如观察中心频率对比深度的降档(downshift)。可使用两个或多个不同发射频率下的测量。
在动作36中,确定用于非线性参数的值。生物组织或流体中的波传播是非线性的。非线性传播引起波形随着波形传播通过介质而失真。该失真随着波行进的距离而增加。由于失真是非线性效应,失真也随着振幅而增加。
还可将失真描述为具有作为未失真波的频率(即,基波频率)的整数倍的频率的谐波的生成。谐波的能量由于能量守恒而来自基波。然而,诊断超声波中的条件使得谐波比基波弱得多,并且使得基波的振幅中的损耗(即,基波耗尽)可能是不可测量的。
基波和谐波被介质之间的界面反射或者被散射体散射并返回朝着换能器传播,形成回波信号。这些回波一般地足够弱,使得非线性传播是可忽略的。散射是频率相关的。然而,散射并不取决于振幅(即,散射是线性的),并且可由超声波测量(例如用具有不同频率的两个发射信号进行的脉冲回波测量)来确定对于某个类型的组织的频率相关性。
可将非线性性表示为作为局部压力的函数的声音速度。该非线性性常常以比B/A来表征,其中,B和A是描述声音速度对局部压力的相关性的泰勒级数的前两个系数。B/A取决于组织类型,并且谐波的振幅取决于B/A。使用脉冲回波技术来测量B/A。使用了响应或其它超声波数据。替换地,在文献中找到用于某些器官的值。由用户来选择或基于所检测的组织类型来使用预定非线性参数值。该值可基于正在被扫描的组织的类型。
在动作38中,将发射功率设置为不同频带中的组织响应的函数。也可基于处于不同功率和/或发射焦点的响应来设置发射功率。还可将衰减和/或非线性参数值用于设置发射功率。可使用附加、不同或较少的输入。
基于声压来设置发射功率。可使用指示声学功率的输入(诸如不同频率下的响应)来设置发射功率。替换地且如动作40所反映的,可估计声压,并且使用压力的估计来设置发射功率。作为不同频率下的组织响应、衰减以及超声波非线性参数的函数来估计声压。
采取另一示例,替代响应或除响应之外,可使用响应的特性。可计算作为位置的函数的振幅的斜率、平均值、差或其它特性。可针对每个响应来计算特性,或者特性可以是相对的(例如,表示响应之间的类似性或差异)。作为另一示例,使用发射波束成形参数。可使焦点深度、切趾特性、孔径尺寸、发射和接收时的元件接受角、发射和/或接收波束对阵列的角度和/或表示波束成形的其它参数与压力的量相关联。
针对对于其而言响应被采样的位置中的至少一个来估计声压。例如,针对发射焦点来估计声压,但是可使用较浅或较深的位置。可针对多个位置估计声压。
输入被提供给模型。该模型响应于输入而输出声压。该模型是正向模型,但是可使用机器学习模型。例如,通过使模拟结果与测量响应匹配而使用实时模拟来搜索声压B/A和衰减值。
在另一实施例中,由查找表或数据库来表示模型。基于不同频率下的响应、从响应导出的参数值、衰减值和/或非线性参数值来查找声压。与不同频率相关联的两个或更多输入(例如,不同频率下的响应)被用于查找输出估计。通过将振幅轮廓或从此类轮廓导出的参数与查找表或数据库相比较,可将绝对声压确定为深度或其它位置分布的函数。
数据库或查找表中的值是实验地确定和/或模拟的。例如,通过使用KZK方程式或具有非线性传播附加的阵列模拟器(例如,Field II)的模拟来生成查找表。
模拟或实验地确定的查找表可包括用于非线性参数和/或衰减的输入。替换地,为不同的非线性参数和/或衰减值提供分离的表格。可使用表格的组合或组合表格的任何组合。
在又一实施例中,模型表示用于给定发射配置和组织的不同频率下的预期响应。该响应是在不同发射功率下获取的。改变发射功率直至响应最佳地或充分地与预期响应匹配为止。与该匹配相关联的发射功率提供声压与发射水平之间的预期关系。
该模型输出声压和/或发射功率。基于不同位置和/或相同位置处的响应来估计一位置处的声压。
由于衰减,不同位置处的声压是不同的。可由模型来估计不同位置处的声压。例如,将用于该位置的输入提供给模型。替换地,基于另一估计声压来确定其它位置处的声压。给定衰减值,可使用一个位置处的声压来确定对于其它位置的压力。该衰减为估计由相同的发射事件引起的压力提供比例因子。
在动作38中,基于声压来设置发射功率。一旦确定了声压,就设置用于造影剂医疗诊断成像、治疗或其它使用的发射功率。查找表、函数或其它模型使声压与发射功率水平相关。该设置可以是绝对值或变化量。该设置可发生一次,或者可以是其中使用不同的发射水平直至得到期望压力为止的迭代过程。
在一个示例中,设置用于发射操作的声学输出以在尝试具有基本上最大输出的同时避免造影剂的毁坏。基本上解决了偏差或公差以为误差提供裕度,诸如10%。由于造影剂的变化和/或压力估计中的不准确,提供公差以避免毁坏所关心区域中的大多数造影剂。通过估计所关心区域处的声压,可提供导致非毁坏性且良好的SNR压力的发射功率。该估计适于患者、系统和/或波束路径变化。
在另一示例中,设置用于发射操作的声学输出以提供期望的剂量。该热剂量对应于在所关心区域处施加的声压。通过估计所关心区域处的声压,该发射功率导致不太可能损害健康组织但仍能够充分地治疗不健康组织的期望热剂量。该估计适于患者、系统和/或波束路径变化。
在一个实施例中,除被估计的压力之外,还基于输入来设置发射功率。为了在不分离地估计压力的情况下直接设置发射功率,可在模型中包括附加输入以输出发射水平。否则,执行进一步计算或查找功能。该设置解决了用来生成压力和发射水平与压力的关系的发射水平。可使用任何附加输入。
一个输入是噪声。可测量该噪声。例如,在没有任何发射事件的情况下测量接收信号。可使用任何噪声测量。
另一输入是造影剂的类型、给定时间内的期望剂量或为其设置了水平的发射的其它使用。例如,可以声能或压力的不同阈值量来破坏不同的造影剂。与毁坏相关联的水平是相对的,因为微泡的变化可导致至多任何水平的某些毁坏。在其处10%、1/3或其它量的造影剂被破坏的压力和对应的发射水平对不同类型的造影剂而言可以是不同的。通过指示造影剂的类型,指示了对应的期望压力水平。作为另一示例,被施加的热剂量可取决于治疗的类型、要治疗的区域的尺寸以及任何散热器(例如,血管)的接近性。所得到的压力和对应的发射水平对于不同的治疗状况而言可以是不同的。
基于成像平面中的声压分布的分析、造影剂的类型或特性以及本底噪声测量,选择发射功率设置。可以由临床医生来设置和/或基于应用而预定附加优化标准。例如,可基于压力来设置功率水平以满足特定标准,诸如最佳SNR、最小毁坏或最大穿透度。
另外,该算法可优化发射脉冲的其它参数,诸如中心频率和/或带宽。可优化接收处理参数,诸如接收带宽和/或增益,因为可基于测量的参数来预测接收频谱和SNR或对其进行建模。知道声压和衰减,还可修改发射聚焦。
在动作42中,执行成像或治疗。成像或治疗包括发射声能。该声能是以在动作38中设置的功率水平发射的。例如,为成像肝脏中造影剂执行一次或多次发射。使用动脉相和门静脉相中的增强样式来诊断病变。基于对处于设置功率水平的发射的响应来检测或观察该样式。向用户显示造影剂的一个或多个图像。在一个实施例中,在显示器上向用户指示机械指数或功率水平。
所使用的发射功率水平被设置成针对给定情况来提供期望的压力水平。发射事件的信号振幅、孔径、f数、切趾法或其它特性被配置成提供发射功率水平。
图4示出了用于医疗超声波中的自适应声压估计的超声波系统。该系统实施了图1的方法或另一方法。该系统包括换能器12、发射放大器14、接收波束成形器16、图像处理器18、处理器22、显示器20以及存储器22。可提供附加、不同或较少的部件,诸如与B模式和流成像相关联的多个检测器。在一个实施例中,该系统是医疗诊断超声波系统。在另一实施例中,该系统是治疗超声波系统。在其它实施例中,该系统是计算机或服务器。
换能器12是压电材料的单个元件或多个元件。在替换实施例中,换能器12包括电容性薄膜结构。对于多个元件而言,换能器12是线性、弯曲线性或多维阵列。可以使用用于在电能与声能之间进行转换的其它换能器。换能器12响应于电波形而生成声能。在由发射放大器14设置的功率下输出声学波形。换能器12响应于声学回波而生成接收信号,诸如用于测量对发射声能的响应。
发射放大器14与换能器12相连,并且是可变放大器、数模转换器或用于改变或增加发射波形的功率、峰值电压或其它功率特性的其它模拟或数字器件。在替换实施例中,发射放大器14包括分压器或用于减少与发射波形相关联的功率的其它器件。为每个系统通道或换能器元件提供了分离的发射放大器14,但是可将一个发射放大器14用于多个通道或元件。在一个实施例中,发射放大器14应用切趾法以便沿着波束进行发射,并被包括作为发射波束成形器的一部分。
从发射放大器14输出的发射波形被换能器12转换成声能。响应于声能和任何造影剂的回波信号被换能器12接收到。换能器12将回波信号转换成电信号或数据。
该回波信号由接收波束成形器16进行波束成形。相对延迟或相位和切趾法被应用于用于接收孔径的信号。随时间推移或者针对不同的深度,相对延迟和/或定相改变,提供动态接收聚焦。将信号组合,诸如求和以进行波束成形。所得到的样本或数据表示沿着线的不同深度处的位置。数据包括一个或多个数字样本或模拟信息。
为了隔离或辨别不同频率下的信息,接收波束成形器16包括滤波器。使用存储器,可将相同的数据滤波多次。在其它实施例中提供了并行滤波。滤波器可以是可编程滤波器或固定滤波器。例如,提供高通和低通滤波器以便将基波和谐波信息分离。作为另一示例,使用一个或多个带通滤波器或串联低通和高通滤波器。替换地,滤波器包括有或没有乘法器的存储器或缓冲器以便对与发射时的不同相位、振幅和/或极性相关联的接收数据求和或求差。执行脉冲反相、振幅编码或其它频率隔离组合。可使用用于识别处于基波发射频带(f)的包括分谐波(例如1/2f)、分数谐波(3/2f)和/或整次谐波(例如2f)的谐波的信息的其它技术和对应的滤波器。
在接收波束成形或其它接收技术之后,数据被提供给图像处理器18或处理器22。在用于估计声压的响应使用强度值的情况下,在估计之前检测数据。替换地,处理器22在检测之前基于波束成形数据来估计压力。
图像处理器18是造影剂检测器、B模式检测器、多普勒检测器、流检测器或用于检测返回信号的一个或多个特性的其它检测器。图像处理器18可以是被用于成像造影剂的相同的检测器。替换地,将图像处理器18与处理器22组合并用于设置发射功率,并且其与检测器分离以用于成像。
处理器22是专用集成电路、通用处理器、数字信号处理器、控制处理器、模拟电路、数字电路、现场可编程门阵列、图形处理单元、其组合或可操作用于设置发射放大器14的发射功率的其它器件中的一个或多个。处理器22控制发射波束成形器14、接收波束成形器16和/或图像处理器18的操作以获得不同频率下的响应和用来估计压力和/或设置发射水平的其它信息。处理器22由软件、硬件、设计或其组合配置。
处理器22被配置成作为两个或更多频率下的振幅轮廓的函数来确立发射放大器14的发射功率,所述振幅轮廓为深度的函数。在其它实施例中,也可使用与不同发射焦点和/或不同发射功率水平相关联的轮廓。使用存储衰减、测量衰减、存储非线性参数值、和/或测量非线性参数值,处理器22通过模型的应用而用振幅轮廓来估计声压。针对一个或多个位置估计声压。可对其它位置处的压力进行估计、外推、内插和/或计算。
使用优化程序或函数,基于声压来设置供在更多发射中使用的发射功率水平。处理器22被配置成设置发射功率水平。可将其它信息用于设置该水平,诸如噪声和应用信息。
存储器24存储查找表、数据库或其它模型信息。存储了输出压力、发射水平、发射特性(例如,孔径、切趾法、焦点区、延迟轮廓和/或f数)或其它输出。基于存储在存储器24中的输入值来查找输出值。发射功率水平和/或压力估计可在有或没有其它信息的情况下基于输入振幅轮廓来选择。替换地,存储器24存储用于基于输入来计算压力或发射水平的算法或其它函数。
在非临时计算机可读存储介质或存储器24(诸如缓存器、缓冲器、RAM、可移动介质、硬盘驱动器或其它计算机可读存储介质)上提供用于实施上面所讨论的过程、方法和/或技术的指令。计算机可读存储介质包括各种类型的易失性和非易失性存储介质。在图中示出或在本文中描述的功能、动作或任务是响应于存储在计算机可读存储介质中或上的一个或多个指令集而执行的。该功能、动作或任务独立于特定类型的指令集、存储介质、处理器或处理策略,并且可由单独地或组合地操作的软件、硬件、集成电路、固件、微代码等来执行。同样地,处理策略可以包括多处理、多任务、并行处理等。
在一个实施例中,指令被存储在可移动介质器件上以供本地或远程系统读取。在其它实施例中,指令被存储在远程位置中以便通过计算机网络或通过电话线传输。在又一实施例中,指令被存储在给定计算机、CPU、GPU或系统内。
显示器20是监视器、液晶显示器、等离子体屏幕、投影仪、触摸屏、平板屏幕、打印机或用于向用户输出信息的其它装置。显示器20可生成图像,诸如对超声波的造影剂响应的图像。该图像可包括估计压力或设置发射水平。
虽然上面已通过参考各种实施例描述了本发明,但应理解的是在不脱离本发明的范围的情况下可以作出许多改变和修改。因此本文旨在应将前述详细描述理解为本发明的当前优选实施例的说明并且不作为对本发明的限定。本文旨在限定本发明的范围的仅仅是包括所有等同物的以下权利要求。

Claims (20)

1.一种用于医疗超声波中的自适应声压估计的方法,所述方法包括:
获取(32)用于多个位置的第一频率下的组织的第一响应;
获取(32)用于所述多个位置的第二频率下的组织的第二响应;
由处理器(22)来估计(40)位置中的至少一个处的声压,所述声压是作为第一和第二响应的函数而估计的;以及
将换能器(12)的声学输出设置(38)为声压的函数。
2.权利要求1的方法,其中,获取(32)第一和第二响应包括分别地在基波和二次谐波频率下的获取(32)。
3.权利要求1的方法,其中,获取(32)第一和第二响应包括沿着线的获取(32),所述位置包括沿着所述线的采样位置。
4.权利要求1的方法,其中,获取(32)第一和第二响应每个包括用相同的发射焦点位置作为深度的函数来获取(32)振幅。
5.权利要求1的方法,其中,获取(32)第一响应包括响应于不同发射功率的获取(32),并且其中,获取(32)第二响应包括响应于不同发射功率的获取(32)。
6.权利要求1的方法,其中,估计(40)包括作为衰减的函数的估计(40)。
7.权利要求6的方法,进一步包括:
根据所述第一响应的斜率来确定(34)衰减。
8.权利要求1的方法,其中,估计(40)包括作为用于各位置的组织的超声波非线性参数的函数的估计(40)。
9.权利要求1的方法,其中,估计(40)包括将第一和第二响应输入到模型中并从模型输出声压。
10.权利要求1的方法,其中,估计(40)包括基于第一响应、第二响应、从第一响应导出的第一参数值或从第二响应导出的第二参数值中的两个或更多的组合来查找声压。
11.权利要求1的方法,其中,估计(40)包括估计(40)所述多个中的第一位置处的声压;
进一步包括:
确定除所述多个中的第一位置之外的其它位置处的声压,所述确定是第一位置处的声压的函数。
12.权利要求1的方法,其中,设置(38)包括设置(38)声学输出以在具有基本上最大输出的同时避免造影剂的毁坏。
13.权利要求1的方法,其中,设置(38)包括作为声压和热剂量的函数来设置(38)声学输出。
14.权利要求1的方法,其中,设置(38)包括:
测量噪声;以及
将发射功率设置(38)为声压和噪声的函数。
15.在具有存储在其中的数据的非临时计算机可读存储介质中,所述数据表示可由已编程处理器(22)执行以用于医疗超声波中的自适应声压估计的指令,所述存储介质包括用于以下各项的指令:
确定(30)对声能的第一和第二频带中的组织响应;
将发射功率设置(38)为第一和第二频带中的组织响应的函数;以及
使用发射功率来执行(42)造影剂成像或声学治疗。
16.权利要求15的非临时计算机可读存储介质,其中,确定(30)组织响应包括针对第一和第二频带中的每一个响应于不同的发射水平而确定(30)不同深度处的组织响应。
17.权利要求15的非临时计算机可读存储介质,其中,设置(38)包括:
作为组织响应、衰减以及超声波非线性参数的函数来估计(40)声压;以及
将发射功率设置(38)为声压的函数。
18.一种用于医疗超声波中的自适应声压估计的超声波系统,所述系统包括:
换能器(12);
发射放大器(14),其与所述换能器(12)相连;以及
处理器(22),其被配置成作为第一频率下的作为深度的函数的第一振幅轮廓和不同于第一频率的第二频率下的作为深度的函数的第二振幅轮廓的函数来确立发射放大器(14)的发射功率。
19.权利要求18的系统,进一步包括存储器(24),存储器(24)存储发射功率,所述发射功率包括可选择作为第一和第二振幅轮廓的函数或作为从第一和第二振幅轮廓导出的参数值的函数的已确立发射功率。
20.权利要求18的系统,其中,所述处理器(22)被配置成根据声压来确立发射功率,所述声压是根据第一和第二振幅轮廓而估计的。
CN201410058929.XA 2013-02-22 2014-02-21 医疗超声波中的自适应声压估计 Active CN104013419B (zh)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
US13/774,885 US9168027B2 (en) 2013-02-22 2013-02-22 Adaptive acoustic pressure estimation in medical ultrasound
US13/774885 2013-02-22
US13/774,885 2013-02-22

Publications (2)

Publication Number Publication Date
CN104013419A true CN104013419A (zh) 2014-09-03
CN104013419B CN104013419B (zh) 2019-06-28

Family

ID=51349550

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410058929.XA Active CN104013419B (zh) 2013-02-22 2014-02-21 医疗超声波中的自适应声压估计

Country Status (5)

Country Link
US (1) US9168027B2 (zh)
JP (1) JP6272707B2 (zh)
KR (1) KR102223052B1 (zh)
CN (1) CN104013419B (zh)
DE (1) DE102014001994B4 (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107205662A (zh) * 2015-02-20 2017-09-26 美国西门子医疗解决公司 医学超声成像中的基于谐波与基波关系的发射功率
CN113117257A (zh) * 2019-12-30 2021-07-16 重庆融海超声医学工程研究中心有限公司 用于检测空化效应的方法及装置、超声治疗设备

Families Citing this family (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10130342B2 (en) 2007-12-28 2018-11-20 Bracco Suisse Sa Initialization of fitting parameters for perfusion assessment based on bolus administration
CA2937467C (en) * 2014-04-07 2022-05-03 Bracco Suisse Sa Estimation of acoustic level in-situ with non-fundamental analysis
GB201513024D0 (en) 2015-07-23 2015-09-09 Univ Heriot Watt A method of, and apparatus for, determination of position in ultrasound imaging
JP6635766B2 (ja) * 2015-11-26 2020-01-29 キヤノンメディカルシステムズ株式会社 超音波診断装置、信号処理装置及び解析プログラム
KR20180091025A (ko) 2015-12-10 2018-08-14 브라코 스위스 에스.에이. 동적 역치화에 의한 고정된 조영제의 검출
WO2017198878A1 (de) * 2016-05-20 2017-11-23 Universität des Saarlandes Automatisierte empfehlung zur gabe von echokontrastmittel mittels machine-learning-algorithmen
WO2019234495A2 (en) * 2018-06-06 2019-12-12 Insightec, Ltd. Improved reflection autofocusing
JP6865810B1 (ja) * 2019-12-25 2021-04-28 ゼネラル・エレクトリック・カンパニイ 超音波診断装置及びその制御プログラム
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
KR102525784B1 (ko) 2020-10-14 2023-04-25 동명대학교산학협력단 휴대용 음압기가 구비된 음압복
CN116251305B (zh) * 2023-05-10 2023-09-01 深圳半岛医疗有限公司 超声换能单元输出功率控制方法、设备和可读存储介质

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1454570A (zh) * 2002-02-08 2003-11-12 株式会社东芝 超声波诊断装置
US20060030779A1 (en) * 2004-08-09 2006-02-09 Siemens Medical Solutions Usa, Inc. Intensity level control for contrast agent imaging
CN101277650A (zh) * 2005-08-05 2008-10-01 美国西门子医疗解决公司 利用医学超声成像的造影剂操控
US20100036292A1 (en) * 2008-08-06 2010-02-11 Mirabilis Medica Inc. Optimization and feedback control of hifu power deposition through the analysis of detected signal characteristics
US20130023767A1 (en) * 2011-05-12 2013-01-24 Mammone Richard J Low-cost, high fidelity ultrasound system

Family Cites Families (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH02193651A (ja) * 1989-01-20 1990-07-31 Fujitsu Ltd 超音波減衰係数の推定方法
JP2735266B2 (ja) * 1989-02-15 1998-04-02 株式会社日立製作所 パルスエコー型超音波装置
JP3862838B2 (ja) * 1997-11-26 2006-12-27 株式会社東芝 超音波診断装置
US6413218B1 (en) 2000-02-10 2002-07-02 Acuson Corporation Medical diagnostic ultrasound imaging system and method for determining an acoustic output parameter of a transmitted ultrasonic beam
JP3538414B2 (ja) * 2001-12-10 2004-06-14 ジーイー横河メディカルシステム株式会社 超音波診断装置
US8021303B2 (en) * 2003-06-12 2011-09-20 Bracco Research Sa System for extracting morphological information through a perfusion assessment process
US7591788B2 (en) * 2003-08-19 2009-09-22 Siemens Medical Solutions Usa, Inc. Adaptive contrast agent medical imaging
US7367943B2 (en) * 2004-01-06 2008-05-06 General Electric Company Systems and methods for determining the electrical activity of a material
US8043219B2 (en) 2004-09-17 2011-10-25 Siemens Medical Solutions Usa, Inc. Automated power level for contrast agent imaging
JP5308674B2 (ja) * 2004-12-23 2013-10-09 ブラッコ・シュイス・ソシエテ・アノニム ボーラス投与に基づく灌流評価法およびシステム
CN101179998B (zh) * 2005-05-20 2010-12-08 株式会社日立医药 图像诊断装置
JP5213019B2 (ja) * 2007-12-14 2013-06-19 ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー 超音波撮像装置
US20100016719A1 (en) 2008-07-16 2010-01-21 Siemens Medical Solutions Usa, Inc. Adaptive regulation of acoustic output power in medical ultrasound imaging
WO2010017419A2 (en) 2008-08-06 2010-02-11 Mirabilis Medica Inc. Optimization and feedback control of hifu power deposition through the analysis of detected signal characteristics
JP5765332B2 (ja) * 2010-03-16 2015-08-19 コニカミノルタ株式会社 超音波診断装置

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1454570A (zh) * 2002-02-08 2003-11-12 株式会社东芝 超声波诊断装置
US20060030779A1 (en) * 2004-08-09 2006-02-09 Siemens Medical Solutions Usa, Inc. Intensity level control for contrast agent imaging
CN101277650A (zh) * 2005-08-05 2008-10-01 美国西门子医疗解决公司 利用医学超声成像的造影剂操控
US20100036292A1 (en) * 2008-08-06 2010-02-11 Mirabilis Medica Inc. Optimization and feedback control of hifu power deposition through the analysis of detected signal characteristics
US20130023767A1 (en) * 2011-05-12 2013-01-24 Mammone Richard J Low-cost, high fidelity ultrasound system

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107205662A (zh) * 2015-02-20 2017-09-26 美国西门子医疗解决公司 医学超声成像中的基于谐波与基波关系的发射功率
CN107205662B (zh) * 2015-02-20 2021-06-15 美国西门子医疗解决公司 医学超声成像中的基于谐波与基波关系的发射功率
CN113117257A (zh) * 2019-12-30 2021-07-16 重庆融海超声医学工程研究中心有限公司 用于检测空化效应的方法及装置、超声治疗设备
CN113117257B (zh) * 2019-12-30 2023-07-11 重庆融海超声医学工程研究中心有限公司 用于检测空化效应的方法及装置、超声治疗设备

Also Published As

Publication number Publication date
DE102014001994B4 (de) 2019-01-17
JP2014161735A (ja) 2014-09-08
US20140243667A1 (en) 2014-08-28
KR102223052B1 (ko) 2021-03-03
KR20140105410A (ko) 2014-09-01
JP6272707B2 (ja) 2018-01-31
US9168027B2 (en) 2015-10-27
CN104013419B (zh) 2019-06-28
DE102014001994A1 (de) 2014-08-28

Similar Documents

Publication Publication Date Title
CN104013419A (zh) 医疗超声波中的自适应声压估计
Oelze Bandwidth and resolution enhancement through pulse compression
US6827686B2 (en) System and method for improved harmonic imaging
KR102170630B1 (ko) 초음파의 음향 흡수 또는 감쇠의 측정
KR101868381B1 (ko) 의료용 초음파 이미징에서의 전단파 정보의 해석
JP5490198B2 (ja) 幅広ビーム映像化
KR102332464B1 (ko) 분석 데이터로부터 전단파 추정
KR102051293B1 (ko) 의료 초음파 전단파 이미징에서의 분류 프리프로세싱
KR20150037689A (ko) 의료용 초음파 영상의 전단파 검출
KR102628879B1 (ko) 기본 및 고조파 신호들을 사용하는 정량적 초음파
US7056290B2 (en) Continuous depth harmonic imaging using transmitted and nonlinearly generated second harmonics
KR20170094521A (ko) 전단파들을 사용하는 사운드 스피드 이미징
US20230404537A1 (en) Ultrasound medical imaging with optimized speed of sound based on fat fraction
KR20220036863A (ko) 매질의 초음파 특성화 방법 및 시스템
JP2003180688A (ja) 幅広ビーム映像化
Kanzler et al. Improved scatterer size estimation using backscatter coefficient measurements with coded excitation and pulse compression
KR20230109615A (ko) 매질의 초음파 특성화 방법 및 시스템
Pinton et al. Characteristics of the spatial coherence function from backscattered ultrasound with phase aberration and reverberation clutter
Clement et al. The feasibility of pulse compression by nonlinear effective bandwidth extension
JP5489758B2 (ja) 幅広ビーム映像化
CN107205662A (zh) 医学超声成像中的基于谐波与基波关系的发射功率
Merks Instantaneous ultrasonic assessment of urinary bladder volume
Pereira et al. Ultrasonic velocity mapping of multilayered media
Umeki Applying synthetic aperture, coded excitation, and tissue harmonic imaging techniques to allow ultrasound imaging with a virtual source
Toulemonde et al. Speckle reduction with multitaper approach to improve B/A imaging

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant