用于利用持续换能器振动进行超声弹性成像的系统和方法

摘要

提供了用于处理使用超声弹性成像获取的数据的系统和方法，在所述系统和方法中，使用超声换能器的持续振动在对象中生成剪切波。在此描述的系统和方法可以有效地去除与所述超声换能器的振动相关联的运动伪影，并且还可以去除逐行成像模式用于获取数据(如许多常规超声扫描器完成的那样)时所造成的数据采样未对准。因此，在此描述的系统和方法提供用于换能器运动校正以及用于对准由逐行扫描超声系统检测到的运动信号的技术。

说明书

本申请是申请日为2016年10月06日、申请号为201680058937.1、题为“用于利用持续换能器振动进行超声弹性成像的系统和方法”的分案申请。

相关-申请的交叉引用

本申请要求于2015年10月8日提交的标题为“SYSTEMS AND METHODS FORULTRASOUND ELASTOGRAPHY WITH CONTINUOUS TRANSDUCER VIBRATION(用于利用持续换能器振动进行超声弹性成像的系统方法)”的美国临时专利申请序列号62/238,891的权益。

关于联邦政府赞助的研究的声明

本发明是在美国国家卫生研究院授予的DK106957下利用政府支持完成的。美国政府享有本发明的某些权利。

技术领域

本发明的领域是用于超声弹性成像的系统和方法。更具体地，本发明涉及用于处理使用超声弹性成像获取的数据的系统和方法。

背景技术

超声剪切波弹性成像(“SWE”)已经作为新的超声成像技术出现，所述技术可以非侵入性且量化地估定组织机械特性，所述组织机械特性是组织健康状态的较强生物标记。通常，在SWE中，将剪切波引入组织中，并且使用脉冲-回波超声检测剪切波的传播参数。然后，剪切波的检测用于计算与组织机械特性相关的参数，包括剪切波传播速度、扩散(即，频率相依性)、剪切波衰减、剪切模量、剪切粘度、杨氏模量、存储模量、损耗模量、损耗正切角和机械弛豫时间。

常规超声SWE使用声学辐射力(“ARF”)来在组织中生成剪切波。ARF需要从超声换能器传输长持续时间的推动脉冲，这要求在下一次传输之前很长的冷却时间，从而避免可能的探测和组织加热。这从根本上限制了超声SWE的帧速率(例如，限制为大约1Hz)。ARF还具有对超声系统的高电源要求，这使得在中端和低端超声扫描器中实现有挑战性。

为了解决这些限制，利用超声换能器的持续振动的超声弹性成像技术(例如，在共同待决美国专利申请序列号62/072,167中所描述的技术)。这种技术通过换能器的持续振动来生成剪切波，并且利用相同的换能器来检测已生成的剪切波信号。因为这种技术不使用ARF来进行剪切波生成，所以这种技术允许持续的高帧速率剪切波成像以及利用中端和低端超声系统的便利实现方式。

然而，换能器的持续振动还引入了剪切波信号处理的挑战。一个主要挑战是针对换能器的运动校正所获取数据，并且另一个主要挑战是利用逐行扫描超声系统成像时的运动信号对准。

发明内容

本发明通过提供一种用于使用具有换能器的超声系统来测量对象的机械特性的方法克服了前述缺点。持续振动被提供至所述超声换能器，由此所述超声换能器的振动在所述对象中引入至少一个剪切波。然后，使用超声换能器从所述对象处获取所述运动数据。所述运动数据指示在所述对象内传播的所述至少一个剪切波。然后指示由所述超声换能器的持续振动造成的所述对象的变形的压缩简档被估计并且被用于通过将所述压缩简档从所述所获取运动数据中解调、分离或以其他方式去除来产生校正数据。然后对所述校正数据进行处理以计算所述对象的机械特性。

