一种基于多频剪切波幅值分析的超声弹性成像技术的制作方法

本发明涉及医学成像领域和人机接口领域，尤其涉及一种基于多频剪切波幅值分析的超声弹性成像技术。

背景技术：

超声弹性成像为近年来新兴的组织成像技术，其可以用来获得人体组织的弹性信息，且人体的正常组织和病变组织表现为不同的弹性，因而该技术在医疗诊断领域中得到了广泛的研究和应用。

另一方面，人机接口能够采集人体的生理信号，通过处理能够获得人体的运动意图，进一步可用于控制外部设备。而骨骼肌等组织在不同的运动状态下也具有不同的弹性系数。因此，超声弹性成像技术可以为人机接口提供可以用于运动意图识别的组织弹性信息，该技术在人机接口领域有着巨大的潜在应用。

无论在医学成像领域，还是在人机接口领域，对生物组织弹性模量测量精度的提升一直是技术研究和产业化的热门。特别地，人机接口的使用人群包括可穿戴设备爱好者及截肢患者，该人群在使用过程中通常对设备的便携性有着较高的要求。具有高精度且可易于被简单设备实现的组织弹性测量方法将有着更加广泛的市场需求。

由于软生物组织中含有固体分量(如细胞、肌纤维等)，人体软组织既可以传递纵波，也可以传递低速的剪切波。现有超声弹性成像技术的分支包括对组织的杨氏模量进行成像和对组织的剪切模量进行成像。因为组织中传递的剪切波波速和组织的剪切模量有关(ct为剪切波波速，μ为剪切模量，ρ为组织密度)，现针对剪切模量的测量多是根据剪切波波速测量实现的，并且基于该方式的剪切波弹性成像技术已经实现了商业化。

然而，测量剪切波波速的实现过程中存在着诸多问题。比如，测量剪切波首先需要对剪切波进行追踪，追踪方式大多为先获取剪切波的相位等信息，根据剪切波相位，计算剪切波传递的时间，结合剪切波传递的距离，计算剪切波的波速。获取剪切波相位信息的准确度受超声信号的脉冲重复率、采样率和噪声等影响。除此之外，剪切波波速的计算需要剪切波传递一定的距离，该过程对剪切波弹性成像的方式有很强的限制，探头声场需要足够大，所覆盖剪切波传递的范围足以用于计算其波速。换而言之，该种测量方式只能通过线阵等多振元超声(b超)探头实现，在不依靠移动平台的情况下，该种测量方式无法通过单(双)振元探头实现。基于剪切波波速测量的弹性模量计算方法，导致所需的设备通常价格昂贵、复杂且沉重。因此，在获取生物组织的剪切模量过程中，采用测量剪切波波速的方式不仅限制其测量精度，而且限制剪切模量测量的应用场景。

因此，本领域的技术人员致力于开发一种测量组织剪切模量的新方法，以改善剪切模量的测量精度，同时降低剪切模量测量过程中所需设备的复杂度，推进人体组织弹性测量设备的便携化和穿戴化实现。

技术实现要素：

有鉴于现有技术的上述缺陷，本发明所要解决的技术问题是提高测量生物组织的剪切模量的精度和抗干扰能力，降低设备复杂度。

为实现上述目的，本发明提供了一种基于多频剪切波幅值分析的超声弹性成像技术，包括以下步骤：

步骤1、对剪切波的激励脉冲进行调制，使其产生不同的角频率的所述剪切波；

步骤2、测量不同所述角频率的所述剪切波所引起的组织中质点振动幅值；

步骤3、对在不同所述角频率下所述剪切波的所述质点振动幅值进行对比，得到最大所述质点振动幅值所对应的所述剪切波的所述角频率；

步骤4、根据所述角频率获得所述剪切波的波速或所述组织的剪切模量或所述剪切波的波速和所述组织的剪切模量。

进一步地，所述步骤1的所述剪切波的激励脉冲包括连续激励脉冲。

进一步地，所述步骤1包括软件方式和硬件方式调制。

进一步地，所述步骤1包括声辐射力调制方式和扫频方式调制。

进一步地，所述扫频方式包括线性扫频和对数扫频。

进一步地，所述步骤1所采用的探头包括声辐射力激励振元和超声回波检测振元。

进一步地，所述声辐射力激励振元和所述超声回波检测振元并列放置。

进一步地，所述声辐射力激励振元和所述超声回波检测振元同轴放置。

进一步地，所述步骤3建立所述剪切波的波速和最大所述质点振动幅值所对应所述角频率的关系模型。

进一步地，所述步骤3建立所述组织的剪切模量和最大所述质点振动幅值所对应所述角频率的关系模型。

本发明与现有技术相比较，具有如下显而易见的实质性特点和显著优点：

本发明通过对比不同频率剪切波下的质点振动幅值，获取最大振动幅值所对应的最大角频率；通过剪切波最大幅值所对应的激励频率和利用该频率与剪切模量或者剪切波波速间的关系模型，获取剪切模量或者剪切波波速。本发明所提出方法无需对幅值进行绝对值的测量，而是进行相对值的测量，幅值信息相对于现有技术中常用的相位信息在被测量时具有更高的抗干扰性。所提出方法可通过两个超声振元实现，易于实现便携化，在人机借口领域有着极大的应用潜力。

