CN110622034B - 超声诊断图像中的混响伪影消除 - Google Patents

Abstract

一种超声诊断成像系统采集响应于多个发射事件而产生的回波信号的接收到的波束。接收到的波束与重新聚焦进行组合以考虑接收波束到发射事件位置中的差异。在重新聚焦中使用的延迟和权重被补充有校正混响伪影的延迟和权重。接收到的回波信号被处理以检测混响伪影的存在，并且混响信号分量到图像场中的虚拟点源的模拟发射被计算。该模拟产生用于混响信号补偿的延迟和权重、或能够从接收到的回波信号减去以减少混响伪影的估计的混响信号。

Description

超声诊断图像中的混响伪影消除

相关申请的交叉引用

本申请要求2017年5月11日递交的美国临时US 62/504681的权益和优先权，通过引用将其整体并入本文。

技术领域

本发明涉及医学诊断超声系统，并且具体地涉及减少或消除由于超声图像中的混响回波的图像伪影的超声系统。

背景技术

医学超声图像能够被源于各种源的图像伪影污染。除了能够在令人不满意的信噪比尚存时出现的简单噪声伪影之外，其他伪影对于超声的现象是特定的。斑点伪影由于超声信号的相干性质而出现。在图像中能够作为轻微水印出现的这些伪影能够通过信号处理技术(诸如频率复合和空间复合)来减少。像差伪影能够由于去往和来自换能器的声音通过超声路径中的不同组织和物质速度的差异而出现，并且能够通过波束形成过程中的延迟补偿来减少。特定于超声的另一图像伪影是混响伪影。当发射的超声波被图像区域的近场中的强反射器反射回来并且行进回到充当再次向外反弹返回波的反射器的换能器的面时，混响发生，由此在回波接收期间将第二向外波引入到图像场。该混响波将导致其自己的回波返回，其将与从发射波返回的回波混在一起。尽管混响回波处于比从发射波返回的那些更低的水平，但是它们常常具有足够的幅度来产生被叠加在主要期望图像上的它们自己的部分体模图像。

美国专利US 6905465(Angelsen等人)描述了一种用于通过发射两次(对混响信号伪影的返回信号进行采样的第一次，然后发射被调节以减少混响的效应的第二次)来校正超声成像中的混响像差的技术。然而，该方法需要两个发射事件，这增加了采集图像数据所需的时间并且因此降低了显示的帧率。

发明内容

本发明旨在期望地从超声图像减少或消除混响伪影，而无需降低显示帧率的多次发射。

根据本发明的原理，描述了减少超声图像中的混响伪影的出现的诊断超声系统和方法。方法和系统通过首先通过操作在一组回波信号数据上来检测超声图像数据中的混响伪影的存在和位置来操作。产生混响伪影的接收到的回波信号数据的信号分量被估计，优选地使用追溯性动态发射聚焦的原理。估计的混响信号从实际接收到的信号减去，或偏移相位和权重调节在波束形成过程中使用以从用于产生超声图像的图像数据减少或消除混响伪影。

附图说明

在附图中：

图1a)至e)图示了超声图像中的混响伪影的现象。

图2以方框图形式图示了根据本发明的原理构建的超声系统。

图3图示了追溯性动态发射聚焦原理。

图4a)至c)图示了在混响伪影信号的估计中被分解到虚拟点源的复波前的模拟。

图5图示了由图2的混响信号处理器执行的混响信号检测和处理。

具体实施方式

图1a)至e)的附图图示了超声成像中的混响伪影的问题。在图1a)中，示出了朝向目标解剖结构10发射超声能量20的平面波的换能器阵列104。在换能器与目标之间的路径中的是多个镜面反射器12，其会将发射的超声能量中的一些朝向换能器反射回来。这些反射器可以是不同的组织、骨骼(诸如肋骨或颅骨)、气腔或存在显著声阻抗不连续的其他结构。图1b)示出了正从反射器12朝向换能器被反射回来的回波22。在图1c)中，回波22已经被反射离开换能器阵列104的面，并且反映之前被反射器12返回的二次(混响)能量波22’现在向后传播通过图像场。同时，原始发射波已经从目标解剖结构10被反射回来，并且来自那里的回波正返回到换能器阵列，如通过波前24示出的。当混响能量22’到达反射器12时，能量中的一些第二次朝向换能器被反射回。在图1d)中的26处示出了返回到换能器阵列的这些次级混响回波，其中，它们与从目标解剖结构返回的回波混合。