本发明的另一方面是提供一种用于使用具有换能器的超声系统来测量对象的机械特性的方法。持续且周期性振动被提供用于所述超声换能器，由此所述超声换能器的振动在所述对象中引入至少一个剪切波。然后，使用超声换能器从所述对象处获取所述运动数据。所述运动数据指示在所述对象内传播的所述至少一个剪切波，并且所述运动数据是在被选择用于减轻运动误差的时间点获取的，所述误差可归因于由所述持续且周期性振动造成的所述对象的变形。然后对所述运动数据进行处理以计算所述对象的机械特性。

本发明的另一方面是提供一种用于使用具有换能器的超声系统来测量对象的机械特性的方法。剪切波包括在对象中并且使用采用脉冲-回波模式的所述超声换能器从所述对象处获取运动数据。所述运动数据指示在所述对象内传播的所述剪切波。然后，由于所述对象中不同位置处的数据获取之间的时间延迟造成的误差而对所述运动数据进行校正，并且处理校正数据以计算所述对象的机械特性。

本发明的前述及其他方面和优点从以下说明书中将变得明显。在说明书中，参照在此构成其一部分的附图，并且在附图中通过图示的方式示出了本发明的优选实施例。然而，这样的实施例并不一定表示本发明的全部范围，并且因此参考权利要求书并在此用于解释本发明的范围。

附图说明

图1是示例超声系统的框图，所述超声系统实现了超声换能器的持续振动以将剪切波引入对象中。

图2A至2C展示了使用超声换能器的持续振动在对象中引入剪切波和附加变形的示例。

图3是流程图，阐述了用于响应于连续振动的超声换能器而获取运动数据并且针对换能器振动的效应来校正所获取数据的示例方法的步骤。

图4A至4B展示了估计来自从所获取运动数据生成的k空间数据的压缩简档的示例。

图5展示了实现在表示连续换能器振动的正弦信号的波峰(圆形)和波谷(方形)附近的脉冲-回波运动检测的示例。

图6展示了运动信号对准方法的示意曲线。运动检测的时间对准展示在区域1内。针对区域2描绘所提出的区域间对准方法。实心方形指示实际脉冲回波事件。空心圆形指示适时插入的数据点。具有虚线轮廓的灰色方形指示基于区域间对准的时移运动信号数据点。

图7展示了使用相移技术的剪切波信号对准的示例。

图8展示了用于低PRF

图9描绘了传播剪切波和非传播压缩简档的k空间表示。

具体实施方式

在此描述的是用于处理使用超声弹性成像获取的数据的系统和方法，在所述系统和方法中，使用超声换能器的持续振动在对象中生成剪切波。在此描述的系统和方法可以有效地去除与所述超声换能器的振动相关联的运动伪影，并且还可以去除逐行成像模式用于获取数据(如许多常规超声扫描器完成的那样)时所造成的数据采样未对准。

因此，在此描述的系统和方法提供用于换能器运动校正以及用于对准由逐行扫描超声系统检测的运动信号的技术。

首先参照图1，展示了用于将剪切波引入对象20的示例系统10。在这一系统中，换能器12通过致动器14而机械振动，这使得换能器12在轴向方向(例如，图1中的z方向)上振荡。作为一个示例，致动器14可以是机械致动器，如，音圈致动器。当换能器12在轴向方向上移动时，剪切波16被引入到对象18之内。还可以从压缩波的模式变换中产生剪切波16。然后由在脉冲-回波模式中操作的相同超声换能器12检测所产生的剪切波16以提供对所述对象的机械特性的定量测量。

致动器14耦合到超声换能器12。作为一个示例，致动器14可以直接附接到换能器12的外表面。出于说明目的，致动器14附接到图1中的换能器12的一侧。然而，在一些申请中，可以优选地将致动器与换能器12同轴地对齐，从而使得换能器运动主要与最小横向和纵向运动同轴。这一设置消除了对于单独振动源的需要，并且因此允许进行方便的单手操作。所述振动优选地是持续的以允许测量的持续更新。

超声换能器12可以沿超声波轴向或者取决于所期望的成像申请在其他方向上轴向地进行振动。可以操作用于剪切波检测的超声系统以通过平行波形成来检测单个A线、多个A线，或者利用平面波成像和软件波形成来检测整个2D面积或3D体积，如在

