一种基于组织中微泡动力学模型的微弹性成像方法

摘要

一种基于组织中微泡动力学模型的微弹性成像方法，根据微泡的振动情况会受到周围组织特性的影响，组织中的微泡动力学模型可建立微泡振动信号与组织弹性的关系，建立一种与微泡信号相关性强而与组织信号相关性弱的母小波，通过脉冲逆转与小波变换相结合的成像算法对微泡信号进行检测成像，当检测信号与模型信号最接近时，可得到最大组织造影比，从而反推回微泡周围十几至几十微米组织的弹性参数；该方法可应用于高强度聚焦超声治疗过程实时监控及生物薄层组织弹性检测，可克服一般弹性成像需要外部施压、易受边界条件影响等局限，有效将成像分辨率从毫米级提高至微米级。

说明书

技术领域

本发明属于生物医学超声成像技术领域，特别涉及到一种基于组 织中微泡动力学模型的微泡振动信号检测、微泡周围组织参数反求及 微弹性成像的系统与方法。

背景技术

弹性成像主要是根据不同组织弹性系数的差异，在受到应力作用 后组织应变程度的不同来显示及定位病变引起的弹性改变，对乳腺结 节、甲状腺瘤、肝硬化等疾病的诊断以及微波、射频、高强度聚焦超 声等热消融引起的组织固化监控有着重要的价值。超声弹性成像是目 前最主要的弹性成像方法，同时利用核磁共振方法也可进行弹性成 像。超声弹性成像主要包括：（1）静态/准静态压缩弹性成像，（2） 声辐射力动态压缩弹性成像，（3）依靠自身血压变化作为激励的心脏 /血管弹性成像。现有弹性成像方法的局限性主要包括：（1）静态压 缩弹性成像中，组织内部的形变主要由外力作用的边界条件与实际弹 性分布决定，外力边界条件的改变会造成应变估计结果与实际弹性分 布不相符，导致明显的应变估计误差。（2）动态压缩弹性成像中，可 通过局部施加激励来计算小范围的弹性分布，但是成像分辨率会受激 励源的尺寸与成像速度的限制。（3）现有弹性成像技术的分辨率基本 在毫米级别，暂时无法满足区分体内尺寸较小的薄层组织，如血管壁、 上皮组织、胃粘膜等弹性变化的需求。

因此，本发明在研究弹性成像、微泡动力学及造影成像种不同理 论方法与实验的基础上，提出了一种将组织中微泡动力学模型与组织 中微泡信号检测相结合的微弹性成像方法。该方法可通过微泡振动信 号反推微泡周围组织弹性参数，克服现有弹性成像需要外部施压、易 受边界条件影响等局限，有效将成像分辨率从毫米级提高至微米级。

发明内容

针对现有超声弹性成像的现状，本发明的目的在于提出一种基于 组织中微泡动力学模型的微弹性成像方法，以组织中的微泡动力学模 型为基础，采用基于模型的微泡成像及信号处理方法，利用微泡采集 信息与模型理论信号的相关性反推振动微泡周围的组织弹性参数；该 方法可与治疗超声空化物理结合，研究和发展超声治疗过程实时微弹 性监控成像技术；与相变微泡超声导入和相变控制结合，研究和发展 体内超薄层组织微弹性成像技术。

为了达到上述目的，本发明的技术方案为：

一种基于组织中微泡动力学模型的微弹性成像方法，包括如下步 骤：

（1）将设置在探测目标上的线阵换能器通过具有射频数据采集接 口的可编程超声设备与主控计算机连接；

（2）选择能反应组织弹性参数的微泡动力学模型，通过实验手段 对模型进行检验，并对其进行灵敏度分析和可靠性评估，最终选定对 组织弹性参数敏感，可信度高的模型构建母小波，对其进行求解，得 到微泡振动半径随时间的变化曲线及径向振动速度和加速度，根据体 积随时间的变化率计算振动微泡辐射出的声压力：其中P为振动微泡辐射出的声压，ρ为液体密度，r为微泡信号检测 换能器与微泡之间的距离，R和分别微泡半径及径向加速度；