图1e中的屏幕124上示出了由这些不同回波产生的图像。图像的左侧上的镜面反射器12的清晰图像由被反射器返回的第一组回波22的接收构建。目标区域10的图像由被目标解剖结构返回的回波24的接收构建。但是在该范例中，目标解剖结构由被反射器12第二次返回的混响回波的伪影部分地遮掩，由于已经使它们与来自目标解剖结构的返回回波混合的其行进时间的原因，其产生叠加目标解剖结构的图像的重建的体模图像。本发明的目的是检测混响回波并且消除其在图像重建中的效应，但是在不使用通过阵列换能器进行的任何另外的“询问”发射的情况下这样做，这会具有降低图像的显示的帧率的不期望效应。

图2以方框图形式图示了根据本发明的原理构建的超声成像系统。超声探头102包括换能器元件的换能器阵列104。选定的换能器元件组分别在相位延迟时刻由发射波束形成器106致动，以在期望方向上并且从沿着阵列的(一个或多个)期望原点发射在选定焦点区域处操纵并聚焦的波束。发射波束形成器通过发射/接收开关被耦合到换能器元件，所述发射/接收开关可以包括保护接收器输入端免受所施加的高电压发射脉冲伤害的交叉点开关。响应于每个发射波束而被阵列104的每个换能器元件接收的回波被应用于多线处理器110a-110n的输入端。每个多线处理器包括接收波束形成器，所述接收波束形成器应用其自己的一组延迟和(如果需要)对从阵列元件接收的回波进行加权的切趾权重，以根据相同发射波束形成不同操纵和聚焦的接收波束。例如，在美国专利US 6695783(Henderson等人)和美国专利US 5318033(Savord)中描述了用于多线处理器110a-110n的合适的多线波束形成器。多线处理器110a-110n的扫描线输出被耦合到存储接收到的多线的线存储设备112，至少直到形成显示数据的扫描线的所有多线已经被采集。

接收的多线被组合器90组合，组合器90在其扫描转换之前执行处理接收到的信号。组合器可以包括若干个单元，诸如乘法器116、加权电路114、延迟118和加法器120。用于形成显示数据的具体线的该组多线被应用于乘法器116a-116n中的相应乘法器，以针对对应扫描线位置产生显示数据。如果期望，来自每条线的回波数据可以通过切趾权重114a-114n来进行加权。一般来说，这些权重将根据其往返(round-trip)脉冲响应对每条线进行加权。合适的加权算法能够通过让术语幅度(x，y)为通过发射波前的图像场中的位置(x，y)处的点的声穿透幅度来导出，方位角位置x＝0对应于发射波束的中心轴。让X为接收到的多线相对于发射波束轴的方位角。被应用于该接收到的多线以在深度Y处形成图像的点的权重为：

权重(X，Y)＝幅度(X，Y)

为了适当延迟特性的确定，让传播_时间(x，y)为发射波前到达位置(x，y)处的点所需的传播时间，方位角x＝0再次对应于发射波束的中心轴。让X为接收到的线相对于发射波束轴的方位角。被应用于该接收到的多线以在深度Y处形成图像的点的延迟为：

延迟(X，Y)＝传播_时间(X，Y)-传播_时间(0，Y)

其中，传播_时间(0，Y)为到达在相同深度处但是在轴上的点的时间。

函数幅度(X，Y)和传播_时间(X，Y)可以例如从发射场的模拟获得。计算传播时间的适当方式是使用来自在若干频率下的单色模拟的场的相位延迟。幅度可以通过平均在若干频率下的场的幅度来计算。此外，依赖深度的归一化能够被应用于权重。这将给定深度处的所有权重乘以公共因子。例如，归一化能够被选择为使得斑点区域具有随深度的一致亮度。通过根据深度改变权重，能够随着深度动态地改变孔的尺寸和形状(切趾)。