通过致动器14施加至超声换能器12的持续振动可以包含多种频率，并且因此可以处理所检测的剪切波以解析所述对象的取决于频率的特性。例如，所述处理可以沿时间维度使用带通滤波器以便仅选择某一时间处的一种频率，并且如果收集具有单一振动频率的数据，则后续处理将与其相同。多频率振动可以加速采集以便进行分散分析。利用持续振动和持续剪切波检测和处理，可以采用基本上实时的方式来持续地更新弹性成像测量。

当换能器12在轴向方向上振动时，如当振动对于换能器12的有源表面20是正常的时，换能器12的运动将损害在对象18中所检测的剪切波信号。由于超声运动检测将换能器12用作非移动参考坐标，因此这一信号损害是存在的，但是当换能器12由于外部振动而震荡时，则违背了这个假设。因此，为了适当地测量来自所检测剪切波的机械特性，需要校正由致动器14造成的换能器12的运动。

在本公开的一方面中，提供了用于在连续换能器振动期间校正换能器运动的系统和方法。在超声中，通过比较两个脉冲回波事件之间的超声回波信号的时移τ来进行组织运动的测量。因为在软组织中超声传播速度c是恒定的(通常假设c＝1540m/s),所以这种时移可以被转换至组织位移为：

其中，因子二说明在脉冲回波超声检测中的往返距离。平均组织粒子速度v还可以通过以下进行计算：

其中，δ是两个脉冲回波事件之间的时间间隔。有时，多个脉冲回波事件的复合用于形成单一回波集合，以便提高信噪比(“SNR”)。在这种情况下，复合的脉冲回波事件各自可以由多个发送-接收过程组成。如在此使用的，术语“运动”可以包括位移、速度、加速度等。

如以下在图2A至2C中展示的，并且如上所描述的，超声中的运动检测使用超声换能器作为参考坐标。然而，当以持续换能器振动来振动换能器自身时，换能器在运动以便产生剪切波。这具有若干效应。首先，运动检测的参考坐标正在移动，这在给定时刻为所有成像的像素或体元添加了恒定运动偏移量。第二，换能器在物理上推动组织并且在组织中造成压缩变形。由于这种效应引起的组织运动随着每个成像的像素或体元甚至在相同时刻的位置而变化。第三，振动的换能器在组织中产生剪切波，所述剪切波是用于估计组织机械特性的信号。因此，重要的是去除由于前两种效应引起的运动，以便恢复由于用于准确计算组织机械特性的剪切波而引起的运动。

图2A示出了由超声换能器12在组织对象18的顶表面上施加的均匀压缩dz。由组织变形引起的位移取决于组织的位置和材料特性两者。举例来讲，图2B中的黑线50示出了在时间点t

在图2C处示出了在第二时间点t

现在参照图3，流程图展示为阐述了用于在持续换能器振动期间校正换能器运动的示例方法的步骤。所述方法总体上包括使用正在振动以在对象中生成剪切波的换能器从对象处获取数据，而数据是从所述对象中获取的，如在步骤302中所指示的。所获取数据指示在对象中总的所观察运动。在获取数据之后，估计表示数据获取期间换能器振动的压缩或变形简档，在步骤304中所指示的。根据多个不同的过程可以估计压缩简档，这在下面进行更详细地描述。使用所估计的压缩简档，将换能器运动的效应从所获取数据中解调、分离或以其他方式去除，如在步骤306中所指示的。作为一个示例，通过从所获取数据中减去压缩简档来从数据中去除换能器运动。根据校正数据，可以估计对象的机械特性，如步骤308处所指示的。

作为一种用于估计压缩简档的示例方法，可以实现曲线拟合过程。在此方法中，通过将总运动(如图2B中虚线所描绘的)拟合至已知函数可以对压缩简档进行估计。举例来讲，已知函数可以包括以下各项中的一项或多项：指数、多项式、幂律、样条、Flamant解或Boussinesq解。在一些实施例中，为了说明材料特性的局部变化，多个函数可以用于在不同深度对多个空间窗口的压缩简档进行估计。空间窗口可以是一维、二维或三维窗口。