（3）母小波参数确定：在计算作用于微泡的激励声压时，需考虑 不同组织的声衰减及微泡所在位置深度，利用公式对到达 不同位置的微泡声压进行衰减估计，其中z为声传播的距离，μ₀为幅 度衰减系数，P₀为输出激励声压；在需要带入组织弹性参数时，根据 正常组织弹性参数经验值及普通弹性成像方法实验值确定组织弹性 参数范围，并根据计算量和计算时间确定步长，再将组织弹性参数及 其对应的步长代入母小波；

（4）利用可编程线阵换能器发射脉冲逆转序列，即沿一条扫描线 相继发射和接受0相与180相的同幅度脉冲，后续扫描线重复操作直 至完成二维平面扫查，形成单帧射频数据；

（5）基于步骤（2）、（3）所得到的母小波，分别对0相和180相 回波信号进行小波相关性分析，得到一系列小波相关系数，分别选取 最大的小波相关系数所在尺度的所有小波系数去代替原始的射频回 波信号；

（6）将0相和180想回波信号中表征原始射频回波信号的小波相 关系数进行相加得到一条波束并进行后续处理成像，包括对数压缩、 坐标变换和灰度转换，并进行图像显示；

（7）使用造影组织比（CTR）对不同组织弹性参数条件下的图像 进行评价，选取相同面积同一扫描深度的感兴趣区域计算其对应的 CTR值：当CTR值最大时，认为相对应的组织弹性参数即为或最接近 感兴趣区域的平均组织弹性参数；

（8）对成像区域继续进行分区，分区内含有有限个可分辨微泡， 以步骤（7）得到的组织弹性参数为中心以尽可能小的步长上下移动 得到对应母小波，对每一分区继续进行信号检测成像，通过遗传算法 进行相关反馈，计算图像CTR值，即可得到微泡周围组织的弹性值。

所述的组织弹性参数范围，其提取参数的操作步骤如下：

（1）通过对相变微泡浓度的精确控制，使得单个相变微泡可以 缓慢流过超声探头聚焦区域，使得所得回波信号为单个微泡背向散射 信号；通过控制发射声波的相位、输入功率和输入时间，控制相变微 泡在不相变的条件下导入薄层组织内部；

（2）延长辐照时间或加大辐照功率，使得目标区域温度达到相 变微泡的相变阈值，在薄层组织内部产生稀疏微泡，利用可编程超声 设备发射微泡振动激励信号，使微泡产生振动；

（3）通过微泡信号检测换能器采集微泡回波信号，并利用基于 组织中微泡动力学模型的微弹性成像方法对治疗区域进行微弹性成 像及微泡周围组织弹性估计。

本发明与现有技术相比，具有下列优点：