以下将结合附图对本发明的构思、具体结构及产生的技术效果作进一步说明，以充分地了解本发明的目的、特征和效果。

附图说明

图1是本发明的一个较佳实施例的方法流程图；

图2是本发明的一个较佳实施例的特定剪切模量或剪切波速下，质点振动幅值a随频率ω的变化；

图3是本发明的一个较佳实施例的不同剪切波波速下的最大幅值a所对应的角频率ωm；

图4是本发明的一个较佳实施例的声辐射力激励振元和超声回波检测振元并列放置；

图5是本发明的一个较佳实施例的声辐射力激励振元和超声回波检测振元同轴放置。

具体实施方式

以下参考说明书附图介绍本发明的多个优选实施例，使其技术内容更加清楚和便于理解。本发明可以通过许多不同形式的实施例来得以体现，本发明的保护范围并非仅限于文中提到的实施例。

在附图中，结构相同的部件以相同数字标号表示，各处结构或功能相似的组件以相似数字标号表示。附图所示的每一组件的尺寸和厚度是任意示出的，本发明并没有限定每个组件的尺寸和厚度。为了使图示更清晰，附图中有些地方适当夸大了部件的厚度。

图1是本发明的一个较佳实施例的方法流程图。本发明所提出的方法可分为四步进行实现：

1.本发明所提出的方法需要对剪切波的激励脉冲包括但不限于连续激励脉冲进行调制，使其产生不同的角频率ω剪切波；

2.测量不同频率剪切波所引起的组织中质点的振动幅值a；

3.对在不同角频率ω下剪切波中的质点振动幅值进行对比，得到最大质点振动幅值am所对应的剪切波角频率ωm；

4.根据该角频率ωm获得剪切波的波速ct或组织的剪切模量μ或剪切波的波速ct和组织的剪切模量μ。

本发明的理论原理基于剪切波的产生与传递方程。以圆形活塞辐射声场为例，在该声场中，剪切波所引起的质点振动位移sx如下式所示：

其中，α为声衰减系数，a圆形活塞换能器的半径，c为纵波声速，ρ为组织密度，f为声场描述，i为声场的发射声强，j0为零阶的一类贝塞尔函数，v为组织粘度，x为质点到活塞换能器发射平面的距离，r为质点到平面活塞换能器轴线的距离，剪切波激励信号的包络，可对声辐射力进行调制。

当调制声辐射力激励信号以达到时，剪切波所引起的质点振动位移sx如下式所示：

在实际测量过程中，选定测量所用换能器和换能器的激励强度后，并主要考虑待测组织中剪切模量μ或横波波速ct的变化，特定位置质点由剪切波所引起的运动将由剪切模量μ或横波波速ct及角频率ω决定。其中，角频率ω可以人为控制。

因此，可以通过软件或硬件等方式，对激励波形(声辐射力)进行调制，产生不同角频率ω的剪切波。

特定参数下，检测特定剪切模量μ或横波波速ct下，质点振动位移sx的归一化幅值a与角频率ω的变化关系如图2所示。

可根据图2进一步获得最大幅值am所对应的角频率ωm，在所选定的参数下，不同剪切波波速ct所对应的频率ωm如图3所示。可见，剪切波波速ct和频率ωm有着很强的关系，而在本发明所选定的参数下表现为线性。

根据剪切波波速ct和频率ωm的关系模型，可以由已知量频率ωm获得剪切波波速ct或剪切模量μ。

在探头布置方面，上述方式即可以通过将声辐射力激励振元5和超声回波检测振元6进行并列放置，如图4，也可以将声辐射力激励振元5和超声回波检测振元6进行同轴放置，如图5。示意图只是提供振元放置方式的基本解释说明，并不将振元数目限制为2个，即使本发明所提出的方法可以仅通过2个振元实现。

本发明提出了一种新的检测组织剪切弹性的方式，与现有技术检测相位信息进行剪切波的波速计算不同，本发明是检测剪切波在不同激励条件下的幅值变化，确定最大幅值所对应的激励频率，并根据该激励频率与剪切模量或者剪切波波速间的关系模型获取剪切模量或剪切波波速。进一步地，本发明所提出方法并非是检测剪切波所引起质点位移的绝对值，而是在不同激励频率下变化的相对值，相对检测剪切波的相位信息而言，所提出方法具有更强的抗干扰能力。

本发明通过剪切波最大幅值所对应的激励频率和利用该频率与剪切模量或者剪切波波速间的关系模型，获取剪切模量或者剪切波波速。

所提出方法不限于上述的声辐射力调制方式，还包括但不限于频率的扫频方式，如线性扫频和对数扫频。

所提出方法不限于应用于上述的探头或振元布置方式，包括但不限于不同的振元数目，振元的谐振频率和振元大小。

所提出方法中不仅限于建立剪切波波速和最大幅值所对应频率的关系模型，还包括剪切模量和最大幅值所对应频率的关系模型。因为剪切波波速与剪切模量在一定范围内等效。

以上详细描述了本发明的较佳具体实施例。应当理解，本领域的普通技术无需创造性劳动就可以根据本发明的构思作出诸多修改和变化。因此，凡本技术领域中技术人员依本发明的构思在现有技术的基础上通过逻辑分析、推理或者有限的实验可以得到的技术方案，皆应在由权利要求书所确定的保护范围内。