幅度和传播时间不需要根据在系统中使用的精确发射特性的模拟导出。设计者例如可以选择使用不同的孔尺寸或不同的切趾。

来自每条线的回波由乘法器116a-116n进行加权，并由延迟线118a-118n进行延迟。通常，如上所示，这些延迟与发射波束中心到接收线位置的位置相关。延迟被用于平衡具有不同发射-接收波束位置组合的多线间从线到线存在的相移变化，使得不会由从不同发射孔组合的信号的相位差引起信号消除。

将意识到在数字系统中，能够通过在存储器中存储经加权的多线回波数据并且稍后读出数据(这实现了必要的延迟)来实现延迟线。也可以使用不同长度的移位寄存器和时钟信号来实现数字延迟，或者可以使用如前述美国专利US 6695783中描述的内插波束形成器。延迟的信号由加法器120进行组合，并且得到的信号被耦合到图像处理器122。图像处理器可以执行扫描转换或其他处理以改进所显示的图像。在图像显示器124上显示得到的图像。

在图5的系统中，延迟线118和加法器120实现了对从若干接收多线接收的信号进行重新聚焦，所述接收多线在给定方向上共对准。重新聚焦操作为从使用针对每条多线的不同发射波束位置引起的相位差进行调节，从而防止组合信号中的非期望相位消除。权重114相对于发射波束与多线位置的接近度对多线的贡献进行加权，从而将更高的权重给予具有更高信噪比的接收波束。这产生沿每条接收线的扩展的场深以及由于在每个接收线方向上多个采样的组合的增强的穿透性(改进的信噪比)。

共对准的接收到的多线的这种重新聚焦还引起如参考图3解释的追溯性动态发射聚焦效应。在该附图中，四个发射波束34由换能器阵列104发射，发射事件被称为Tx1、Tx2、Tx3和Tx4。每个发射波束在图像场中在不同的方向上被发射；在该范例中，发射和接收波束平行，并且因此每个发射波束源于沿着换能器阵列的不同位置(x维度)。每个发射的能量具有通过沙漏形状的线30大致指示的发射波束轮廓，沙漏形状的线30在通过箭头32指示的焦点处收敛。因此，每个发射波束通过阵列的其自己的发射子孔来生成。在每次发射后，与发射线共定位的接收线被接收并且被波束形成，并且额外的接收波束40由阵列换能器接收，每次在沿着阵列的相同位置(x方向)处。两个回波信号位置通过每个接收波束40上的相同深度y处的圆圈来标记，一个在浅深度处而另一个在更深的深度处。在42处示出了针对第一发射事件Tx1的接收波束。应看出对于第一发射事件，接收波束40被偏移到发射波束34的右侧，并且由于距离和时间在波束形成中是相当的，从画圈回波信号位置接收的回波如图所示被定位于在波束形成后的接收波束42上。

在Tx2的第二发射-接收循环中，发射波束34横向地更靠近接收波束40。该更小的偏移得到两个画圈的回波信号位置被定位为更靠近接收波束的焦点，如在第二接收波束42’上示出的。在Tx3的第三发射-接收循环中，发射波束被横向地定位到接收波束40的右侧，被分开与Tx2循环相同的横向距离。画圈的回波信号位置的回波被定位，如通过接收波束42”上的信号示出的。在Tx4的第四发射-接收循环中，发射波束还被定位到接收波束中心的右侧，以与在第一Tx1循环的情况下相同的距离偏移。画圈的回波信号位置的回波被定位，如通过接收波束42’”上的信号示出的，与接收波束焦点分开更大的距离。画圈的回波信号位置的变化可以通过跨过四个接收波束42、42’、42”和42’”绘制的曲线44和46来表示。该变化在追溯性动态发射重新聚焦中通过在对它们进行组合之前将延迟校正应用于相应的接收多线来校正。例如在美国专利US 8137272(Cooley等人)中示出了示范性延迟曲线轮廓。当这种补偿延迟借助于图2中的延迟线118a-118n被应用时，画圈的回波位置全部被重现在沿着多线42-42’”的相同的相应深度处，如通过画圈的深度45和47示出的，画圈的深度45和47在如通过笔直曲线54和56示出的对应接收线深度处。四个校正的多线然后可以与对应于发射到接收线偏移的加权进行组合，以产生具有组合的回波信号52和58的最终显示扫描线50。