(1）与一般弹性成像方法比较，本发明提出的基于微泡动力学模 型的造影成像技术，检测灵敏度更高，可有效抑制组织信号，提高造 影组织比。

（2）本发明提出的基于组织中微泡动力学模型的微弹性成像方 法，可克服一般弹性成像方法需要外加施压、易受边界条件影响等局 限。

（3）本发明利用微泡信号反求周围微米级别组织弹性特性，有 效将成像分辨率从毫米级别提高至了微米级别。

（4）本发明提出的生物薄层组织弹性实时监控成像方法，提高 了研究生物薄层组织弹性的实验方法，除了具有上述优点外，还为生 物薄层组织病变显示及定位提供了手段。

附图说明

图1是本发明提出的基于组织中微泡动力学的造影成像方法的 方法流程图。

图2是组织弹性参数方法操作流程图。

图3是成像区域分区及参数反求。

图4是相变微泡相变过程示意图。

具体实施方式

参照图1，一种基于组织中微泡动力学模型的微弹性成像方法， 包括如下步骤：

（1）将设置在探测目标上的线阵换能器通过具有射频数据采集接 口的可编程超声设备与主控计算机连接；

（2）选择能反应组织弹性参数的微泡动力学模型，通过实验手段 对模型进行检验，并对其进行灵敏度分析和可靠性评估，最终选定对 组织弹性参数敏感，可信度高的模型构建母小波，对其进行求解，得 到微泡振动半径随时间的变化曲线及径向振动速度和加速度，根据体 积随时间的变化率计算振动微泡辐射出的声压力：其中P为振动微泡辐射出的声压，ρ为液体密度，r为微泡信号检测 换能器与微泡之间的距离，R和分别微泡半径及径向加速度；

（3）母小波参数确定：在计算作用于微泡的激励声压时，需考虑 不同组织的声衰减及微泡所在位置深度，利用公式对到达 不同位置的微泡声压进行衰减估计，其中z为声传播的距离，μ₀为幅 度衰减系数，P₀为输出激励声压；在需要带入组织弹性参数时，根据 正常组织弹性参数经验值及普通弹性成像方法实验值确定组织弹性 参数范围，并根据计算量和计算时间确定步长，再将组织弹性参数及 其对应的步长代入母小波；

（4）利用可编程线阵换能器发射脉冲逆转序列，即沿一条扫描线 相继发射和接受0相与180相的同幅度脉冲，后续扫描线重复操作直 至完成二维平面扫查，形成单帧射频数据；

（5）基于步骤（2）、（3）所得到的母小波，分别对0相和180相 回波信号进行小波相关性分析，得到一系列小波相关系数，分别选取 最大的小波相关系数所在尺度的所有小波系数去代替原始的射频回 波信号；

（6）将0相和180想回波信号中表征原始射频回波信号的小波相 关系数进行相加得到一条波束并进行后续处理成像，包括对数压缩、 坐标变换和灰度转换，并进行图像显示；

（7）使用造影组织比（CTR）对不同组织弹性参数条件下的图像 进行评价，选取相同面积同一扫描深度的感兴趣区域计算其对应的 CTR值：当CTR值最大时，认为相对应的组织弹性参数即为或最接近 感兴趣区域的平均组织弹性参数。

（8）对成像区域继续进行分区，参照图3，分区内含有有限个可 分辨微泡，以步骤（7）得到的组织弹性参数为中心以尽可能小的步 长上下移动得到对应母小波，对每一分区继续进行信号检测成像，通 过遗传算法进行相关反馈，计算图像CTR值，即可得到微泡周围组织 的弹性值。

参照图2，所述的组织弹性参数范围，其提取参数的操作步骤如 下：

（1）通过对相变微泡浓度的精确控制，使得单个相变微泡可以 缓慢流过超声探头聚焦区域，使得所得回波信号为单个微泡背向散射 信号，通过控制发射声波的相位、输入功率和输入时间，控制相变微 泡在不相变的条件下导入薄层组织内部。

（2）延长辐照时间或加大辐照功率，使得目标区域温度达到相 变微泡的相变阈值，在薄层组织内部产生稀疏微泡，利用可编程超声 设备发射微泡振动激励信号，使微泡产生振动，参照图4；

（3）通过微泡信号检测换能器采集微泡回波信号，并利用基于 组织中微泡动力学模型的微弹性成像方法对治疗区域进行微弹性成 像及微泡周围组织弹性估计，参照图5。

本发明实施时需要成像系统，该系统由微泡产生子系统和微弹性 成像子系统构成：

所述的微泡产生子系统包括双通道任意波形发生器、功率放大 器、、可编程超声设备、微泡控制换能器；其中双通道任意波形发生 器用以实现整个系统的时序控制；可编程超声设备用以控制华能器输 出声压、能量、频率等参数。通过换能器的不同激励，实现相变微泡 的组织导入及控制相变，在生物薄层组织中产生空化微泡。

微弹性成像子系统包括微泡信号检测换能器、宽带增益接收器、 高速数据采集卡和主控计算机；微泡信号检测换能器采集得到的数据 经过宽带增益接收器放大经高速数据采集卡输出到主控计算机，并用 数据处理软件对数据进行实时处理。