为了改进对压缩简档的估计，在曲线拟合过程之前，所测量的总运动可以被去噪。举例来讲，可以使用滤波或正则化方法实现去噪。

在此描述的方法可以使用各种方法扩展到更高维度的图像。在一个示例中，之前描述的拟合方法可以扩展到多维对应物，从而允许针对一个或多个帧执行对压缩的轴向、横向和时间估计。术语“帧”可以指在给定时间处获得的二维(“2D”)超声回波数据集合，并且随着时间推移可以在相同的2D平面处获得多个帧。在另一个示例中，沿每个横向位置处重复的单一1D简档执行拟合，这允许从整个成像帧中去除换能器运动效应。在其他实现方式中，单一变形简档可以被应用于图像内的每个横向位置，以便进一步缩短计算时间。这种过程还可以在每个获取帧上重复(即，应用于不同时刻)，以便随着时间推移校正来自整个获取的换能器运动。

在上述曲线拟合示例中，从单一横向位置和单一帧沿着压缩方向估计换能器运动。像任何获取技术一样，所有的测量将包含一些误差。如此，利用来自多个空间位置和多个帧的信息将减少随机误差并且提供对真正压缩简档的更好估计。这可以通过至少两种方式来完成。

在一种方法中，使用平均数、加权平均数、中值或类似技术来组合在若干相邻横向位置处测量的多个总观察运动信号，以便获得用于曲线拟合和减法的噪声较少的总运动测量结果。

在另一种方法中，对来自每个帧的压缩简档进行估计，所估计的简档被组合成单一压缩简档，并且组合的简档从每个帧中减去以便去除压缩效应。尽管每个帧将在不同的时间点处获得，并且将组织压缩至不同的程度，但是作为第一阶近似值，可以假设的是压缩简档将是相互幅度按比例缩放的版本。

因此，可以对在每个横向位置处的不同帧的所有单独的压缩简档进行归一化，从而使得简档进将被组合(例如，使用取平均值)成单一压缩简档以用于随机噪声被抑制的每个横向位置。组合的压缩简档可以被缩放，并拟合至各个帧，并且然后被减去，以便获得在所述横向位置处的真正剪切波运动。对于所有横向位置可以重复相同的过程，以便获得跨2D区域的剪切波信号以用于进一步处理。需注意，这两种技术不互斥并且可以彼此结合使用。

作为另一种用于估计压缩简档的示例方法，可以获取并实现参考压缩简档。在此方法中，准静态压缩应用于对象并且通过使用相同的超声换能器的脉冲回波检测来对参考压缩简档进行估计。可以利用手动压缩或者通过以远低于在剪切波成像中典型使用的频率(例如，1Hz)来振动换能器从而实现准静态压缩。

可以假设，由于剪切波的运动在此情况中是可以忽略不计的；因此，所测量的运动简档应该只是由换能器压缩引起的。作为第一阶近似值，在不同压缩级别dz的压缩简档应该是彼此的缩放版本。因此，在单一压缩级别处获得的一个参考压缩简档应该是充足的。替代性地，多个压缩简档可以在不同压缩级别处获得，并且可以使用如上关于用于曲线拟合的组合压缩简档所描述的那些类似过程被缩放且被组合以形成具有更高信噪比(“SNR”)的单一参考压缩简档。

参考压缩简档可以被缩放、被拟合成所测得的总组织运动，并且从总组织运动中减去，以便获得真正的剪切波运动。以上关于曲线拟合所描述的空间和时间求平均技术还可以用在参考压缩方法中，以便提高SNR。