根据本发明的原理，图2的超声系统还包括混响信号处理器100。混响信号处理器的目的是识别从图像场接收的回波信号中的混响信号伪影的存在，然后通过接收到的混响信号分量回到位于发射波束焦点的虚拟点源的时间反演的过程，估计将用于生成混响信号的定相和加权，如果实际上它们要被发射和接收的话。但是没有额外的发射和接收被采用。替代地，估计的相位和权重被用于补充在追溯性动态发射重新聚焦中使用的相位和权重，由此针对用于成像的接收到的回波信号中的混响信号伪影获得消除信号。

该过程以检测接收到的回波信号中的混响信号伪影的存在开始。这通过操作在接收到的多线的包络检测的信号上来完成。在图2的实施方式中，通过多线处理器产生的初始波束形成信号通过包络检测器92进行包络检测。作为混响的源的强反射器能够被定义为在近场中并且靠近换能器的任何图像像素p(x，y)，诸如在图像的前一半内的像素，即，y≤y_max/2，其中，y指示图像深度。从该强反射器接收的信号也应当具有大于强度范围的阈值的值，例如，最大像素强度的至少80％。例如，如果像素值范围为0至255，80％阈值将会为204，确保怀疑包含引起混响的信号的任何像素具有最够大的幅度。怀疑的像素的强度值应当为比由相邻像素组成的感兴趣区域的平均强度大至少两个标准差，以确保总体增加的图像增益不引起高强度。一旦这些亮的高强度像素在图像信号中被识别，然后对应于这些高强度点的波束形成的r.f.线(A线)被研究以找到强回波。在该范例中，波束形成的A线是检测到的多线信号，并且像素位置p(x，y)对应于A_i(t)处的信号，其中，A_i是第i条A线。对于线性阵列几何结构，图像深度y与当回波发生时的时间(t_e)直接相关，y＝c₀t_e/2，其中，c₀是声音在成像介质中的速度，并且线数i是到p的横向x-位置的最靠近的超声波束。这可以被表达为i＝((x-x₁)/Δx)+1的积分值，其中，x₁是第一A线的横向位置，并且Δx是两个相继的A线之间的距离。对于扇形扫描几何结构，其中，波束被角度地散开并且空间地起源于公共顶点位置处，扫描转换需要首先被反转，使得相对于通过(x₀，y₀)给出的扫描原点的深度r与回波的信号时间相关，并且i是i＝[(tan((x-x₀)/(y-y₀))-θ₁)/Δθ]+1的积分值，其中，θ₁是第一A线的角度，并且Δθ是两条相继的A线之间的角度。A_j(t)附近的rf-信号的包络被计算以找到A_j(t)的峰周围的信号，直至包络下降至峰值幅度的一半。具有这些特性的回波被识别为次级发射的回波，即，混响伪影。

让S_ij(t)指示在第j个发射(Tx_j)后在第i个接收器元件处接收到的信号。A_j通过针对阵列的所有接收元件对信号s_ij进行波束形成来获得，如通过图2中的多线处理器110a-110n指示的。被用于构建A_j的波束形成延迟被用于识别接收到的信号，使得如果来自接收器i的信号被延迟φ，并且回波时间被计算为t_e，那么t_e+φ周围的信号对混响负责。为此目的，来自探头和混响信号处理器的未处理的回波数据之间存在直接通信，如在图2中示出的。波束形成延迟能够从多线处理器实时地读取，或备选地正在多线处理器中使用的相同延迟也能够被预先加载到混响信号处理器。根据检测到的信号包络，峰t_e+φ的两侧上的信号被检查，直到包络幅度下降至最大幅度的一半。该回波位置和信号水平被保存用于混响移除处理。这些信号在下文中被称为s^rev(t)。来自所有接收元件的信号s^rev _ij(t)定义在第j个发射后从换能器表面反射的混响波前。该波前通常是复数的(非平面的、非收敛的、非均匀切趾的)，因为它们从其反射的解剖表面很少是平坦的。

接下来，混响信号处理器使用复波前作为来自信号路径的换能器侧的输入来执行混响波传播的模拟。模拟的波22’从换能器阵列104朝向如在图4a)和b)中针对扇形几何结构示出的波束的焦点130被向外传播，以将波前分解成多个虚拟点源。传播的波前22”的幅度和到达时间针对如在图4b)中示出的波束的焦点130来计算。这通过使用第i个元件和第j个聚焦发射的信号s^rev _ij(t)构建发射-接收变换矩阵K_TxRx(t)来完成：

其中，M是换能器元件的数量，并且N是聚焦发射的数量。K_TxRx矩阵的每个列表示在聚焦发射后每个元件的接收到的混响信号数据。类似地，聚焦发射矩阵K_焦点(t)被构建为：

其中，s_ij(t)指代从用于N次聚焦发射的M个换能器元件中的每一个发射的信号。该计算通过发射波束形成器106和混响信号处理器之间的数据线来促进，如在图2中示出的。该矩阵包括与发射波束形成相关联的延迟和切趾权重，其能够被预加载到混响信号处理器。在时间维度上的傅里叶变换之后这两个变换矩阵也能够在频域中进行表示，K_TxRx(ω)和K_焦点(ω)。

能够使矩阵K_TxRx的任何列与矩阵K_焦点相乘，以模拟混响信号到焦点(虚拟源130)的传播。因此，

V(ω)＝K_焦点(ω)K_TxRx(ω)

其中，V(ω)是其列表示反射的混响回波到N个虚拟源130的分解的N乘N矩阵。结果是N个元素的复(相位和幅度)向量。该计算能够针对每个频率ω进行重复，通常通过对应于信号采样频率的奈奎斯特范围的n个不同频率的傅里叶变换，以将来自每个发射事件的混响波前分解成其虚拟源分量。V(ω)的逆傅里叶变换被执行以回朔到焦点处的到达时间和波幅度，其将会充当用于通过上面描述的追溯性动态发射聚焦调节执行的校正中的虚拟源的延迟时间和权重。图4c)图示了根据延迟时间t_延迟被重新定位的虚拟源点130，所述延迟时间t_延迟表示在虚拟源点处用于校正混响伪影所需的延迟时间。

已经被存储在线存储设备112中的来自个体发射事件Tx_j的接收到的A线在应用标准追溯性动态发射聚焦校正之后被组合(求和)，以估计用于以下方程中的次级(混响)发射A^rev的A线：

A^rev(ω)＝[V(ω)]^HA(ω)

其中，上标H指示矩阵的厄米算符。在该表达式中，A(ω)指示在傅里叶变换之后的A线矩阵表达式的列向量，并且A^rev指示在次级发射实际已经被发射的情况下到接收到的且波束形成的A线的近似。然而，实际上没有物理波束发射被完成，并且对混响伪影负责的信号通过A^rew(ω)的逆傅里叶变换被近似地计算。备选地，来自个体发射的每信道数据(代替于波束形成的A线)用于追溯性动态发射重新聚焦，并且被组合(求和)为包含所计算的权重和延迟以针对次级(混响)发射估计接收到的信号。

最后通过减法从实际接收到的信号消除虚拟次级发射的估计的混响回波信号：

A^校正(ω)＝A(ω)-A^rev(ω)

经校正的A线通过逆傅里叶变换A^校正(ω)来获得。尽管混响消除过程已经在频域中进行图示，但是其也能够在时域中被执行。在时间延迟和权重的计算之后，延迟能够使用位移位技术被应用于波束形成的A线，并且权重能够与延迟的A线相乘，并且求和能够在时域中被执行以获得A^rev(t)。最后，从A(t)减去A^rev(t)以获得A^校正(ω)。