作为另一种用于估计压缩简档的示例方法，可以实现压缩简档建模。在此方法中，压缩简档可以从有限元方法(“FEM”)模拟或分析解决方案中获取。一旦已知已建模的压缩简档，其就可以被缩放、拟合并且从所测得的总组织运动中减去，如以上所述。作为第一近似值，对象可以被假设为均匀的。针对均匀介质，来自平坦表面换能器的压缩简档不应该随着介质的剪切模量而变化。因此，典型的剪切模量(比如，1kPa)可以用于这种建模。针对包含均匀材料或组织的对象，来自均匀假设的压缩简档可用于获得对象的2D弹性图像的一阶解决方案。然后，这种图像可用于运行另一种FEM模拟，以便获得更准确的压缩简档以用于更好地重建对象的真正2D弹性图像。

作为另一种用于估计压缩简档的示例方法，可以实现自适应估计方法。在此方法中，通过将空间平均值、加权平均值、中值或类似技术应用于沿深度轴(即，z轴)在一系列小空间窗口中所测得的总运动，可以对压缩简档进行估计。空间窗口可以是一维、二维或三维窗口。设想的是，剪切波将在深度方向上循环并且当应用求平均过程时将减小。因此，在求平均过程之后，设想压缩简档将保持。类似于以上所描述的曲线拟合方法，这些自适应方法可以结合来自多个空间位置和时间实例的信息，从而提高所估计压缩简档的准确度和精度。这可以通过应用具有专门内核(比如，高斯或拉普拉斯内核)的多维卷积技术，或其他多维滤波器(比如，中值或双边滤波器)来完成。

在一些情况中，不必估计压缩，而是压缩可以直接与剪切波解耦。如果跨换能器的全运动路径获取多个帧，从而使得运动在多个时间点(例如，帧)在深度方向(例如，z方向)以及横向方向(例如，x方向)两者上被获取，则传播剪切波信号可以从压缩中分离。这可以通过利用传播波和非传播运动的k空间表示的差异来完成。

如图9所示，对于在深度方向(例如，z轴)上传播的波，k空间k

因为压缩简档跨k