在图5的流程图中顺序地图示了前述的处理和混响伪影校正。方框70图示了使用一系列N个发射事件(Tx#1、Tx#2、…、Tx#N)的B模式发射-接收序列产生超声图像的方法，来自所述一系列N个发射事件的回波数据被用于产生超声图像。每个发射事件导致回波信号由换能器阵列中的每个元件(信道)的接收，该每信道(每元件)数据被示为N个发射事件的Rcv#1…Rcv#N。来自每个发射事件的回波数据被波束形成为A线，因此形成A-L#1…A-L#N。在步骤72中，来自Rcv#1的回波数据集被分析以识别对混响负责的强回波，其能够包括rf信号的包络检测，并且搜索大于阈值并且比那些周围像素大多个(例如，两个)标准差的幅度峰值。在步骤74中，所识别的混响信号波前从换能器朝向虚拟点源130的传播的模拟使用上面给出的信号矩阵方程(例如，V(ω))来执行。V(ω)的逆傅里叶变换产生焦点处的到达时间和信号幅度，其充当用于混响校正的延迟时间和权重，如在步骤76中指示的。在步骤78中，权重w_i和延迟t_i被应用于来自N个发射的波束形成的A线A-L#1…A-L#N，并且在步骤80中，加权的且延迟的A线被求和以估计接收到的混响信号(步骤82)。该方法还可以包括从A线消除混响信号的步骤84。该过程针对Rcv#2回波数据集进行重复(步骤88)以从该数据集消除混响信号，并且在所有数据集都已经被处理之后，图像可以使用经校正的A线来形成(步骤86)。

应当注意，适合于在本发明的实施方式中使用的超声系统和具体地图2的超声系统的部件结构可以以硬件、软件或其组合来实施。超声系统的各种实施例和/或部件(例如，多线处理器、加权和延迟电路、包括检测器和混响信号处理器、和图像处理器、以及其中的部件和控制器)也可以被实施为一个或多个计算机或微处理器的部分。计算机或处理器可以包括计算设备、输入设备、显示单元和例如用于访问因特网的接口。计算机或处理器可以包括微处理器。微处理器可以连接到通信总线，例如，以访问PACS系统或用于输入训练图像的数据网络。计算机或处理器还可以包括存储器。存储器设备(诸如线存储设备112)可以包括随机存取存储器(RAM)和只读存储器(ROM)。计算机或处理器还可以包括存储设备，其可以是硬盘驱动器或可移除存储驱动器，诸如软盘驱动器、光盘驱动器、固态拇指驱动器等。存储设备也可以是用于将计算机程序或其他指令加载到计算机或处理器中的其他类似模块。

如本文所使用的，术语“计算机”或“模块”或“处理器”或“工作站”可以包括任何基于处理器或基于微处理器的系统，包括使用微控制器、简化指令集计算机(RISC)、ASIC、逻辑电路和能够执行本文所描述的功能的任何其它电路或处理器的系统。上述范例仅是示范性的，并且因此不旨在以任何方式限制这些术语的定义和/或含义。

计算机或处理器执行存储在一个或多个存储元件中的指令集，以便处理输入数据。存储元件还可以根据需要或期望来存储数据或其他信息。存储元件可以采取处理机器内的信息源或物理存储器元件的形式。

超声系统的指令集(包括控制如上面描述的超声图像的采集和处理的那些)可以包括各种命令，其指导计算机或处理器作为处理机器以执行诸如本发明的各种实施例的方法和过程的特定操作。指令集可以采取软件程序的形式。软件可以采取诸如系统软件或应用软件的各种形式，并且其可以被实现为有形和非瞬态计算机可读介质。此外，软件可以采取单独的程序或模块的集合的形式，诸如执行上面描述的混响信号处理器的方程的模拟和处理的程序或模块。软件还可以包括面向对象编程形式的模块化编程。由处理机器对输入数据的处理可以响应于操作者命令，或者响应于先前处理的结果，或者响应于由另一处理机器进行的请求。

此外，对以下权利要求的限制不以模块加功能的格式写出，并且不旨在基于35U.S.C.112第六段来进行解释，除非和直到这样的权利要求限制明确地使用短语“用于…的模块”，之后是没有进一步结构的对功能的说明。

Claims (15)

1.一种用于减少超声图像中的混响信号伪影的方法，包括：

接收包含混响信号伪影的超声回波信号；

通过包络检测来处理所述超声回波信号；

根据包络检测的超声信号的幅度来在所述超声回波信号中检测混响信号分量的存在；

根据对混响信号到虚拟点源的发射的模拟来估计用于混响信号校正的时间延迟和幅度；

使用所估计的时间延迟和幅度来组合所述超声回波信号的多条A线以减少所组合的A线中的混响信号分量；并且

使用所组合的A线来形成超声图像。

2.根据权利要求1所述的方法，其中，接收超声回波信号还包括发射多个发射事件，并且响应于每个发射波束而接收多个多线。

3.根据权利要求2所述的方法，其中，组合多条A线还包括组合来自多个发射事件的多个多线。

4.根据权利要求3所述的方法，其中，组合多个多线还包括使用相对于彼此重新聚焦多个多线的时间延迟来组合所述多线。

5.根据权利要求1-4中的任一项所述的方法，其中，检测混响信号的存在还包括识别大于特定阈值并且比周围像素大至少两个标准差的超声信号像素幅度。

6.根据权利要求1-5中的任一项所述的方法，其中，所述虚拟点源被定位于发射的或接收到的超声波束的焦点深度处。

7.根据权利要求1-6中的任一项所述的方法，其中，估计用于混响信号校正的时间延迟和幅度还包括根据接收到的回波信号来估计混响信号分量。

8.根据权利要求7所述的方法，其中，所述虚拟点源被定位于发射的或接收到的超声波束的焦点深度处；并且

其中，所述时间延迟和幅度是根据所述虚拟点源处的模拟的混响信号分量的到达时间和幅度来估计的。

9.根据权利要求7或8所述的方法，其中，使用所估计的时间延迟和幅度来组合所述超声回波信号的多条A线以减少混响信号分量还包括从接收到的超声回波信号减去所估计的混响信号分量。

10.根据权利要求8或9所述的方法，其中，组合多条A线还包括使用重新聚焦延迟来组合多条A线以估计包含混响伪影的A线。

11.一种用于减少超声图像中的混响信号伪影的诊断超声系统，包括：

多线处理器(110a、...、110n)，其被布置为通过对包含混响信号伪影的多个接收到的超声回波信号进行至少部分波束形成来提供多条A线；

包络检测器(92)，其响应于所述多条A线，并且被布置为检测每条A线的包络，其中，每个包络检测包括包络幅度检测；以及

混响信号处理器(100)，其被耦合到所述包络检测器，并且被布置为根据所述包络的幅度来在所述超声回波信号中检测混响信号分量的存在，并且根据对存在的混响信号到虚拟点源的发射的模拟来估计用于混响信号校正的时间延迟和幅度。

12.根据权利要求11所述的诊断超声系统，还包括：

组合器(90)，其被耦合到所述多线处理器和所述混响信号处理器两者，并且被布置为使用所估计的时间延迟和幅度来组合所述超声回波信号的所述多条A线以减少所组合的A线中的混响信号分量；以及

图像处理器(122)，其被布置为使用所组合的A线来形成超声图像。

13.根据权利要求11或12所述的诊断超声系统，其中，所述混响信号处理器被布置为通过还计算聚焦的发射信号矩阵与发射-接收变换矩阵的乘积来估计用于混响信号校正的时间延迟和幅度。

14.根据权利要求13所述的诊断超声系统，其中，所述混响信号处理器被布置为通过还计算复混响波前矩阵的逆傅里叶变换来估计用于混响信号校正的时间延迟和幅度。

15.根据权利要求12、13或14中的任一项所述的诊断超声系统，其中，所述图像处理器被布置为通过根据多个多线形成图像来形成超声图像，所述多个多线已经根据其发射波束的轴被重新聚焦。