作为另一种用于估计压缩简档的示例方法，可以从k空间对压缩简档进行估计。在此方法中，在给定的帧中沿所观察信号u(z)的深度轴(例如，z轴)应用傅立叶变换，以便获得频域表示U(k

可以假设，压缩简档正在缓慢增加并且平滑。因此，k空间中的变形信号还将保持平滑。然而，剪切波的k空间分量将位于k空间中与图像空间中波长对应的不同点处。这导致了与这些剪切波分量的频率相关联的k

这种过程总体上在图4A中展示，其中，正弦剪切波信号已经被添加至使用弹性材料(黑色点)的FEM模拟所模拟的压缩简档。使用傅立叶变换来转变至k空间，所产生的频谱具有与50m

在一些实施例中，不对压缩简档进行估计，而是对数据获取进行更改以最小化换能器运动的效应。具体地，对于正弦振动的换能器，在正弦信号的波峰和波谷附近的时刻处进行的脉冲回波检测可以抑制由于换能器压缩引起的组织变形。

如图5所示，当在对称于正弦波波峰的时间t

类似地，在正弦波波谷处的检测(图5中的方形)还可以被实现为抑制由于换能器压缩引起的组织变形。这种方法对于脉冲回波事件之间的不同时间间隔有效，只要时间间隔被置于对称于正弦波的波峰或波谷。然而，因为所检测的运动是这两个脉冲回波事件之间的平均运动，所以时间间隔Δt＝t

在当换能器处于相同位置时不发生检测的情况下，随着换能器运动效应被最小化可以获得图像。这可以通过当检测关于运动简档的波峰和波谷对称时利用其他检测来恢复或估计运动来完成。这可以通过使用插值法、参数拟合或相移方法来执行。

作为另一种用于估计压缩简档的示例方法，可以通过利用耦合至或集成于超声换能器的运动传感器所获取的数据来估计压缩简档。在此方法中，用于测量加速度、速度、位移等的运动传感器可以耦合至超声换能器或集成在超声换能器内以测量其振动。这种方法提供了某些附加优势。作为一个示例，换能器的振动响应可能具有相较于用于驱动振动换能器的致动器的正弦信号的相位延迟。在这些实例中，测量换能器运动的传感器可以提供准确的同步以用于在如上述的正弦波的波峰或波谷处的超声检测。作为另一个示例，由运动传感器测量的换能器运动可用于在上述的运动减法方法中适当地缩放变形简档。因此，设想的是，使用运动传感器来测量换能器的运动可单独地或结合上述方法使用。

作为使用运动传感器的替代方案，由移动的换能器检测到的静止目标的运动还可用于估计换能器的位置，并且用于提高用于在如上述的正弦波的波峰或波谷处的超声检测的同步的准确度。举例来讲，静止目标可以是在换能器视场中在较深位置(在所述位置中，剪切波被完全衰减)处的非移动骨架或组织目标。

在本公开的一方面中，提供了用于对超声运动检测中时间延迟进行校正的系统和方法。在剪切波在组织中生成之后，可以使用脉冲回波运动检测来检测剪切波，如上所述。为了产生组织的机械特性的2D图像，期望对跨具有高时间分辨率的较大2D区域的组织运动进行同步检测。这可以通过“平面波”成像器来实现，在所述成像器中，来自2D检测区域内的每个像素的回波可以由平面超声波的单一传输进行重构。

然而，大多数商业超声扫描器不使用平面波成像，而是仍然使用顺序逐行扫描方法，其中，多个脉冲回波事件需要覆盖2D检测区域。因此，逐行扫描器具有比平面波成像器显著更低的成像帧速率。此外，利用逐行扫描器，需要考虑2D成像区域内每个成像行之间的时间延迟，以便根据所检测的剪切波正确地计算组织的机械特性。下面描述用于解决利用逐行扫描器检测剪切波的这种挑战的若干技术。

如图6所示，在逐行扫描器中实现的换能器12可以遵循某个顺序V

其中，V是成像矢量，每一个成像矢量包含在一个脉冲回波事件期间可以被平行地波束成形的n成像A行；N是每个成像区域内成像矢量的总数；并且M是在每个成像矢量位置处获取的脉冲回波事件的数量。图6中的实黑色方形表示脉冲回波事件的时刻，并且每个黑色方形附近的数字指示每个脉冲回波事件的时间序列。

假设脉冲回波事件的脉冲重复频率是PRF

经常期望保持较高的PRF

如图6所示，在已经针对前一个区域(例如，区域1)中的每个矢量收集了所有M样本之后，顺序跟踪进行至换能器12的另一个区域(例如，区域2)。重复相同的跟踪顺序，直到已经针对区域2中的每个矢量收集了所有M样本。针对区域2的检测事件还由区域2下的黑色实方形表示。可以重复这个过程，直到所有区域中的数据都被收集。

使用如图6所示的顺序跟踪方法所检测的运动信号存在两个问题。首先，在每个区域内检测事件时间上没有对准。第二，存在时间延迟，

在时间上相邻的区域之间(例如，区域1和区域2)。如以上所解释的，对组织机械特性的适当计算要求跨整个2D成像区域的运动信号应该具有相同的时间网格(即，应该是时间对准的)。

用于在每个区域内对准剪切波信号的方法，比如，在共同待决美国临时专利申请序列号62/072,167中所描述的时间插值法，可用于校正非对准运动信号。在时间插值法中，来自跨不同脉冲回波事件的每个矢量的回波首先用于计算组织运动，包括由于剪切波和换能器压缩造成的变形所引起的运动。因此，在时间插值法技术中，组织运动在黑色实方形处测量。在每个矢量位置(被指示为图6的区域1中的白色圆形)处的时间插值可以对准每个区域内的组织运动信号。

作为一个示例，通过将合适的相移应用于未对准的运动信号，运动信号可以被对准至公共时间网格。因为组织运动具有已知频率的正弦波，所以其允许附加方法来对准运动信号。假设由在区域1中的矢量V

M

其中，D

其中，D

这种相移方法可用于对准每个区域内的矢量，并且用于对准跨区域的矢量。

在图7中示出了如上所述的区域对区域时间对准方法的示例。在此示例中，50Hz的持续振动曾用于在仿组织均匀体模内产生持续剪切波。用于剪切波检测的PRF

作为对准运动信号的另一个示例，跨多个时间点(帧)在给定像素处所检测的组织运动可以拟合至正弦时间函数，以便根据时间估计正弦波信号的振幅和相位。一旦已知正弦时间函数的振幅和相位，则可以计算在任何时间点的运动信号。因此可以估计每个像素的正弦信号的振幅和相位参数，并且然后可以计算在公共对准的时间网格处所有像素的运动信号。这种方法可用于利用每个区域内或跨区域的不同矢量所检测的时间对准像素。

作为又另一个示例，运动信号可以基于由换能器的持续正弦振动产生的组织运动的循环性质来进行对准。在这些实例中，在每个空间位置处的运动信号在时间上以周期T重复。周期T是由持续振动的角频率ω确定的，

因此，正弦运动信号自身以相位差2π或时间周期T重复。因此，通过在图6中仔细设计时间延迟Δt可以时间上对准来自两个成像区域的运动信号，使得，

因此，通过根据等式(9)选择时间延迟，来自不同成像区域的运动信号被“自动地”对准，而不需要进一步时间对准。

基于等式(4和(9)，振动频率ω；检测PRF

替代性地，“等待时间”ε以添加至检测区域之间，使得

因此，等待时间ε可以被方便地调节以满足等式(9)的要求，或每个区域的检测时间可以由外部触发信号启动，所述外部触发信号在正弦振动信号的不同周期上被同步至固定相位。需注意，区域内对准对于去除由区域内顺序跟踪引入的时间延迟而言仍然是有必要的。

如在图7中可以看到的，在区域内对准之后，每个成像区域内的组织运动信号被对准并且准备好用于机械特性计算。然而，在一些实例中，可以在每个单独的区域内计算组织机械特性，而不是跨不同成像区域对准运动信号，并且然后可以将来自不同区域的估计组合成一个最终图。以此方式，运动信号跨区域的不连续性将无关紧要。如果在区域之间最终机械特性图中存在较小的不连续性，则空间滤波器(例如，中值滤波器)可用于使最终机械特性图平滑。这种方法适用于持续换能器振动剪切波，因为主要剪切波传播方向是从上到下远离换能器。这种传播方向主要是在具有很少跨区传播的区域内，并且因此计算每个可行区域内的组织机械特性。

用于运动信号对准问题的另一种解决方案是避免使用多个成像区域(即，仅使用单一成像区域)。使用单一成像区域一般需要大量的成像矢量(即，N的较大值)以便覆盖足够大的成像区域。然而，这种方法的挑战是如果N很大，则PRF

然而，对于由换能器振动产生的持续组织运动，可以通过利用用于驱动换能器的持续振动的正弦波的循环性质对混叠进行校正。如图8中所展示的，检测PRF

在图8中示出的示例中，PRF

在(例如，使用上述方法)对换能器运动的效应以及在超声脉冲回波运动检测中的时间延迟进行校正之后，在公共时间网格处的2D剪切波信号可用于使用标准弹性成像处理方法(比如，本振频率估计“LFE”、飞行时间、直接反演以及其他方法)来计算组织的机械特性。

在以上给出的示例中，假设超声换能器以单一频率进行振动。在此描述的方法还可以轻易地扩展至持续换能器振动包含多个正弦频率或啁啾信号的情况。如此，可以在多个频率上测量机械特性。

已经针对使用1D线性超声阵列换能器的2D弹性成像描述了以上技术。这些技术还可适用于单元素、1D弯曲阵列、1.5D阵列、1.75D阵列和2D阵列换能器。对于单元素换能器，所述方法可以向下缩放至1D弹性成像。对于2D阵列，所述方法可以扩展至3D弹性成像。还可以将对换能器运动的校正与超声运动检测中的延迟组合在一起。在此描述的方法还可用于测量组织和非组织材料(例如，聚合物)的机械特性。

已经根据一个或多个优选实施例对本发明进行了描述，并且应当理解的是，除明确陈述的以外，许多等价物、备选方案、变化和修改都是可能的并且在本发明的范围之内。